1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155
|
#!/bin/env python
from __future__ import print_function
from Xdmf import *
txt = """<?xml version="1.0" ?>
<!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []>
<Xdmf xmlns:xi="http://www.w3.org/2001/XInclude" Version="2.0">
<Domain>
<Grid Name="Polyvertex" GridType="Uniform">
<Topology TopologyType="Polyvertex"
Dimensions="3"
NodesPerElement="3">
<DataItem Format="XML"
Dimensions="3 3"
NumberType="Float">
0 1 2
4 5 6
8 9 10
</DataItem>
</Topology>
<Geometry Type="XYZ">
<DataItem Format="XML" Dimensions="4 4 3">
0.0 0.0 0.0
1.0 0.0 0.0
1.0 1.0 0.0
0.0 1.0 0.0
0.0 0.0 2.0
1.0 0.0 2.0
1.0 1.0 2.0
0.0 1.0 2.0
0.0 0.0 3.0
1.0 0.0 3.0
1.0 1.0 3.0
0.0 1.0 3.0
0.0 0.0 4.0
1.0 0.0 4.0
1.0 1.0 4.0
0.0 1.0 4.0
</DataItem>
</Geometry>
<Attribute Name="Node Centered Values" Center="Node">
<DataItem Format="XML" Dimensions="4 4">
100 200 300 400
500 600 600 700
800 900 1000 1100
1200 1300 1400 1500
</DataItem>
</Attribute>
<Attribute Name="Cell Centered Values" Center="Cell">
<DataItem Format="XML" Dimensions="3">
100 200 400
</DataItem>
</Attribute>
<Set SetType="Node">
<DataItem NumberType="Int" Dimensions="3" Format="XML">
0 2 4
</DataItem>
<Attribute Name="Node Centered Set Values" Center="Node">
<DataItem Format="XML" Dimensions="3">
100 200 400
</DataItem>
</Attribute>
</Set>
<Set SetType="Face">
<DataItem NumberType="Int" Dimensions="2" Format="XML">
0 1
</DataItem>
<DataItem NumberType="Int" Dimensions="2" Format="XML">
0 0
</DataItem>
</Set>
<Set SetType="Edge">
<DataItem NumberType="Int" Dimensions="4" Format="XML">
0 0 1 1
</DataItem>
<DataItem NumberType="Int" Dimensions="4" Format="XML">
0 0 0 0
</DataItem>
<DataItem NumberType="Int" Dimensions="4" Format="XML">
0 1 0 2
</DataItem>
<Attribute Name="Set Values" Center="Edge">
<DataItem Format="XML" Dimensions="4">
100 200 400 600
</DataItem>
</Attribute>
<Attribute Name="Other Set Values" Center="Edge">
<DataItem Format="XML" Dimensions="3">
1000 2000 4000 6000
</DataItem>
</Attribute>
</Set>
</Grid>
</Domain>
</Xdmf>
"""
d = XdmfDOM()
d.Parse(txt)
dm = d.FindElement('Domain')
gnode = d.FindElement('Grid', 0, dm)
g = XdmfGrid()
g.SetDOM(d)
g.SetElement(gnode)
g.UpdateInformation()
g.Update()
print ('Grid Sets :', g.GetNumberOfSets())
for i in range(g.GetNumberOfSets()) :
s = g.GetSets(i)
print ('Set #', i)
print (' Type: ', s.GetSetTypeAsString())
print (' Size: ', s.GetSize())
s.Update()
ids = s.GetIds()
print (' # Ids : ', ids.GetNumberOfElements())
print (' Ids : ', ids.GetValues())
cellids = s.GetCellIds()
print (' # Cell Ids : ', cellids.GetNumberOfElements())
print (' Cell Ids : ', cellids.GetValues())
faceids = s.GetFaceIds()
print (' # Face Ids : ', faceids.GetNumberOfElements())
print (' Face Ids : ', faceids.GetValues())
print (' #Attributes: ', s.GetNumberOfAttributes())
for j in range(s.GetNumberOfAttributes()) :
a = s.GetAttribute(j)
print (' Attribute #', j)
print (' Attribute values: ', a.GetAttributeTypeAsString())
a.Update()
v = a.GetValues()
print (' Attribute type: ', v.GetValues())
ns = XdmfSet()
ns.SetDOM(d)
ns.SetSetType(XDMF_SET_TYPE_NODE)
a = XdmfArray()
a.SetNumberOfElements(4)
a.SetNumberType(XDMF_INT32_TYPE)
a.SetValues(0, "100 200 300 400")
ns.SetIds(a)
g.Insert(ns)
ns.Build()
print (d.Serialize())
|