[Paraview] Quadrature in Paraview
Burlen Loring
bloring at lbl.gov
Fri Mar 29 12:43:45 EDT 2013
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.paraview.org/pipermail/paraview/attachments/20130329/cdeed8f1/attachment-0001.htm>
More information about the ParaView
mailing list