Read Xdmf

From XdmfWeb
Revision as of 13:59, 20 May 2016 by Crayzeewulf (talk | contribs) (→‎Reading XDMF Data: Placed code in source tags so it get syntax highlighting. Cleaned up indentation of the code.)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigationJump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

Reading XDMF Data

TwoHex.jpg

The following Xdmf XML file is a simple example of a Uniform Grid that contains two Hexahedron that share a face. There are values centered at the nodes and at the cell centers. The values for geometry, connectivity, and scalars are all stored directly in the XML file.

<source lang="xml" line="1"> <?xml version="1.0" ?> <!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []> <Xdmf>

   <Domain>
       <Grid Name="TestGrid">
           <Topology Type="Hexahedron" NumberOfElements="2" >
               <DataItem Format="XML" DataType="Float"
                         Dimensions="2 8">
                   0 1 7 6   3 4 10 9
                   1 2 8 7   4 5 11 10
               </DataItem>
           </Topology>
           <Geometry Type="XYZ">
               <DataItem Format="XML" DataType="Float" Precision="8"
                         Dimensions="4 3 3">
                   0.0   0.0   1.0
                   1.0   0.0   1.0
                   3.0   0.0   2.0
                   
                   0.0   1.0   1.0
                   1.0   1.0   1.0
                   3.0   2.0   2.0
                   
                   0.0   0.0   -1.0
                   1.0   0.0   -1.0
                   3.0   0.0   -2.0
                   
                   0.0   1.0   -1.0
                   1.0   1.0   -1.0
                   3.0   2.0   -2.0
               </DataItem>
           </Geometry>
           <Attribute Name="NodeValues" Center="Node">
               <DataItem Format="XML" DataType="Float" Precision="8"
                         Dimensions="4 3" >
                   100 200 300
                   300 400 500
                   300 400 500
                   500 600 700
               </DataItem>
           </Attribute>
           <Attribute Name="CellValues" Center="Cell">
               <DataItem Format="XML" DataType="Float" Precision="8"
                         Dimensions="2" >
                   100 200
               </DataItem>
           </Attribute>
       </Grid>
   </Domain>

</Xdmf> </source>

The XML is stored in the file "MyGrid.xmf". The following Python program demonstrated parsing the XML and retrieving the data values.

<source lang="python" line="1">

  1. ! /usr/bin/env python

from Xdmf import * reader = XdmfReader.New() dom = XdmfReader.read('MyGrid.xmf')

  1. We now have a tree. Find the one and only Grid element

grid = dom.getUnstruturedGrid(0)

top = grid.GetTopology() print 'Values = ', top.getValuesString()

  1. Release values from data when done

top.release()

geo = grid.GetGeometry() print 'Geo Type = ', geo.getType().getName(), ' # Points = ', geo.getNumberPoints() print 'Points = ', geo.getValuesString() geo.release()

  1. for each Attribute, print Light Data and Values

for i in range(grid.getNumberAttributes()):

   attr = grid.getAttribute(i)
   print 'Attribute ', i, ' Name = ', attr.getName()
   # Attribute HeavyData is not Updated by default
   # there could potentially be many causing huge IO
   attr.read()
   print 'Values ', attr.getValuesString()
   attr.release()

</source>