[Paraview] Quadrature in Paraview
Burlen Loring
bloring at lbl.gov
Fri Mar 29 12:56:35 EDT 2013
That's good to know, however, you should follow my instructions for now
since you have to apply the patch I sent.
On 03/29/2013 09:49 AM, Andy Bauer wrote:
> 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
> <mailto: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
> <http://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
> <http://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 <tel:4%205%206%207%208%209%2010%2011>
>> </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 bywww.kitware.com <http://www.kitware.com>
>>
>> Visit other Kitware open-source projects athttp://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 <http://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/dd8d8f24/attachment-0001.htm>
More information about the ParaView
mailing list