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
|
#include <Xdmf.h>
#include <XdmfDiff.h>
XdmfDOM * createGrid(int * connections, double * points, double * nodeValues, double * cellValues)
{
XdmfRoot * myRoot = new XdmfRoot();
XdmfDomain * myDomain = new XdmfDomain();
XdmfDOM * myDOM = new XdmfDOM();
myRoot->SetDOM(myDOM);
myRoot->Build();
myRoot->Insert(myDomain);
XdmfGrid * myGrid = new XdmfGrid();
myGrid->SetName("test");
// Write Topology
XdmfTopology * myTopology = myGrid->GetTopology();
myTopology->SetTopologyType(XDMF_HEX);
myTopology->SetNumberOfElements(2);
XdmfArray * myConnections = myTopology->GetConnectivity();
myConnections->SetHeavyDataSetName("output.h5:/Connections");
myConnections->SetNumberType(XDMF_INT32_TYPE);
myConnections->SetNumberOfElements(16);
myConnections->SetValues(0, connections, 16);
// Write Geometry
XdmfGeometry * myGeometry = myGrid->GetGeometry();
myGeometry->SetGeometryType(XDMF_GEOMETRY_XYZ);
myGeometry->SetNumberOfPoints(12);
XdmfArray * myPoints = myGeometry->GetPoints();
myPoints->SetHeavyDataSetName("output.h5:/XYZ");
myPoints->SetNumberType(XDMF_FLOAT64_TYPE);
myPoints->SetNumberOfElements(36);
myPoints->SetValues(0, points, 36);
myDomain->Insert(myGrid);
XdmfAttribute * currAttribute = new XdmfAttribute();
currAttribute->SetName("NodeValues");
currAttribute->SetAttributeCenter(XDMF_ATTRIBUTE_CENTER_NODE);
currAttribute->SetAttributeType(XDMF_ATTRIBUTE_TYPE_SCALAR);
currAttribute->SetDeleteOnGridDelete(true);
XdmfArray * array = currAttribute->GetValues();
array->SetHeavyDataSetName("output.h5:/NodeValues");
// For now only support the data we know we are writing
array->SetNumberType(XDMF_FLOAT64_TYPE);
array->SetNumberOfElements(12);
array->SetValues(0,nodeValues,12);
myGrid->Insert(currAttribute);
XdmfAttribute * cellAttribute = new XdmfAttribute();
cellAttribute->SetName("CellValues");
cellAttribute->SetAttributeCenter(XDMF_ATTRIBUTE_CENTER_CELL);
cellAttribute->SetAttributeType(XDMF_ATTRIBUTE_TYPE_SCALAR);
cellAttribute->SetDeleteOnGridDelete(true);
XdmfArray * cellArray = cellAttribute->GetValues();
cellArray->SetHeavyDataSetName("output.h5:/CellValues");
// For now only support the data we know we are writing
cellArray->SetNumberType(XDMF_FLOAT64_TYPE);
cellArray->SetNumberOfElements(2);
cellArray->SetValues(0,cellValues,2);
myGrid->Insert(cellAttribute);
myGrid->Build();
delete myGrid;
delete myDomain;
delete myRoot;
return myDOM;
}
int main()
{
int connections[16] = {0,1,7,6,3,4,10,9,1,2,8,7,4,5,11,10};
double points[36] = {0,0,1,1,0,1,3,0,2,0,1,1,1,1,1,3,2,2,0,0,-1,1,0,-1,3,0,-2,0,1,-1,1,1,-1,3,2,-2};
double nodeValues[12] = {100, 200, 300, 300, 400, 500, 300, 400, 500, 500, 600, 700};
double cellValues[2] = {100, 200};
XdmfDOM * myDOM1 = createGrid(connections, points, nodeValues, cellValues);
nodeValues[0] = 110;
XdmfDOM * myDOM2 = createGrid(connections, points, nodeValues, cellValues);
XdmfDiff * myDiff = new XdmfDiff(myDOM1, myDOM2);
if (myDiff->AreEquivalent())
{
return -1;
}
myDiff->SetAbsoluteError(5);
if (myDiff->AreEquivalent())
{
return -1;
}
myDiff->SetAbsoluteError(11);
if (!myDiff->AreEquivalent())
{
return -1;
}
myDiff->SetRelativeError(.05);
if (myDiff->AreEquivalent())
{
return -1;
}
myDiff->SetRelativeError(.1);
if (!myDiff->AreEquivalent())
{
return -1;
}
myDiff->SetRelativeError(0);
myDiff->SetIgnoreAllAttributes(true);
if (!myDiff->AreEquivalent())
{
return -1;
}
delete myDiff;
return 0;
}
|
