[Paraview] Quadrature in Paraview
Andy Bauer
andy.bauer at kitware.com
Fri Mar 29 12:49:27 EDT 2013
Hi,
FYI: We now generate installers for ParaView nightly for git HEAD. They're
available at http://www.paraview.org/paraview/resources/software.php and by
choosing Nightly for Version of ParaView.
Andy
On Fri, Mar 29, 2013 at 12:43 PM, Burlen Loring <bloring at lbl.gov> wrote:
> Hi afad sfgsfdgsfd,
>
> Fixes for bugs introduced by changes in VTK's xml file format are in the
> pipeline, so you'll have to checkout the latest git HEAD and build ParaView
> from source to get the bug fixes. If that's something you're willing to do,
> here are first steps
>
> # get sources
> $ git clone --recursive git://www.paraview.org/ParaView.git
> # apply the patch
> # patch is attached and I just now submitted to PV, will show up in PV
> after the next review.
> $ cd ParaView
> $ git am -3 0001-QuadraturePoints-any-integer-type-for-offsets.patch
> # build pv
>
> That will get you the bug fixes. Build it, use default build opts is fine.
>
> Then if you like you can get the test data, and open the file
> ParaViewData/Data/blow.vtk, construct the pipeline showed in the post (
> http://markmail.org/message/6qz3ebszvny4vwxp#query:+page:1+mid:6qz3ebszvny4vwxp+state:results).
> That shows how it's supposed to work.
>
> $ git clone www.paraview.org/ParaViewData.git
>
> Now onto your data, the xml you gave is incorrect , it's missing the
> offset array. let me see if I can get your example working, I'll post a
> follow up, hopefully in a few minutes.
>
> Burlen
>
>
> On 03/28/2013 12:43 AM, afad sfgsfdgsfd wrote:
>
> I am trying to understand quadrature in paraview. I am running finite
> element simulations in which I have displacements at the nodes of the
> element, and stress/strain at the quadrature points. I have linear and
> quadratic hexahedral or tetrahedral elements. I have been spending days
> trying to figure out how to implement data on the quadrature points. I
> have downloaded just about every document I can find.
>
> First: I downloaded CylinderQuadratic.vtk as shown in
> VTK-Quadrature-Point-Design-Doc. I was able to generate the quadrature
> scheme dictioanary and generate quadrature points. But, when I tried to
> use the Interpolate to Quadrature Points filter I could not select
> QuadratureOffset in the Select Source Array drop down menu. I was able to
> do this when using the Generate Quadrature Points filter.
>
> Second: Considering I couldn't figure out how to interpolate quadrature
> points using CylinderQuadratic.vtk, I decided to investigate just using a
> single element. I took quadraticTetra01.vtu and modified it to only have
> on quadratic tetrahedron. Then I used paraview to generate the quadrature
> scheme dictionary. I saved this into a VTU file in order to see how to
> implement the dictionary on my own. When I tried to reopen this new VTU
> file, paraview was unable to open the file. Specifically, it choked on the
> QuadratureOffset data array.
> This is the file that paraview cannot open
>
> <VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
> <UnstructuredGrid>
> <Piece NumberOfPoints="10" NumberOfCells="1">
> <PointData Scalars="scalars">
> <DataArray type="Float32" Name="scalars" format="ascii"
> RangeMin="0" RangeMax="1">
> 1 1 1 1 0 0
> 0 0 0 0
> </DataArray>
> </PointData>
> <CellData>
> <DataArray type="Int64" Name="QuadratureOffset" format="ascii"
> RangeMin="0" RangeMax="0">
> <InformationKey name="DICTIONARY"
> location="vtkQuadratureSchemeDefinition">
> <vtkQuadratureSchemeDefinition>
> <CellType value="24">
>
> </CellType>
> <NumberOfNodes value="10">
>
> </NumberOfNodes>
> <NumberOfQuadraturePoints value="4">
>
> </NumberOfQuadraturePoints>
> <ShapeFunctionWeights>
> 1.5625000000000000e-001 -9.3750000000000000e-002
> -9.3750000000000000e-002 -9.3750000000000000e-002
> 3.1250000000000000e-001 6.2500000000000000e-002
> 3.1250000000000000e-001 3.1250000000000000e-001
> 6.2500000000000000e-002 6.2500000000000000e-002
> -9.3750000000000000e-002 7.0312500000000000e-002
> -1.1718750000000000e-001 -9.3750000000000000e-002
> 2.8125000000000000e-001 4.2187500000000000e-001
> 9.3750000000000000e-002 6.2500000000000000e-002
> 2.8125000000000000e-001 9.3750000000000000e-002
> -9.3750000000000000e-002 -1.1718750000000000e-001
> 7.0312500000000000e-002 -9.3750000000000000e-002
> 9.3750000000000000e-002 4.2187500000000000e-001
> 2.8125000000000000e-001 6.2500000000000000e-002
> 9.3750000000000000e-002 2.8125000000000000e-001
> -9.3750000000000000e-002 -5.4687500000000000e-002
> -5.4687500000000000e-002 3.7500000000000000e-001
> 3.1250000000000000e-002 1.5625000000000000e-002
> 3.1250000000000000e-002 3.7500000000000000e-001
> 1.8750000000000000e-001 1.8750000000000000e-001
> </ShapeFunctionWeights>
> <QuadratureWeights>
> 2.1219962781139991e-314 2.1219957914593380e-314
> 1.1013158155129526e-311 1.1501217187461613e-311
> </QuadratureWeights>
>
> </vtkQuadratureSchemeDefinition>
>
> </InformationKey>
> 0
> </DataArray>
> </CellData>
> <Points>
> <DataArray type="Float32" Name="Points" NumberOfComponents="3"
> format="ascii" RangeMin="0" RangeMax="1.1874342107">
> 0 0 0 1 0 0
> 0.5 0.80000001192 0 0.5 0.40000000596 1
> 0.5 0 -0.20000000298 0.60000002384 0.60000002384 0
> 0.30000001192 0.40000000596 0 0.40000000596 0.20000000298 0.5
> 0.85000002384 0.30000001192 0.5 0.5 0.60000002384 0.44999998808
> </DataArray>
> </Points>
> <Cells>
> <DataArray type="Int64" Name="connectivity" format="ascii"
> RangeMin="0" RangeMax="9">
> 0 1 2 3 4 5
> 6 7 8 9
> </DataArray>
> <DataArray type="Int64" Name="offsets" format="ascii"
> RangeMin="10" RangeMax="10">
> 10
> </DataArray>
> <DataArray type="UInt8" Name="types" format="ascii" RangeMin="24"
> RangeMax="24">
> 24
> </DataArray>
> </Cells>
> </Piece>
> </UnstructuredGrid>
> </VTKFile>
>
>
> Third: I investigated quadrature for hexahedron. I made a similar file
> to quadraticTetra01.vtu and tried to generate a quadrature dictionary.
> Upon hitting apply for this, paraview crashes. It crashes for both linear
> and quadratic hexahedron. This is the file it crashes on,
>
> <?xml version="1.0"?>
> <VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian"
> compressor="vtkZLibDataCompressor">
> <UnstructuredGrid>
> <Piece NumberOfPoints="12" NumberOfCells="2">
> <PointData Scalars="scalars">
> <DataArray type="Float32" Name="scalars" format="ascii">
> 1 1 1 1 0 0 0 0 0 0 1 1
> </DataArray>
> </PointData>
>
> <Points>
> <DataArray type="Float32" NumberOfComponents="3" format="ascii">
> -1.0 -1.0 1.0
> -1.0 -1.0 -1.0
> -1.0 1.0 -1.0
> -1.0 1.0 1.0
> 1.0 -1.0 1.0
> 1.0 -1.0 -1.0
> 1.0 1.0 -1.0
> 1.0 1.0 1.0
> 3.0 -1.0 1.0
> 3.0 -1.0 -1.0
> 3.0 1.0 -1.0
> 3.0 1.0 1.0
> </DataArray>
> </Points>
> <Cells>
> <DataArray type="Int32" Name="connectivity" format="ascii">
> 0 1 2 3 4 5 6 7
> 4 5 6 7 8 9 10 11
> </DataArray>
> <DataArray type="Int32" Name="offsets" format="ascii">
> 8 16
> </DataArray>
> <DataArray type="UInt8" Name="types" format="ascii">
> 12 12
> </DataArray>
> </Cells>
> </Piece>
> </UnstructuredGrid>
> </VTKFile>
>
>
>
>
> My question is can anyone give me a quick example to implement
> quadrature in paraview. I am at wits end trying to figure this out.
>
>
>
>
> _______________________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the ParaView Wiki at: http://paraview.org/Wiki/ParaView
>
> Follow this link to subscribe/unsubscribe:http://www.paraview.org/mailman/listinfo/paraview
>
>
>
> _______________________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the ParaView Wiki at:
> http://paraview.org/Wiki/ParaView
>
> Follow this link to subscribe/unsubscribe:
> http://www.paraview.org/mailman/listinfo/paraview
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.paraview.org/pipermail/paraview/attachments/20130329/4d296bf0/attachment.htm>
More information about the ParaView
mailing list