/[escript]/trunk-mpi-branch/finley/test/python/run_visualization_interface.py
ViewVC logotype

Contents of /trunk-mpi-branch/finley/test/python/run_visualization_interface.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 150 - (show annotations)
Thu Sep 15 03:44:45 2005 UTC (14 years, 2 months ago) by jgs
Original Path: trunk/esys2/finley/test/python/test_visualization_interface.py
File MIME type: text/x-python
File size: 72848 byte(s)
Merge of development branch dev-02 back to main trunk on 2005-09-15

1 # $Id$
2
3 import sys
4 import unittest
5 from esys.escript import *
6 from esys.finley import ReadMesh
7
8 REERFENCE_FILE_DIR="test_meshes"
9
10 class Test_VisualizationInterface(unittest.TestCase):
11 def check_vtk(self,f,reference_f):
12 out_string=open(f).read().splitlines()
13 ref_string=open(REERFENCE_FILE_DIR+"/"+reference_f).read().splitlines()
14 c=0
15 for l in range(0,len(ref_string)):
16 if not ref_string[l].strip()[:2]=="<!":
17 self.failUnlessEqual(out_string[c].strip().replace("e-00","e+00"),ref_string[l].strip(),"line %d (%s) in vtk files does not match reference (%s)"%(c,out_string[c].strip(),ref_string[l].strip()))
18 c+=1
19
20 def check_dx(self,f,reference_f):
21 out_string=open(f).read().splitlines()
22 ref_string=open(REERFENCE_FILE_DIR+"/"+reference_f).read().splitlines()
23 c=0
24 for l in range(0,len(ref_string)):
25 if not ref_string[l].strip()[0]=="#":
26 self.failUnlessEqual(out_string[c].strip().replace("e-00","e+00"),ref_string[l].strip(),"line %d (%s) in dx file does not match reference (%s)"%(c,out_string[c].strip(),ref_string[l].strip()))
27 c+=1
28
29 class Test_VTKFiles(Test_VisualizationInterface):
30 # ======================================================================================================================
31 def test_hex_contact_2D_order1_ContinuousFunction_Scalar_vtk(self):
32 reference="hex_2D_o1_node_s.xml"
33 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
34 x=ContinuousFunction(dom).getX()
35 saveVTK("hex_contact_2D_order1_ContinuousFunction_Scalar.xml",data=x[0])
36 self.check_vtk("data.hex_contact_2D_order1_ContinuousFunction_Scalar.xml",reference)
37 def test_hex_contact_2D_order1_ContinuousFunction_Vector_vtk(self):
38 reference="hex_2D_o1_node_v.xml"
39 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
40 x=ContinuousFunction(dom).getX()
41 saveVTK("hex_contact_2D_order1_ContinuousFunction_Vector.xml",data=x[0]*[1.,2.])
42 self.check_vtk("data.hex_contact_2D_order1_ContinuousFunction_Vector.xml",reference)
43 def test_hex_contact_2D_order1_ContinuousFunction_Tensor_vtk(self):
44 reference="hex_2D_o1_node_t.xml"
45 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
46 x=ContinuousFunction(dom).getX()
47 saveVTK("hex_contact_2D_order1_ContinuousFunction_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
48 self.check_vtk("data.hex_contact_2D_order1_ContinuousFunction_Tensor.xml",reference)
49 def test_hex_contact_2D_order1_Solution_Scalar_vtk(self):
50 reference="hex_2D_o1_node_s.xml"
51 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
52 x=Solution(dom).getX()
53 saveVTK("hex_contact_2D_order1_Solution_Scalar.xml",data=x[0])
54 self.check_vtk("data.hex_contact_2D_order1_Solution_Scalar.xml",reference)
55 def test_hex_contact_2D_order1_Solution_Vector_vtk(self):
56 reference="hex_2D_o1_node_v.xml"
57 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
58 x=Solution(dom).getX()
59 saveVTK("hex_contact_2D_order1_Solution_Vector.xml",data=x[0]*[1.,2.])
60 self.check_vtk("data.hex_contact_2D_order1_Solution_Vector.xml",reference)
61 def test_hex_contact_2D_order1_Solution_Tensor_vtk(self):
62 reference="hex_2D_o1_node_t.xml"
63 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
64 x=Solution(dom).getX()
65 saveVTK("hex_contact_2D_order1_Solution_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
66 self.check_vtk("data.hex_contact_2D_order1_Solution_Tensor.xml",reference)
67 def test_hex_contact_2D_order1_ReducedSolution_Scalar_vtk(self):
68 reference="hex_2D_o1_node_s.xml"
69 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
70 x=ReducedSolution(dom).getX()
71 saveVTK("hex_contact_2D_order1_ReducedSolution_Scalar.xml",data=x[0])
72 self.check_vtk("data.hex_contact_2D_order1_ReducedSolution_Scalar.xml",reference)
73 def test_hex_contact_2D_order1_ReducedSolution_Vector_vtk(self):
74 reference="hex_2D_o1_node_v.xml"
75 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
76 x=ReducedSolution(dom).getX()
77 saveVTK("hex_contact_2D_order1_ReducedSolution_Vector.xml",data=x[0]*[1.,2.])
78 self.check_vtk("data.hex_contact_2D_order1_ReducedSolution_Vector.xml",reference)
79 def test_hex_contact_2D_order1_ReducedSolution_Tensor_vtk(self):
80 reference="hex_2D_o1_node_t.xml"
81 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
82 x=ReducedSolution(dom).getX()
83 saveVTK("hex_contact_2D_order1_ReducedSolution_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
84 self.check_vtk("data.hex_contact_2D_order1_ReducedSolution_Tensor.xml",reference)
85 def test_hex_contact_2D_order1_Function_Scalar_vtk(self):
86 reference="hex_2D_o1_cell_s.xml"
87 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
88 x=Function(dom).getX()
89 saveVTK("hex_contact_2D_order1_Function_Scalar.xml",data=x[0])
90 self.check_vtk("data.hex_contact_2D_order1_Function_Scalar.xml",reference)
91 def test_hex_contact_2D_order1_Function_Vector_vtk(self):
92 reference="hex_2D_o1_cell_v.xml"
93 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
94 x=Function(dom).getX()
95 saveVTK("hex_contact_2D_order1_Function_Vector.xml",data=x[0]*[1.,2.])
96 self.check_vtk("data.hex_contact_2D_order1_Function_Vector.xml",reference)
97 def test_hex_contact_2D_order1_Function_Tensor_vtk(self):
98 reference="hex_2D_o1_cell_t.xml"
99 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
100 x=Function(dom).getX()
101 saveVTK("hex_contact_2D_order1_Function_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
102 self.check_vtk("data.hex_contact_2D_order1_Function_Tensor.xml",reference)
103 def test_hex_contact_2D_order1_FunctionOnBoundary_Scalar_vtk(self):
104 reference="hex_2D_o1_boundary_s.xml"
105 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
106 x=FunctionOnBoundary(dom).getX()
107 saveVTK("hex_contact_2D_order1_FunctionOnBoundary_Scalar.xml",data=x[0])
108 self.check_vtk("data.hex_contact_2D_order1_FunctionOnBoundary_Scalar.xml",reference)
109 def test_hex_contact_2D_order1_FunctionOnBoundary_Vector_vtk(self):
110 reference="hex_2D_o1_boundary_v.xml"
111 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
112 x=FunctionOnBoundary(dom).getX()
113 saveVTK("hex_contact_2D_order1_FunctionOnBoundary_Vector.xml",data=x[0]*[1.,2.])
114 self.check_vtk("data.hex_contact_2D_order1_FunctionOnBoundary_Vector.xml",reference)
115 def test_hex_contact_2D_order1_FunctionOnBoundary_Tensor_vtk(self):
116 reference="hex_2D_o1_boundary_t.xml"
117 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
118 x=FunctionOnBoundary(dom).getX()
119 saveVTK("hex_contact_2D_order1_FunctionOnBoundary_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
120 self.check_vtk("data.hex_contact_2D_order1_FunctionOnBoundary_Tensor.xml",reference)
121 def test_hex_contact_2D_order1_onFace_FunctionOnBoundary_Scalar_vtk(self):
122 reference="hex_2D_o1_f_boundary_s.xml"
123 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1_onFace.msh")
124 x=FunctionOnBoundary(dom).getX()
125 saveVTK("hex_contact_2D_order1_onFace_FunctionOnBoundary_Scalar.xml",data=x[0])
126 self.check_vtk("data.hex_contact_2D_order1_onFace_FunctionOnBoundary_Scalar.xml",reference)
127 def test_hex_contact_2D_order1_onFace_FunctionOnBoundary_Vector_vtk(self):
128 reference="hex_2D_o1_f_boundary_v.xml"
129 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1_onFace.msh")
130 x=FunctionOnBoundary(dom).getX()
131 saveVTK("hex_contact_2D_order1_onFace_FunctionOnBoundary_Vector.xml",data=x[0]*[1.,2.])
132 self.check_vtk("data.hex_contact_2D_order1_onFace_FunctionOnBoundary_Vector.xml",reference)
133 def test_hex_contact_2D_order1_onFace_FunctionOnBoundary_Tensor_vtk(self):
134 reference="hex_2D_o1_f_boundary_t.xml"
135 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1_onFace.msh")
136 x=FunctionOnBoundary(dom).getX()
137 saveVTK("hex_contact_2D_order1_onFace_FunctionOnBoundary_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
138 self.check_vtk("data.hex_contact_2D_order1_onFace_FunctionOnBoundary_Tensor.xml",reference)
139 def test_hex_contact_2D_order1_FunctionOnContactZero_Scalar_vtk(self):
140 reference="hex_2D_o1_contact_s.xml"
141 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
142 x=FunctionOnContactZero(dom).getX()
143 saveVTK("hex_contact_2D_order1_FunctionOnContactZero_Scalar.xml",data=x[0])
144 self.check_vtk("data.hex_contact_2D_order1_FunctionOnContactZero_Scalar.xml",reference)
145 def test_hex_contact_2D_order1_FunctionOnContactZero_Vector_vtk(self):
146 reference="hex_2D_o1_contact_v.xml"
147 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
148 x=FunctionOnContactZero(dom).getX()
149 saveVTK("hex_contact_2D_order1_FunctionOnContactZero_Vector.xml",data=x[0]*[1.,2.])
150 self.check_vtk("data.hex_contact_2D_order1_FunctionOnContactZero_Vector.xml",reference)
151 def test_hex_contact_2D_order1_FunctionOnContactZero_Tensor_vtk(self):
152 reference="hex_2D_o1_contact_t.xml"
153 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
154 x=FunctionOnContactZero(dom).getX()
155 saveVTK("hex_contact_2D_order1_FunctionOnContactZero_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
156 self.check_vtk("data.hex_contact_2D_order1_FunctionOnContactZero_Tensor.xml",reference)
157 def test_hex_contact_2D_order1_onFace_FunctionOnContactZero_Scalar_vtk(self):
158 reference="hex_2D_o1_contact_s.xml"
159 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1_onFace.msh")
160 x=FunctionOnContactZero(dom).getX()
161 saveVTK("hex_contact_2D_order1_onFace_FunctionOnContactZero_Scalar.xml",data=x[0])
162 self.check_vtk("data.hex_contact_2D_order1_onFace_FunctionOnContactZero_Scalar.xml",reference)
163 def test_hex_contact_2D_order1_onFace_FunctionOnContactZero_Vector_vtk(self):
164 reference="hex_2D_o1_contact_v.xml"
165 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1_onFace.msh")
166 x=FunctionOnContactZero(dom).getX()
167 saveVTK("hex_contact_2D_order1_onFace_FunctionOnContactZero_Vector.xml",data=x[0]*[1.,2.])
168 self.check_vtk("data.hex_contact_2D_order1_onFace_FunctionOnContactZero_Vector.xml",reference)
169 def test_hex_contact_2D_order1_onFace_FunctionOnContactZero_Tensor_vtk(self):
170 reference="hex_2D_o1_contact_t.xml"
171 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1_onFace.msh")
172 x=FunctionOnContactZero(dom).getX()
173 saveVTK("hex_contact_2D_order1_onFace_FunctionOnContactZero_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
174 self.check_vtk("data.hex_contact_2D_order1_onFace_FunctionOnContactZero_Tensor.xml",reference)
175 def test_hex_contact_2D_order1_FunctionOnContactOne_Scalar_vtk(self):
176 reference="hex_2D_o1_contact_s.xml"
177 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
178 x=FunctionOnContactOne(dom).getX()
179 saveVTK("hex_contact_2D_order1_FunctionOnContactOne_Scalar.xml",data=x[0])
180 self.check_vtk("data.hex_contact_2D_order1_FunctionOnContactOne_Scalar.xml",reference)
181 def test_hex_contact_2D_order1_FunctionOnContactOne_Vector_vtk(self):
182 reference="hex_2D_o1_contact_v.xml"
183 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
184 x=FunctionOnContactOne(dom).getX()
185 saveVTK("hex_contact_2D_order1_FunctionOnContactOne_Vector.xml",data=x[0]*[1.,2.])
186 self.check_vtk("data.hex_contact_2D_order1_FunctionOnContactOne_Vector.xml",reference)
187 def test_hex_contact_2D_order1_FunctionOnContactOne_Tensor_vtk(self):
188 reference="hex_2D_o1_contact_t.xml"
189 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
190 x=FunctionOnContactOne(dom).getX()
191 saveVTK("hex_contact_2D_order1_FunctionOnContactOne_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
192 self.check_vtk("data.hex_contact_2D_order1_FunctionOnContactOne_Tensor.xml",reference)
193 def test_hex_contact_2D_order1_onFace_FunctionOnContactOne_Scalar_vtk(self):
194 reference="hex_2D_o1_contact_s.xml"
195 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1_onFace.msh")
196 x=FunctionOnContactOne(dom).getX()
197 saveVTK("hex_contact_2D_order1_onFace_FunctionOnContactOne_Scalar.xml",data=x[0])
198 self.check_vtk("data.hex_contact_2D_order1_onFace_FunctionOnContactOne_Scalar.xml",reference)
199 def test_hex_contact_2D_order1_onFace_FunctionOnContactOne_Vector_vtk(self):
200 reference="hex_2D_o1_contact_v.xml"
201 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1_onFace.msh")
202 x=FunctionOnContactOne(dom).getX()
203 saveVTK("hex_contact_2D_order1_onFace_FunctionOnContactOne_Vector.xml",data=x[0]*[1.,2.])
204 self.check_vtk("data.hex_contact_2D_order1_onFace_FunctionOnContactOne_Vector.xml",reference)
205 def test_hex_contact_2D_order1_onFace_FunctionOnContactOne_Tensor_vtk(self):
206 reference="hex_2D_o1_contact_t.xml"
207 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1_onFace.msh")
208 x=FunctionOnContactOne(dom).getX()
209 saveVTK("hex_contact_2D_order1_onFace_FunctionOnContactOne_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
210 self.check_vtk("data.hex_contact_2D_order1_onFace_FunctionOnContactOne_Tensor.xml",reference)
211 # ======================================================================================================================
212 def test_hex_contact_2D_order2_ContinuousFunction_Scalar_vtk(self):
213 reference="hex_2D_o2_node_s.xml"
214 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
215 x=ContinuousFunction(dom).getX()
216 saveVTK("hex_contact_2D_order2_ContinuousFunction_Scalar.xml",data=x[0])
217 self.check_vtk("data.hex_contact_2D_order2_ContinuousFunction_Scalar.xml",reference)
218 def test_hex_contact_2D_order2_ContinuousFunction_Vector_vtk(self):
219 reference="hex_2D_o2_node_v.xml"
220 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
221 x=ContinuousFunction(dom).getX()
222 saveVTK("hex_contact_2D_order2_ContinuousFunction_Vector.xml",data=x[0]*[1.,2.])
223 self.check_vtk("data.hex_contact_2D_order2_ContinuousFunction_Vector.xml",reference)
224 def test_hex_contact_2D_order2_ContinuousFunction_Tensor_vtk(self):
225 reference="hex_2D_o2_node_t.xml"
226 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
227 x=ContinuousFunction(dom).getX()
228 saveVTK("hex_contact_2D_order2_ContinuousFunction_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
229 self.check_vtk("data.hex_contact_2D_order2_ContinuousFunction_Tensor.xml",reference)
230 def test_hex_contact_2D_order2_Solution_Scalar_vtk(self):
231 reference="hex_2D_o2_node_s.xml"
232 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
233 x=Solution(dom).getX()
234 saveVTK("hex_contact_2D_order2_Solution_Scalar.xml",data=x[0])
235 self.check_vtk("data.hex_contact_2D_order2_Solution_Scalar.xml",reference)
236 def test_hex_contact_2D_order2_Solution_Vector_vtk(self):
237 reference="hex_2D_o2_node_v.xml"
238 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
239 x=Solution(dom).getX()
240 saveVTK("hex_contact_2D_order2_Solution_Vector.xml",data=x[0]*[1.,2.])
241 self.check_vtk("data.hex_contact_2D_order2_Solution_Vector.xml",reference)
242 def test_hex_contact_2D_order2_Solution_Tensor_vtk(self):
243 reference="hex_2D_o2_node_t.xml"
244 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
245 x=Solution(dom).getX()
246 saveVTK("hex_contact_2D_order2_Solution_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
247 self.check_vtk("data.hex_contact_2D_order2_Solution_Tensor.xml",reference)
248 def test_hex_contact_2D_order2_ReducedSolution_Scalar_vtk(self):
249 reference="hex_2D_o1_node_s.xml"
250 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
251 x=ReducedSolution(dom).getX()
252 saveVTK("hex_contact_2D_order2_ReducedSolution_Scalar.xml",data=x[0])
253 self.check_vtk("data.hex_contact_2D_order2_ReducedSolution_Scalar.xml",reference)
254 def test_hex_contact_2D_order2_ReducedSolution_Vector_vtk(self):
255 reference="hex_2D_o1_node_v.xml"
256 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
257 x=ReducedSolution(dom).getX()
258 saveVTK("hex_contact_2D_order2_ReducedSolution_Vector.xml",data=x[0]*[1.,2.])
259 self.check_vtk("data.hex_contact_2D_order2_ReducedSolution_Vector.xml",reference)
260 def test_hex_contact_2D_order2_ReducedSolution_Tensor_vtk(self):
261 reference="hex_2D_o1_node_t.xml"
262 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
263 x=ReducedSolution(dom).getX()
264 saveVTK("hex_contact_2D_order2_ReducedSolution_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
265 self.check_vtk("data.hex_contact_2D_order2_ReducedSolution_Tensor.xml",reference)
266 def test_hex_contact_2D_order2_Function_Scalar_vtk(self):
267 reference="hex_2D_o2_cell_s.xml"
268 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
269 x=Function(dom).getX()
270 saveVTK("hex_contact_2D_order2_Function_Scalar.xml",data=x[0])
271 self.check_vtk("data.hex_contact_2D_order2_Function_Scalar.xml",reference)
272 def test_hex_contact_2D_order2_Function_Vector_vtk(self):
273 reference="hex_2D_o2_cell_v.xml"
274 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
275 x=Function(dom).getX()
276 saveVTK("hex_contact_2D_order2_Function_Vector.xml",data=x[0]*[1.,2.])
277 self.check_vtk("data.hex_contact_2D_order2_Function_Vector.xml",reference)
278 def test_hex_contact_2D_order2_Function_Tensor_vtk(self):
279 reference="hex_2D_o2_cell_t.xml"
280 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
281 x=Function(dom).getX()
282 saveVTK("hex_contact_2D_order2_Function_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
283 self.check_vtk("data.hex_contact_2D_order2_Function_Tensor.xml",reference)
284 def test_hex_contact_2D_order2_FunctionOnBoundary_Scalar_vtk(self):
285 reference="hex_2D_o2_boundary_s.xml"
286 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
287 x=FunctionOnBoundary(dom).getX()
288 saveVTK("hex_contact_2D_order2_FunctionOnBoundary_Scalar.xml",data=x[0])
289 self.check_vtk("data.hex_contact_2D_order2_FunctionOnBoundary_Scalar.xml",reference)
290 def test_hex_contact_2D_order2_FunctionOnBoundary_Vector_vtk(self):
291 reference="hex_2D_o2_boundary_v.xml"
292 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
293 x=FunctionOnBoundary(dom).getX()
294 saveVTK("hex_contact_2D_order2_FunctionOnBoundary_Vector.xml",data=x[0]*[1.,2.])
295 self.check_vtk("data.hex_contact_2D_order2_FunctionOnBoundary_Vector.xml",reference)
296 def test_hex_contact_2D_order2_FunctionOnBoundary_Tensor_vtk(self):
297 reference="hex_2D_o2_boundary_t.xml"
298 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
299 x=FunctionOnBoundary(dom).getX()
300 saveVTK("hex_contact_2D_order2_FunctionOnBoundary_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
301 self.check_vtk("data.hex_contact_2D_order2_FunctionOnBoundary_Tensor.xml",reference)
302 def test_hex_contact_2D_order2_onFace_FunctionOnBoundary_Scalar_vtk(self):
303 reference="hex_2D_o2_f_boundary_s.xml"
304 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2_onFace.msh")
305 x=FunctionOnBoundary(dom).getX()
306 saveVTK("hex_contact_2D_order2_onFace_FunctionOnBoundary_Scalar.xml",data=x[0])
307 self.check_vtk("data.hex_contact_2D_order2_onFace_FunctionOnBoundary_Scalar.xml",reference)
308 def test_hex_contact_2D_order2_onFace_FunctionOnBoundary_Vector_vtk(self):
309 reference="hex_2D_o2_f_boundary_v.xml"
310 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2_onFace.msh")
311 x=FunctionOnBoundary(dom).getX()
312 saveVTK("hex_contact_2D_order2_onFace_FunctionOnBoundary_Vector.xml",data=x[0]*[1.,2.])
313 self.check_vtk("data.hex_contact_2D_order2_onFace_FunctionOnBoundary_Vector.xml",reference)
314 def test_hex_contact_2D_order2_onFace_FunctionOnBoundary_Tensor_vtk(self):
315 reference="hex_2D_o2_f_boundary_t.xml"
316 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2_onFace.msh")
317 x=FunctionOnBoundary(dom).getX()
318 saveVTK("hex_contact_2D_order2_onFace_FunctionOnBoundary_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
319 self.check_vtk("data.hex_contact_2D_order2_onFace_FunctionOnBoundary_Tensor.xml",reference)
320 def test_hex_contact_2D_order2_FunctionOnContactZero_Scalar_vtk(self):
321 reference="hex_2D_o2_contact_s.xml"
322 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
323 x=FunctionOnContactZero(dom).getX()
324 saveVTK("hex_contact_2D_order2_FunctionOnContactZero_Scalar.xml",data=x[0])
325 self.check_vtk("data.hex_contact_2D_order2_FunctionOnContactZero_Scalar.xml",reference)
326 def test_hex_contact_2D_order2_FunctionOnContactZero_Vector_vtk(self):
327 reference="hex_2D_o2_contact_v.xml"
328 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
329 x=FunctionOnContactZero(dom).getX()
330 saveVTK("hex_contact_2D_order2_FunctionOnContactZero_Vector.xml",data=x[0]*[1.,2.])
331 self.check_vtk("data.hex_contact_2D_order2_FunctionOnContactZero_Vector.xml",reference)
332 def test_hex_contact_2D_order2_FunctionOnContactZero_Tensor_vtk(self):
333 reference="hex_2D_o2_contact_t.xml"
334 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
335 x=FunctionOnContactZero(dom).getX()
336 saveVTK("hex_contact_2D_order2_FunctionOnContactZero_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
337 self.check_vtk("data.hex_contact_2D_order2_FunctionOnContactZero_Tensor.xml",reference)
338 def test_hex_contact_2D_order2_onFace_FunctionOnContactZero_Scalar_vtk(self):
339 reference="hex_2D_o2_contact_s.xml"
340 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2_onFace.msh")
341 x=FunctionOnContactZero(dom).getX()
342 saveVTK("hex_contact_2D_order2_onFace_FunctionOnContactZero_Scalar.xml",data=x[0])
343 self.check_vtk("data.hex_contact_2D_order2_onFace_FunctionOnContactZero_Scalar.xml",reference)
344 def test_hex_contact_2D_order2_onFace_FunctionOnContactZero_Vector_vtk(self):
345 reference="hex_2D_o2_contact_v.xml"
346 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2_onFace.msh")
347 x=FunctionOnContactZero(dom).getX()
348 saveVTK("hex_contact_2D_order2_onFace_FunctionOnContactZero_Vector.xml",data=x[0]*[1.,2.])
349 self.check_vtk("data.hex_contact_2D_order2_onFace_FunctionOnContactZero_Vector.xml",reference)
350 def test_hex_contact_2D_order2_onFace_FunctionOnContactZero_Tensor_vtk(self):
351 reference="hex_2D_o2_contact_t.xml"
352 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2_onFace.msh")
353 x=FunctionOnContactZero(dom).getX()
354 saveVTK("hex_contact_2D_order2_onFace_FunctionOnContactZero_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
355 self.check_vtk("data.hex_contact_2D_order2_onFace_FunctionOnContactZero_Tensor.xml",reference)
356 def test_hex_contact_2D_order2_FunctionOnContactOne_Scalar_vtk(self):
357 reference="hex_2D_o2_contact_s.xml"
358 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
359 x=FunctionOnContactOne(dom).getX()
360 saveVTK("hex_contact_2D_order2_FunctionOnContactOne_Scalar.xml",data=x[0])
361 self.check_vtk("data.hex_contact_2D_order2_FunctionOnContactOne_Scalar.xml",reference)
362 def test_hex_contact_2D_order2_FunctionOnContactOne_Vector_vtk(self):
363 reference="hex_2D_o2_contact_v.xml"
364 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
365 x=FunctionOnContactOne(dom).getX()
366 saveVTK("hex_contact_2D_order2_FunctionOnContactOne_Vector.xml",data=x[0]*[1.,2.])
367 self.check_vtk("data.hex_contact_2D_order2_FunctionOnContactOne_Vector.xml",reference)
368 def test_hex_contact_2D_order2_FunctionOnContactOne_Tensor_vtk(self):
369 reference="hex_2D_o2_contact_t.xml"
370 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
371 x=FunctionOnContactOne(dom).getX()
372 saveVTK("hex_contact_2D_order2_FunctionOnContactOne_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
373 self.check_vtk("data.hex_contact_2D_order2_FunctionOnContactOne_Tensor.xml",reference)
374 def test_hex_contact_2D_order2_onFace_FunctionOnContactOne_Scalar_vtk(self):
375 reference="hex_2D_o2_contact_s.xml"
376 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2_onFace.msh")
377 x=FunctionOnContactOne(dom).getX()
378 saveVTK("hex_contact_2D_order2_onFace_FunctionOnContactOne_Scalar.xml",data=x[0])
379 self.check_vtk("data.hex_contact_2D_order2_onFace_FunctionOnContactOne_Scalar.xml",reference)
380 def test_hex_contact_2D_order2_onFace_FunctionOnContactOne_Vector_vtk(self):
381 reference="hex_2D_o2_contact_v.xml"
382 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2_onFace.msh")
383 x=FunctionOnContactOne(dom).getX()
384 saveVTK("hex_contact_2D_order2_onFace_FunctionOnContactOne_Vector.xml",data=x[0]*[1.,2.])
385 self.check_vtk("data.hex_contact_2D_order2_onFace_FunctionOnContactOne_Vector.xml",reference)
386 def test_hex_contact_2D_order2_onFace_FunctionOnContactOne_Tensor_vtk(self):
387 reference="hex_2D_o2_contact_t.xml"
388 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2_onFace.msh")
389 x=FunctionOnContactOne(dom).getX()
390 saveVTK("hex_contact_2D_order2_onFace_FunctionOnContactOne_Tensor.xml",data=x[0]*[[11.,12.],[21.,22.]])
391 self.check_vtk("data.hex_contact_2D_order2_onFace_FunctionOnContactOne_Tensor.xml",reference)
392
393 # ======================================================================================================================
394 def test_hex_contact_3D_order1_ContinuousFunction_Scalar_vtk(self):
395 reference="hex_3D_o1_node_s.xml"
396 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
397 x=ContinuousFunction(dom).getX()
398 saveVTK("hex_contact_3D_order1_ContinuousFunction_Scalar.xml",data=x[0])
399 self.check_vtk("data.hex_contact_3D_order1_ContinuousFunction_Scalar.xml",reference)
400 def test_hex_contact_3D_order1_ContinuousFunction_Vector_vtk(self):
401 reference="hex_3D_o1_node_v.xml"
402 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
403 x=ContinuousFunction(dom).getX()
404 saveVTK("hex_contact_3D_order1_ContinuousFunction_Vector.xml",data=x[0]*[1.,2.,3.])
405 self.check_vtk("data.hex_contact_3D_order1_ContinuousFunction_Vector.xml",reference)
406 def test_hex_contact_3D_order1_ContinuousFunction_Tensor_vtk(self):
407 reference="hex_3D_o1_node_t.xml"
408 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
409 x=ContinuousFunction(dom).getX()
410 saveVTK("hex_contact_3D_order1_ContinuousFunction_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
411 self.check_vtk("data.hex_contact_3D_order1_ContinuousFunction_Tensor.xml",reference)
412 def test_hex_contact_3D_order1_Solution_Scalar_vtk(self):
413 reference="hex_3D_o1_node_s.xml"
414 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
415 x=Solution(dom).getX()
416 saveVTK("hex_contact_3D_order1_Solution_Scalar.xml",data=x[0])
417 self.check_vtk("data.hex_contact_3D_order1_Solution_Scalar.xml",reference)
418 def test_hex_contact_3D_order1_Solution_Vector_vtk(self):
419 reference="hex_3D_o1_node_v.xml"
420 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
421 x=Solution(dom).getX()
422 saveVTK("hex_contact_3D_order1_Solution_Vector.xml",data=x[0]*[1.,2.,3.])
423 self.check_vtk("data.hex_contact_3D_order1_Solution_Vector.xml",reference)
424 def test_hex_contact_3D_order1_Solution_Tensor_vtk(self):
425 reference="hex_3D_o1_node_t.xml"
426 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
427 x=Solution(dom).getX()
428 saveVTK("hex_contact_3D_order1_Solution_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
429 self.check_vtk("data.hex_contact_3D_order1_Solution_Tensor.xml",reference)
430 def test_hex_contact_3D_order1_ReducedSolution_Scalar_vtk(self):
431 reference="hex_3D_o1_node_s.xml"
432 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
433 x=ReducedSolution(dom).getX()
434 saveVTK("hex_contact_3D_order1_ReducedSolution_Scalar.xml",data=x[0])
435 self.check_vtk("data.hex_contact_3D_order1_ReducedSolution_Scalar.xml",reference)
436 def test_hex_contact_3D_order1_ReducedSolution_Vector_vtk(self):
437 reference="hex_3D_o1_node_v.xml"
438 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
439 x=ReducedSolution(dom).getX()
440 saveVTK("hex_contact_3D_order1_ReducedSolution_Vector.xml",data=x[0]*[1.,2.,3.])
441 self.check_vtk("data.hex_contact_3D_order1_ReducedSolution_Vector.xml",reference)
442 def test_hex_contact_3D_order1_ReducedSolution_Tensor_vtk(self):
443 reference="hex_3D_o1_node_t.xml"
444 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
445 x=ReducedSolution(dom).getX()
446 saveVTK("hex_contact_3D_order1_ReducedSolution_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
447 self.check_vtk("data.hex_contact_3D_order1_ReducedSolution_Tensor.xml",reference)
448 def test_hex_contact_3D_order1_Function_Scalar_vtk(self):
449 reference="hex_3D_o1_cell_s.xml"
450 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
451 x=Function(dom).getX()
452 saveVTK("hex_contact_3D_order1_Function_Scalar.xml",data=x[0])
453 self.check_vtk("data.hex_contact_3D_order1_Function_Scalar.xml",reference)
454 def test_hex_contact_3D_order1_Function_Vector_vtk(self):
455 reference="hex_3D_o1_cell_v.xml"
456 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
457 x=Function(dom).getX()
458 saveVTK("hex_contact_3D_order1_Function_Vector.xml",data=x[0]*[1.,2.,3.])
459 self.check_vtk("data.hex_contact_3D_order1_Function_Vector.xml",reference)
460 def test_hex_contact_3D_order1_Function_Tensor_vtk(self):
461 reference="hex_3D_o1_cell_t.xml"
462 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
463 x=Function(dom).getX()
464 saveVTK("hex_contact_3D_order1_Function_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
465 self.check_vtk("data.hex_contact_3D_order1_Function_Tensor.xml",reference)
466 def test_hex_contact_3D_order1_FunctionOnBoundary_Scalar_vtk(self):
467 reference="hex_3D_o1_boundary_s.xml"
468 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
469 x=FunctionOnBoundary(dom).getX()
470 saveVTK("hex_contact_3D_order1_FunctionOnBoundary_Scalar.xml",data=x[0])
471 self.check_vtk("data.hex_contact_3D_order1_FunctionOnBoundary_Scalar.xml",reference)
472 def test_hex_contact_3D_order1_FunctionOnBoundary_Vector_vtk(self):
473 reference="hex_3D_o1_boundary_v.xml"
474 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
475 x=FunctionOnBoundary(dom).getX()
476 saveVTK("hex_contact_3D_order1_FunctionOnBoundary_Vector.xml",data=x[0]*[1.,2.,3.])
477 self.check_vtk("data.hex_contact_3D_order1_FunctionOnBoundary_Vector.xml",reference)
478 def test_hex_contact_3D_order1_FunctionOnBoundary_Tensor_vtk(self):
479 reference="hex_3D_o1_boundary_t.xml"
480 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
481 x=FunctionOnBoundary(dom).getX()
482 saveVTK("hex_contact_3D_order1_FunctionOnBoundary_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
483 self.check_vtk("data.hex_contact_3D_order1_FunctionOnBoundary_Tensor.xml",reference)
484 def test_hex_contact_3D_order1_onFace_FunctionOnBoundary_Scalar_vtk(self):
485 reference="hex_3D_o1_f_boundary_s.xml"
486 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1_onFace.msh")
487 x=FunctionOnBoundary(dom).getX()
488 saveVTK("hex_contact_3D_order1_onFace_FunctionOnBoundary_Scalar.xml",data=x[0])
489 self.check_vtk("data.hex_contact_3D_order1_onFace_FunctionOnBoundary_Scalar.xml",reference)
490 def test_hex_contact_3D_order1_onFace_FunctionOnBoundary_Vector_vtk(self):
491 reference="hex_3D_o1_f_boundary_v.xml"
492 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1_onFace.msh")
493 x=FunctionOnBoundary(dom).getX()
494 saveVTK("hex_contact_3D_order1_onFace_FunctionOnBoundary_Vector.xml",data=x[0]*[1.,2.,3.])
495 self.check_vtk("data.hex_contact_3D_order1_onFace_FunctionOnBoundary_Vector.xml",reference)
496 def test_hex_contact_3D_order1_onFace_FunctionOnBoundary_Tensor_vtk(self):
497 reference="hex_3D_o1_f_boundary_t.xml"
498 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1_onFace.msh")
499 x=FunctionOnBoundary(dom).getX()
500 saveVTK("hex_contact_3D_order1_onFace_FunctionOnBoundary_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
501 self.check_vtk("data.hex_contact_3D_order1_onFace_FunctionOnBoundary_Tensor.xml",reference)
502 def test_hex_contact_3D_order1_FunctionOnContactZero_Scalar_vtk(self):
503 reference="hex_3D_o1_contact_s.xml"
504 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
505 x=FunctionOnContactZero(dom).getX()
506 saveVTK("hex_contact_3D_order1_FunctionOnContactZero_Scalar.xml",data=x[0])
507 self.check_vtk("data.hex_contact_3D_order1_FunctionOnContactZero_Scalar.xml",reference)
508 def test_hex_contact_3D_order1_FunctionOnContactZero_Vector_vtk(self):
509 reference="hex_3D_o1_contact_v.xml"
510 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
511 x=FunctionOnContactZero(dom).getX()
512 saveVTK("hex_contact_3D_order1_FunctionOnContactZero_Vector.xml",data=x[0]*[1.,2.,3.])
513 self.check_vtk("data.hex_contact_3D_order1_FunctionOnContactZero_Vector.xml",reference)
514 def test_hex_contact_3D_order1_FunctionOnContactZero_Tensor_vtk(self):
515 reference="hex_3D_o1_contact_t.xml"
516 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
517 x=FunctionOnContactZero(dom).getX()
518 saveVTK("hex_contact_3D_order1_FunctionOnContactZero_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
519 self.check_vtk("data.hex_contact_3D_order1_FunctionOnContactZero_Tensor.xml",reference)
520 def test_hex_contact_3D_order1_onFace_FunctionOnContactZero_Scalar_vtk(self):
521 reference="hex_3D_o1_contact_s.xml"
522 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1_onFace.msh")
523 x=FunctionOnContactZero(dom).getX()
524 saveVTK("hex_contact_3D_order1_onFace_FunctionOnContactZero_Scalar.xml",data=x[0])
525 self.check_vtk("data.hex_contact_3D_order1_onFace_FunctionOnContactZero_Scalar.xml",reference)
526 def test_hex_contact_3D_order1_onFace_FunctionOnContactZero_Vector_vtk(self):
527 reference="hex_3D_o1_contact_v.xml"
528 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1_onFace.msh")
529 x=FunctionOnContactZero(dom).getX()
530 saveVTK("hex_contact_3D_order1_onFace_FunctionOnContactZero_Vector.xml",data=x[0]*[1.,2.,3.])
531 self.check_vtk("data.hex_contact_3D_order1_onFace_FunctionOnContactZero_Vector.xml",reference)
532 def test_hex_contact_3D_order1_onFace_FunctionOnContactZero_Tensor_vtk(self):
533 reference="hex_3D_o1_contact_t.xml"
534 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1_onFace.msh")
535 x=FunctionOnContactZero(dom).getX()
536 saveVTK("hex_contact_3D_order1_onFace_FunctionOnContactZero_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
537 self.check_vtk("data.hex_contact_3D_order1_onFace_FunctionOnContactZero_Tensor.xml",reference)
538 def test_hex_contact_3D_order1_FunctionOnContactOne_Scalar_vtk(self):
539 reference="hex_3D_o1_contact_s.xml"
540 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
541 x=FunctionOnContactOne(dom).getX()
542 saveVTK("hex_contact_3D_order1_FunctionOnContactOne_Scalar.xml",data=x[0])
543 self.check_vtk("data.hex_contact_3D_order1_FunctionOnContactOne_Scalar.xml",reference)
544 def test_hex_contact_3D_order1_FunctionOnContactOne_Vector_vtk(self):
545 reference="hex_3D_o1_contact_v.xml"
546 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
547 x=FunctionOnContactOne(dom).getX()
548 saveVTK("hex_contact_3D_order1_FunctionOnContactOne_Vector.xml",data=x[0]*[1.,2.,3.])
549 self.check_vtk("data.hex_contact_3D_order1_FunctionOnContactOne_Vector.xml",reference)
550 def test_hex_contact_3D_order1_FunctionOnContactOne_Tensor_vtk(self):
551 reference="hex_3D_o1_contact_t.xml"
552 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
553 x=FunctionOnContactOne(dom).getX()
554 saveVTK("hex_contact_3D_order1_FunctionOnContactOne_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
555 self.check_vtk("data.hex_contact_3D_order1_FunctionOnContactOne_Tensor.xml",reference)
556 def test_hex_contact_3D_order1_onFace_FunctionOnContactOne_Scalar_vtk(self):
557 reference="hex_3D_o1_contact_s.xml"
558 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1_onFace.msh")
559 x=FunctionOnContactOne(dom).getX()
560 saveVTK("hex_contact_3D_order1_onFace_FunctionOnContactOne_Scalar.xml",data=x[0])
561 self.check_vtk("data.hex_contact_3D_order1_onFace_FunctionOnContactOne_Scalar.xml",reference)
562 def test_hex_contact_3D_order1_onFace_FunctionOnContactOne_Vector_vtk(self):
563 reference="hex_3D_o1_contact_v.xml"
564 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1_onFace.msh")
565 x=FunctionOnContactOne(dom).getX()
566 saveVTK("hex_contact_3D_order1_onFace_FunctionOnContactOne_Vector.xml",data=x[0]*[1.,2.,3.])
567 self.check_vtk("data.hex_contact_3D_order1_onFace_FunctionOnContactOne_Vector.xml",reference)
568 def test_hex_contact_3D_order1_onFace_FunctionOnContactOne_Tensor_vtk(self):
569 reference="hex_3D_o1_contact_t.xml"
570 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1_onFace.msh")
571 x=FunctionOnContactOne(dom).getX()
572 saveVTK("hex_contact_3D_order1_onFace_FunctionOnContactOne_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
573 self.check_vtk("data.hex_contact_3D_order1_onFace_FunctionOnContactOne_Tensor.xml",reference)
574 # ======================================================================================================================
575 def test_hex_contact_3D_order2_ContinuousFunction_Scalar_vtk(self):
576 reference="hex_3D_o2_node_s.xml"
577 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
578 x=ContinuousFunction(dom).getX()
579 saveVTK("hex_contact_3D_order2_ContinuousFunction_Scalar.xml",data=x[0])
580 self.check_vtk("data.hex_contact_3D_order2_ContinuousFunction_Scalar.xml",reference)
581 def test_hex_contact_3D_order2_ContinuousFunction_Vector_vtk(self):
582 reference="hex_3D_o2_node_v.xml"
583 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
584 x=ContinuousFunction(dom).getX()
585 saveVTK("hex_contact_3D_order2_ContinuousFunction_Vector.xml",data=x[0]*[1.,2.,3.])
586 self.check_vtk("data.hex_contact_3D_order2_ContinuousFunction_Vector.xml",reference)
587 def test_hex_contact_3D_order2_ContinuousFunction_Tensor_vtk(self):
588 reference="hex_3D_o2_node_t.xml"
589 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
590 x=ContinuousFunction(dom).getX()
591 saveVTK("hex_contact_3D_order2_ContinuousFunction_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
592 self.check_vtk("data.hex_contact_3D_order2_ContinuousFunction_Tensor.xml",reference)
593 def test_hex_contact_3D_order2_Solution_Scalar_vtk(self):
594 reference="hex_3D_o2_node_s.xml"
595 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
596 x=Solution(dom).getX()
597 saveVTK("hex_contact_3D_order2_Solution_Scalar.xml",data=x[0])
598 self.check_vtk("data.hex_contact_3D_order2_Solution_Scalar.xml",reference)
599 def test_hex_contact_3D_order2_Solution_Vector_vtk(self):
600 reference="hex_3D_o2_node_v.xml"
601 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
602 x=Solution(dom).getX()
603 saveVTK("hex_contact_3D_order2_Solution_Vector.xml",data=x[0]*[1.,2.,3.])
604 self.check_vtk("data.hex_contact_3D_order2_Solution_Vector.xml",reference)
605 def test_hex_contact_3D_order2_Solution_Tensor_vtk(self):
606 reference="hex_3D_o2_node_t.xml"
607 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
608 x=Solution(dom).getX()
609 saveVTK("hex_contact_3D_order2_Solution_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
610 self.check_vtk("data.hex_contact_3D_order2_Solution_Tensor.xml",reference)
611 def test_hex_contact_3D_order2_ReducedSolution_Scalar_vtk(self):
612 reference="hex_3D_o1_node_s.xml"
613 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
614 x=ReducedSolution(dom).getX()
615 saveVTK("hex_contact_3D_order2_ReducedSolution_Scalar.xml",data=x[0])
616 self.check_vtk("data.hex_contact_3D_order2_ReducedSolution_Scalar.xml",reference)
617 def test_hex_contact_3D_order2_ReducedSolution_Vector_vtk(self):
618 reference="hex_3D_o1_node_v.xml"
619 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
620 x=ReducedSolution(dom).getX()
621 saveVTK("hex_contact_3D_order2_ReducedSolution_Vector.xml",data=x[0]*[1.,2.,3.])
622 self.check_vtk("data.hex_contact_3D_order2_ReducedSolution_Vector.xml",reference)
623 def test_hex_contact_3D_order2_ReducedSolution_Tensor_vtk(self):
624 reference="hex_3D_o1_node_t.xml"
625 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
626 x=ReducedSolution(dom).getX()
627 saveVTK("hex_contact_3D_order2_ReducedSolution_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
628 self.check_vtk("data.hex_contact_3D_order2_ReducedSolution_Tensor.xml",reference)
629 def test_hex_contact_3D_order2_Function_Scalar_vtk(self):
630 reference="hex_3D_o2_cell_s.xml"
631 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
632 x=Function(dom).getX()
633 saveVTK("hex_contact_3D_order2_Function_Scalar.xml",data=x[0])
634 self.check_vtk("data.hex_contact_3D_order2_Function_Scalar.xml",reference)
635 def test_hex_contact_3D_order2_Function_Vector_vtk(self):
636 reference="hex_3D_o2_cell_v.xml"
637 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
638 x=Function(dom).getX()
639 saveVTK("hex_contact_3D_order2_Function_Vector.xml",data=x[0]*[1.,2.,3.])
640 self.check_vtk("data.hex_contact_3D_order2_Function_Vector.xml",reference)
641 def test_hex_contact_3D_order2_Function_Tensor_vtk(self):
642 reference="hex_3D_o2_cell_t.xml"
643 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
644 x=Function(dom).getX()
645 saveVTK("hex_contact_3D_order2_Function_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
646 self.check_vtk("data.hex_contact_3D_order2_Function_Tensor.xml",reference)
647 def test_hex_contact_3D_order2_FunctionOnBoundary_Scalar_vtk(self):
648 reference="hex_3D_o2_boundary_s.xml"
649 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
650 x=FunctionOnBoundary(dom).getX()
651 saveVTK("hex_contact_3D_order2_FunctionOnBoundary_Scalar.xml",data=x[0])
652 self.check_vtk("data.hex_contact_3D_order2_FunctionOnBoundary_Scalar.xml",reference)
653 def test_hex_contact_3D_order2_FunctionOnBoundary_Vector_vtk(self):
654 reference="hex_3D_o2_boundary_v.xml"
655 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
656 x=FunctionOnBoundary(dom).getX()
657 saveVTK("hex_contact_3D_order2_FunctionOnBoundary_Vector.xml",data=x[0]*[1.,2.,3.])
658 self.check_vtk("data.hex_contact_3D_order2_FunctionOnBoundary_Vector.xml",reference)
659 def test_hex_contact_3D_order2_FunctionOnBoundary_Tensor_vtk(self):
660 reference="hex_3D_o2_boundary_t.xml"
661 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
662 x=FunctionOnBoundary(dom).getX()
663 saveVTK("hex_contact_3D_order2_FunctionOnBoundary_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
664 self.check_vtk("data.hex_contact_3D_order2_FunctionOnBoundary_Tensor.xml",reference)
665 def test_hex_contact_3D_order2_onFace_FunctionOnBoundary_Scalar_vtk(self):
666 reference="hex_3D_o2_f_boundary_s.xml"
667 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2_onFace.msh")
668 x=FunctionOnBoundary(dom).getX()
669 saveVTK("hex_contact_3D_order2_onFace_FunctionOnBoundary_Scalar.xml",data=x[0])
670 self.check_vtk("data.hex_contact_3D_order2_onFace_FunctionOnBoundary_Scalar.xml",reference)
671 def test_hex_contact_3D_order2_onFace_FunctionOnBoundary_Vector_vtk(self):
672 reference="hex_3D_o2_f_boundary_v.xml"
673 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2_onFace.msh")
674 x=FunctionOnBoundary(dom).getX()
675 saveVTK("hex_contact_3D_order2_onFace_FunctionOnBoundary_Vector.xml",data=x[0]*[1.,2.,3.])
676 self.check_vtk("data.hex_contact_3D_order2_onFace_FunctionOnBoundary_Vector.xml",reference)
677 def test_hex_contact_3D_order2_onFace_FunctionOnBoundary_Tensor_vtk(self):
678 reference="hex_3D_o2_f_boundary_t.xml"
679 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2_onFace.msh")
680 x=FunctionOnBoundary(dom).getX()
681 saveVTK("hex_contact_3D_order2_onFace_FunctionOnBoundary_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
682 self.check_vtk("data.hex_contact_3D_order2_onFace_FunctionOnBoundary_Tensor.xml",reference)
683 def test_hex_contact_3D_order2_FunctionOnContactZero_Scalar_vtk(self):
684 reference="hex_3D_o2_contact_s.xml"
685 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
686 x=FunctionOnContactZero(dom).getX()
687 saveVTK("hex_contact_3D_order2_FunctionOnContactZero_Scalar.xml",data=x[0])
688 self.check_vtk("data.hex_contact_3D_order2_FunctionOnContactZero_Scalar.xml",reference)
689 def test_hex_contact_3D_order2_FunctionOnContactZero_Vector_vtk(self):
690 reference="hex_3D_o2_contact_v.xml"
691 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
692 x=FunctionOnContactZero(dom).getX()
693 saveVTK("hex_contact_3D_order2_FunctionOnContactZero_Vector.xml",data=x[0]*[1.,2.,3.])
694 self.check_vtk("data.hex_contact_3D_order2_FunctionOnContactZero_Vector.xml",reference)
695 def test_hex_contact_3D_order2_FunctionOnContactZero_Tensor_vtk(self):
696 reference="hex_3D_o2_contact_t.xml"
697 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
698 x=FunctionOnContactZero(dom).getX()
699 saveVTK("hex_contact_3D_order2_FunctionOnContactZero_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
700 self.check_vtk("data.hex_contact_3D_order2_FunctionOnContactZero_Tensor.xml",reference)
701 def test_hex_contact_3D_order2_onFace_FunctionOnContactZero_Scalar_vtk(self):
702 reference="hex_3D_o2_contact_s.xml"
703 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2_onFace.msh")
704 x=FunctionOnContactZero(dom).getX()
705 saveVTK("hex_contact_3D_order2_onFace_FunctionOnContactZero_Scalar.xml",data=x[0])
706 self.check_vtk("data.hex_contact_3D_order2_onFace_FunctionOnContactZero_Scalar.xml",reference)
707 def test_hex_contact_3D_order2_onFace_FunctionOnContactZero_Vector_vtk(self):
708 reference="hex_3D_o2_contact_v.xml"
709 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2_onFace.msh")
710 x=FunctionOnContactZero(dom).getX()
711 saveVTK("hex_contact_3D_order2_onFace_FunctionOnContactZero_Vector.xml",data=x[0]*[1.,2.,3.])
712 self.check_vtk("data.hex_contact_3D_order2_onFace_FunctionOnContactZero_Vector.xml",reference)
713 def test_hex_contact_3D_order2_onFace_FunctionOnContactZero_Tensor_vtk(self):
714 reference="hex_3D_o2_contact_t.xml"
715 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2_onFace.msh")
716 x=FunctionOnContactZero(dom).getX()
717 saveVTK("hex_contact_3D_order2_onFace_FunctionOnContactZero_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
718 self.check_vtk("data.hex_contact_3D_order2_onFace_FunctionOnContactZero_Tensor.xml",reference)
719 def test_hex_contact_3D_order2_FunctionOnContactOne_Scalar_vtk(self):
720 reference="hex_3D_o2_contact_s.xml"
721 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
722 x=FunctionOnContactOne(dom).getX()
723 saveVTK("hex_contact_3D_order2_FunctionOnContactOne_Scalar.xml",data=x[0])
724 self.check_vtk("data.hex_contact_3D_order2_FunctionOnContactOne_Scalar.xml",reference)
725 def test_hex_contact_3D_order2_FunctionOnContactOne_Vector_vtk(self):
726 reference="hex_3D_o2_contact_v.xml"
727 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
728 x=FunctionOnContactOne(dom).getX()
729 saveVTK("hex_contact_3D_order2_FunctionOnContactOne_Vector.xml",data=x[0]*[1.,2.,3.])
730 self.check_vtk("data.hex_contact_3D_order2_FunctionOnContactOne_Vector.xml",reference)
731 def test_hex_contact_3D_order2_FunctionOnContactOne_Tensor_vtk(self):
732 reference="hex_3D_o2_contact_t.xml"
733 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
734 x=FunctionOnContactOne(dom).getX()
735 saveVTK("hex_contact_3D_order2_FunctionOnContactOne_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
736 self.check_vtk("data.hex_contact_3D_order2_FunctionOnContactOne_Tensor.xml",reference)
737 def test_hex_contact_3D_order2_onFace_FunctionOnContactOne_Scalar_vtk(self):
738 reference="hex_3D_o2_contact_s.xml"
739 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2_onFace.msh")
740 x=FunctionOnContactOne(dom).getX()
741 saveVTK("hex_contact_3D_order2_onFace_FunctionOnContactOne_Scalar.xml",data=x[0])
742 self.check_vtk("data.hex_contact_3D_order2_onFace_FunctionOnContactOne_Scalar.xml",reference)
743 def test_hex_contact_3D_order2_onFace_FunctionOnContactOne_Vector_vtk(self):
744 reference="hex_3D_o2_contact_v.xml"
745 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2_onFace.msh")
746 x=FunctionOnContactOne(dom).getX()
747 saveVTK("hex_contact_3D_order2_onFace_FunctionOnContactOne_Vector.xml",data=x[0]*[1.,2.,3.])
748 self.check_vtk("data.hex_contact_3D_order2_onFace_FunctionOnContactOne_Vector.xml",reference)
749 def test_hex_contact_3D_order2_onFace_FunctionOnContactOne_Tensor_vtk(self):
750 reference="hex_3D_o2_contact_t.xml"
751 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2_onFace.msh")
752 x=FunctionOnContactOne(dom).getX()
753 saveVTK("hex_contact_3D_order2_onFace_FunctionOnContactOne_Tensor.xml",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
754 self.check_vtk("data.hex_contact_3D_order2_onFace_FunctionOnContactOne_Tensor.xml",reference)
755
756 class Test_DXFiles(Test_VisualizationInterface):
757 # ======================================================================================================================
758 def test_hex_contact_2D_order1_ContinuousFunction_Scalar_dx(self):
759 reference="hex_2D_o1_node_s.dx"
760 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
761 x=ContinuousFunction(dom).getX()
762 saveDX("hex_contact_2D_order1_ContinuousFunction_Scalar.dx",data=x[0])
763 self.check_dx("data.hex_contact_2D_order1_ContinuousFunction_Scalar.dx",reference)
764 def test_hex_contact_2D_order1_ContinuousFunction_Vector_dx(self):
765 reference="hex_2D_o1_node_v.dx"
766 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
767 x=ContinuousFunction(dom).getX()
768 saveDX("hex_contact_2D_order1_ContinuousFunction_Vector.dx",data=x[0]*[1.,2.])
769 self.check_dx("data.hex_contact_2D_order1_ContinuousFunction_Vector.dx",reference)
770 def test_hex_contact_2D_order1_ContinuousFunction_Tensor_dx(self):
771 reference="hex_2D_o1_node_t.dx"
772 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
773 x=ContinuousFunction(dom).getX()
774 saveDX("hex_contact_2D_order1_ContinuousFunction_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
775 self.check_dx("data.hex_contact_2D_order1_ContinuousFunction_Tensor.dx",reference)
776 def test_hex_contact_2D_order1_Solution_Scalar_dx(self):
777 reference="hex_2D_o1_node_s.dx"
778 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
779 x=Solution(dom).getX()
780 saveDX("hex_contact_2D_order1_Solution_Scalar.dx",data=x[0])
781 self.check_dx("data.hex_contact_2D_order1_Solution_Scalar.dx",reference)
782 def test_hex_contact_2D_order1_Solution_Vector_dx(self):
783 reference="hex_2D_o1_node_v.dx"
784 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
785 x=Solution(dom).getX()
786 saveDX("hex_contact_2D_order1_Solution_Vector.dx",data=x[0]*[1.,2.])
787 self.check_dx("data.hex_contact_2D_order1_Solution_Vector.dx",reference)
788 def test_hex_contact_2D_order1_Solution_Tensor_dx(self):
789 reference="hex_2D_o1_node_t.dx"
790 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
791 x=Solution(dom).getX()
792 saveDX("hex_contact_2D_order1_Solution_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
793 self.check_dx("data.hex_contact_2D_order1_Solution_Tensor.dx",reference)
794 def test_hex_contact_2D_order1_ReducedSolution_Scalar_dx(self):
795 reference="hex_2D_o1_node_s.dx"
796 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
797 x=ReducedSolution(dom).getX()
798 saveDX("hex_contact_2D_order1_ReducedSolution_Scalar.dx",data=x[0])
799 self.check_dx("data.hex_contact_2D_order1_ReducedSolution_Scalar.dx",reference)
800 def test_hex_contact_2D_order1_ReducedSolution_Vector_dx(self):
801 reference="hex_2D_o1_node_v.dx"
802 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
803 x=ReducedSolution(dom).getX()
804 saveDX("hex_contact_2D_order1_ReducedSolution_Vector.dx",data=x[0]*[1.,2.])
805 self.check_dx("data.hex_contact_2D_order1_ReducedSolution_Vector.dx",reference)
806 def test_hex_contact_2D_order1_ReducedSolution_Tensor_dx(self):
807 reference="hex_2D_o1_node_t.dx"
808 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
809 x=ReducedSolution(dom).getX()
810 saveDX("hex_contact_2D_order1_ReducedSolution_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
811 self.check_dx("data.hex_contact_2D_order1_ReducedSolution_Tensor.dx",reference)
812 def test_hex_contact_2D_order1_Function_Scalar_dx(self):
813 reference="hex_2D_o1_cell_s.dx"
814 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
815 x=Function(dom).getX()
816 saveDX("hex_contact_2D_order1_Function_Scalar.dx",data=x[0])
817 self.check_dx("data.hex_contact_2D_order1_Function_Scalar.dx",reference)
818 def test_hex_contact_2D_order1_Function_Vector_dx(self):
819 reference="hex_2D_o1_cell_v.dx"
820 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
821 x=Function(dom).getX()
822 saveDX("hex_contact_2D_order1_Function_Vector.dx",data=x[0]*[1.,2.])
823 self.check_dx("data.hex_contact_2D_order1_Function_Vector.dx",reference)
824 def test_hex_contact_2D_order1_Function_Tensor_dx(self):
825 reference="hex_2D_o1_cell_t.dx"
826 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
827 x=Function(dom).getX()
828 saveDX("hex_contact_2D_order1_Function_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
829 self.check_dx("data.hex_contact_2D_order1_Function_Tensor.dx",reference)
830 def test_hex_contact_2D_order1_FunctionOnBoundary_Scalar_dx(self):
831 reference="hex_2D_o1_boundary_s.dx"
832 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
833 x=FunctionOnBoundary(dom).getX()
834 saveDX("hex_contact_2D_order1_FunctionOnBoundary_Scalar.dx",data=x[0])
835 self.check_dx("data.hex_contact_2D_order1_FunctionOnBoundary_Scalar.dx",reference)
836 def test_hex_contact_2D_order1_FunctionOnBoundary_Vector_dx(self):
837 reference="hex_2D_o1_boundary_v.dx"
838 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
839 x=FunctionOnBoundary(dom).getX()
840 saveDX("hex_contact_2D_order1_FunctionOnBoundary_Vector.dx",data=x[0]*[1.,2.])
841 self.check_dx("data.hex_contact_2D_order1_FunctionOnBoundary_Vector.dx",reference)
842 def test_hex_contact_2D_order1_FunctionOnBoundary_Tensor_dx(self):
843 reference="hex_2D_o1_boundary_t.dx"
844 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order1.msh")
845 x=FunctionOnBoundary(dom).getX()
846 saveDX("hex_contact_2D_order1_FunctionOnBoundary_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
847 self.check_dx("data.hex_contact_2D_order1_FunctionOnBoundary_Tensor.dx",reference)
848 # ======================================================================================================================
849 def test_hex_contact_2D_order2_ContinuousFunction_Scalar_dx(self):
850 reference="hex_2D_o2_node_s.dx"
851 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
852 x=ContinuousFunction(dom).getX()
853 saveDX("hex_contact_2D_order2_ContinuousFunction_Scalar.dx",data=x[0])
854 self.check_dx("data.hex_contact_2D_order2_ContinuousFunction_Scalar.dx",reference)
855 def test_hex_contact_2D_order2_ContinuousFunction_Vector_dx(self):
856 reference="hex_2D_o2_node_v.dx"
857 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
858 x=ContinuousFunction(dom).getX()
859 saveDX("hex_contact_2D_order2_ContinuousFunction_Vector.dx",data=x[0]*[1.,2.])
860 self.check_dx("data.hex_contact_2D_order2_ContinuousFunction_Vector.dx",reference)
861 def test_hex_contact_2D_order2_ContinuousFunction_Tensor_dx(self):
862 reference="hex_2D_o2_node_t.dx"
863 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
864 x=ContinuousFunction(dom).getX()
865 saveDX("hex_contact_2D_order2_ContinuousFunction_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
866 self.check_dx("data.hex_contact_2D_order2_ContinuousFunction_Tensor.dx",reference)
867 def test_hex_contact_2D_order2_Solution_Scalar_dx(self):
868 reference="hex_2D_o2_node_s.dx"
869 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
870 x=Solution(dom).getX()
871 saveDX("hex_contact_2D_order2_Solution_Scalar.dx",data=x[0])
872 self.check_dx("data.hex_contact_2D_order2_Solution_Scalar.dx",reference)
873 def test_hex_contact_2D_order2_Solution_Vector_dx(self):
874 reference="hex_2D_o2_node_v.dx"
875 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
876 x=Solution(dom).getX()
877 saveDX("hex_contact_2D_order2_Solution_Vector.dx",data=x[0]*[1.,2.])
878 self.check_dx("data.hex_contact_2D_order2_Solution_Vector.dx",reference)
879 def test_hex_contact_2D_order2_Solution_Tensor_dx(self):
880 reference="hex_2D_o2_node_t.dx"
881 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
882 x=Solution(dom).getX()
883 saveDX("hex_contact_2D_order2_Solution_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
884 self.check_dx("data.hex_contact_2D_order2_Solution_Tensor.dx",reference)
885 def test_hex_contact_2D_order2_ReducedSolution_Scalar_dx(self):
886 reference="hex_2D_o1_node_s.dx"
887 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
888 x=ReducedSolution(dom).getX()
889 saveDX("hex_contact_2D_order2_ReducedSolution_Scalar.dx",data=x[0])
890 self.check_dx("data.hex_contact_2D_order2_ReducedSolution_Scalar.dx",reference)
891 def test_hex_contact_2D_order2_ReducedSolution_Vector_dx(self):
892 reference="hex_2D_o1_node_v.dx"
893 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
894 x=ReducedSolution(dom).getX()
895 saveDX("hex_contact_2D_order2_ReducedSolution_Vector.dx",data=x[0]*[1.,2.])
896 self.check_dx("data.hex_contact_2D_order2_ReducedSolution_Vector.dx",reference)
897 def test_hex_contact_2D_order2_ReducedSolution_Tensor_dx(self):
898 reference="hex_2D_o1_node_t.dx"
899 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
900 x=ReducedSolution(dom).getX()
901 saveDX("hex_contact_2D_order2_ReducedSolution_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
902 self.check_dx("data.hex_contact_2D_order2_ReducedSolution_Tensor.dx",reference)
903 def test_hex_contact_2D_order2_Function_Scalar_dx(self):
904 reference="hex_2D_o2_cell_s.dx"
905 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
906 x=Function(dom).getX()
907 saveDX("hex_contact_2D_order2_Function_Scalar.dx",data=x[0])
908 self.check_dx("data.hex_contact_2D_order2_Function_Scalar.dx",reference)
909 def test_hex_contact_2D_order2_Function_Vector_dx(self):
910 reference="hex_2D_o2_cell_v.dx"
911 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
912 x=Function(dom).getX()
913 saveDX("hex_contact_2D_order2_Function_Vector.dx",data=x[0]*[1.,2.])
914 self.check_dx("data.hex_contact_2D_order2_Function_Vector.dx",reference)
915 def test_hex_contact_2D_order2_Function_Tensor_dx(self):
916 reference="hex_2D_o2_cell_t.dx"
917 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
918 x=Function(dom).getX()
919 saveDX("hex_contact_2D_order2_Function_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
920 self.check_dx("data.hex_contact_2D_order2_Function_Tensor.dx",reference)
921 def test_hex_contact_2D_order2_FunctionOnBoundary_Scalar_dx(self):
922 reference="hex_2D_o2_boundary_s.dx"
923 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
924 x=FunctionOnBoundary(dom).getX()
925 saveDX("hex_contact_2D_order2_FunctionOnBoundary_Scalar.dx",data=x[0])
926 self.check_dx("data.hex_contact_2D_order2_FunctionOnBoundary_Scalar.dx",reference)
927 def test_hex_contact_2D_order2_FunctionOnBoundary_Vector_dx(self):
928 reference="hex_2D_o2_boundary_v.dx"
929 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
930 x=FunctionOnBoundary(dom).getX()
931 saveDX("hex_contact_2D_order2_FunctionOnBoundary_Vector.dx",data=x[0]*[1.,2.])
932 self.check_dx("data.hex_contact_2D_order2_FunctionOnBoundary_Vector.dx",reference)
933 def test_hex_contact_2D_order2_FunctionOnBoundary_Tensor_dx(self):
934 reference="hex_2D_o2_boundary_t.dx"
935 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_2D_order2.msh")
936 x=FunctionOnBoundary(dom).getX()
937 saveDX("hex_contact_2D_order2_FunctionOnBoundary_Tensor.dx",data=x[0]*[[11.,12.],[21.,22.]])
938 self.check_dx("data.hex_contact_2D_order2_FunctionOnBoundary_Tensor.dx",reference)
939 # ======================================================================================================================
940 def test_hex_contact_3D_order1_ContinuousFunction_Scalar_dx(self):
941 reference="hex_3D_o1_node_s.dx"
942 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
943 x=ContinuousFunction(dom).getX()
944 saveDX("hex_contact_3D_order1_ContinuousFunction_Scalar.dx",data=x[0])
945 self.check_dx("data.hex_contact_3D_order1_ContinuousFunction_Scalar.dx",reference)
946 def test_hex_contact_3D_order1_ContinuousFunction_Vector_dx(self):
947 reference="hex_3D_o1_node_v.dx"
948 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
949 x=ContinuousFunction(dom).getX()
950 saveDX("hex_contact_3D_order1_ContinuousFunction_Vector.dx",data=x[0]*[1.,2.,3.])
951 self.check_dx("data.hex_contact_3D_order1_ContinuousFunction_Vector.dx",reference)
952 def test_hex_contact_3D_order1_ContinuousFunction_Tensor_dx(self):
953 reference="hex_3D_o1_node_t.dx"
954 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
955 x=ContinuousFunction(dom).getX()
956 saveDX("hex_contact_3D_order1_ContinuousFunction_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
957 self.check_dx("data.hex_contact_3D_order1_ContinuousFunction_Tensor.dx",reference)
958 def test_hex_contact_3D_order1_Solution_Scalar_dx(self):
959 reference="hex_3D_o1_node_s.dx"
960 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
961 x=Solution(dom).getX()
962 saveDX("hex_contact_3D_order1_Solution_Scalar.dx",data=x[0])
963 self.check_dx("data.hex_contact_3D_order1_Solution_Scalar.dx",reference)
964 def test_hex_contact_3D_order1_Solution_Vector_dx(self):
965 reference="hex_3D_o1_node_v.dx"
966 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
967 x=Solution(dom).getX()
968 saveDX("hex_contact_3D_order1_Solution_Vector.dx",data=x[0]*[1.,2.,3.])
969 self.check_dx("data.hex_contact_3D_order1_Solution_Vector.dx",reference)
970 def test_hex_contact_3D_order1_Solution_Tensor_dx(self):
971 reference="hex_3D_o1_node_t.dx"
972 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
973 x=Solution(dom).getX()
974 saveDX("hex_contact_3D_order1_Solution_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
975 self.check_dx("data.hex_contact_3D_order1_Solution_Tensor.dx",reference)
976 def test_hex_contact_3D_order1_ReducedSolution_Scalar_dx(self):
977 reference="hex_3D_o1_node_s.dx"
978 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
979 x=ReducedSolution(dom).getX()
980 saveDX("hex_contact_3D_order1_ReducedSolution_Scalar.dx",data=x[0])
981 self.check_dx("data.hex_contact_3D_order1_ReducedSolution_Scalar.dx",reference)
982 def test_hex_contact_3D_order1_ReducedSolution_Vector_dx(self):
983 reference="hex_3D_o1_node_v.dx"
984 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
985 x=ReducedSolution(dom).getX()
986 saveDX("hex_contact_3D_order1_ReducedSolution_Vector.dx",data=x[0]*[1.,2.,3.])
987 self.check_dx("data.hex_contact_3D_order1_ReducedSolution_Vector.dx",reference)
988 def test_hex_contact_3D_order1_ReducedSolution_Tensor_dx(self):
989 reference="hex_3D_o1_node_t.dx"
990 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
991 x=ReducedSolution(dom).getX()
992 saveDX("hex_contact_3D_order1_ReducedSolution_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
993 self.check_dx("data.hex_contact_3D_order1_ReducedSolution_Tensor.dx",reference)
994 def test_hex_contact_3D_order1_Function_Scalar_dx(self):
995 reference="hex_3D_o1_cell_s.dx"
996 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
997 x=Function(dom).getX()
998 saveDX("hex_contact_3D_order1_Function_Scalar.dx",data=x[0])
999 self.check_dx("data.hex_contact_3D_order1_Function_Scalar.dx",reference)
1000 def test_hex_contact_3D_order1_Function_Vector_dx(self):
1001 reference="hex_3D_o1_cell_v.dx"
1002 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
1003 x=Function(dom).getX()
1004 saveDX("hex_contact_3D_order1_Function_Vector.dx",data=x[0]*[1.,2.,3.])
1005 self.check_dx("data.hex_contact_3D_order1_Function_Vector.dx",reference)
1006 def test_hex_contact_3D_order1_Function_Tensor_dx(self):
1007 reference="hex_3D_o1_cell_t.dx"
1008 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
1009 x=Function(dom).getX()
1010 saveDX("hex_contact_3D_order1_Function_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
1011 self.check_dx("data.hex_contact_3D_order1_Function_Tensor.dx",reference)
1012 def test_hex_contact_3D_order1_FunctionOnBoundary_Scalar_dx(self):
1013 reference="hex_3D_o1_boundary_s.dx"
1014 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
1015 x=FunctionOnBoundary(dom).getX()
1016 saveDX("hex_contact_3D_order1_FunctionOnBoundary_Scalar.dx",data=x[0])
1017 self.check_dx("data.hex_contact_3D_order1_FunctionOnBoundary_Scalar.dx",reference)
1018 def test_hex_contact_3D_order1_FunctionOnBoundary_Vector_dx(self):
1019 reference="hex_3D_o1_boundary_v.dx"
1020 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
1021 x=FunctionOnBoundary(dom).getX()
1022 saveDX("hex_contact_3D_order1_FunctionOnBoundary_Vector.dx",data=x[0]*[1.,2.,3.])
1023 self.check_dx("data.hex_contact_3D_order1_FunctionOnBoundary_Vector.dx",reference)
1024 def test_hex_contact_3D_order1_FunctionOnBoundary_Tensor_dx(self):
1025 reference="hex_3D_o1_boundary_t.dx"
1026 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order1.msh")
1027 x=FunctionOnBoundary(dom).getX()
1028 saveDX("hex_contact_3D_order1_FunctionOnBoundary_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
1029 self.check_dx("data.hex_contact_3D_order1_FunctionOnBoundary_Tensor.dx",reference)
1030 # ======================================================================================================================
1031 def test_hex_contact_3D_order2_ContinuousFunction_Scalar_dx(self):
1032 reference="hex_3D_o1_node_s.dx"
1033 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1034 x=ContinuousFunction(dom).getX()
1035 saveDX("hex_contact_3D_order2_ContinuousFunction_Scalar.dx",data=x[0])
1036 self.check_dx("data.hex_contact_3D_order2_ContinuousFunction_Scalar.dx",reference)
1037 def test_hex_contact_3D_order2_ContinuousFunction_Vector_dx(self):
1038 reference="hex_3D_o1_node_v.dx"
1039 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1040 x=ContinuousFunction(dom).getX()
1041 saveDX("hex_contact_3D_order2_ContinuousFunction_Vector.dx",data=x[0]*[1.,2.,3.])
1042 self.check_dx("data.hex_contact_3D_order2_ContinuousFunction_Vector.dx",reference)
1043 def test_hex_contact_3D_order2_ContinuousFunction_Tensor_dx(self):
1044 reference="hex_3D_o1_node_t.dx"
1045 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1046 x=ContinuousFunction(dom).getX()
1047 saveDX("hex_contact_3D_order2_ContinuousFunction_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
1048 self.check_dx("data.hex_contact_3D_order2_ContinuousFunction_Tensor.dx",reference)
1049 def test_hex_contact_3D_order2_Solution_Scalar_dx(self):
1050 reference="hex_3D_o1_node_s.dx"
1051 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1052 x=Solution(dom).getX()
1053 saveDX("hex_contact_3D_order2_Solution_Scalar.dx",data=x[0])
1054 self.check_dx("data.hex_contact_3D_order2_Solution_Scalar.dx",reference)
1055 def test_hex_contact_3D_order2_Solution_Vector_dx(self):
1056 reference="hex_3D_o1_node_v.dx"
1057 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1058 x=Solution(dom).getX()
1059 saveDX("hex_contact_3D_order2_Solution_Vector.dx",data=x[0]*[1.,2.,3.])
1060 self.check_dx("data.hex_contact_3D_order2_Solution_Vector.dx",reference)
1061 def test_hex_contact_3D_order2_Solution_Tensor_dx(self):
1062 reference="hex_3D_o1_node_t.dx"
1063 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1064 x=Solution(dom).getX()
1065 saveDX("hex_contact_3D_order2_Solution_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
1066 self.check_dx("data.hex_contact_3D_order2_Solution_Tensor.dx",reference)
1067 def test_hex_contact_3D_order2_ReducedSolution_Scalar_dx(self):
1068 reference="hex_3D_o1_node_s.dx"
1069 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1070 x=ReducedSolution(dom).getX()
1071 saveDX("hex_contact_3D_order2_ReducedSolution_Scalar.dx",data=x[0])
1072 self.check_dx("data.hex_contact_3D_order2_ReducedSolution_Scalar.dx",reference)
1073 def test_hex_contact_3D_order2_ReducedSolution_Vector_dx(self):
1074 reference="hex_3D_o1_node_v.dx"
1075 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1076 x=ReducedSolution(dom).getX()
1077 saveDX("hex_contact_3D_order2_ReducedSolution_Vector.dx",data=x[0]*[1.,2.,3.])
1078 self.check_dx("data.hex_contact_3D_order2_ReducedSolution_Vector.dx",reference)
1079 def test_hex_contact_3D_order2_ReducedSolution_Tensor_dx(self):
1080 reference="hex_3D_o1_node_t.dx"
1081 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1082 x=ReducedSolution(dom).getX()
1083 saveDX("hex_contact_3D_order2_ReducedSolution_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
1084 self.check_dx("data.hex_contact_3D_order2_ReducedSolution_Tensor.dx",reference)
1085 def test_hex_contact_3D_order2_Function_Scalar_dx(self):
1086 reference="hex_3D_o1_cell_s.dx"
1087 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1088 x=Function(dom).getX()
1089 saveDX("hex_contact_3D_order2_Function_Scalar.dx",data=x[0])
1090 self.check_dx("data.hex_contact_3D_order2_Function_Scalar.dx",reference)
1091 def test_hex_contact_3D_order2_Function_Vector_dx(self):
1092 reference="hex_3D_o1_cell_v.dx"
1093 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1094 x=Function(dom).getX()
1095 saveDX("hex_contact_3D_order2_Function_Vector.dx",data=x[0]*[1.,2.,3.])
1096 self.check_dx("data.hex_contact_3D_order2_Function_Vector.dx",reference)
1097 def test_hex_contact_3D_order2_Function_Tensor_dx(self):
1098 reference="hex_3D_o1_cell_t.dx"
1099 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1100 x=Function(dom).getX()
1101 saveDX("hex_contact_3D_order2_Function_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
1102 self.check_dx("data.hex_contact_3D_order2_Function_Tensor.dx",reference)
1103 def test_hex_contact_3D_order2_FunctionOnBoundary_Scalar_dx(self):
1104 reference="hex_3D_o1_boundary_s.dx"
1105 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1106 x=FunctionOnBoundary(dom).getX()
1107 saveDX("hex_contact_3D_order2_FunctionOnBoundary_Scalar.dx",data=x[0])
1108 self.check_dx("data.hex_contact_3D_order2_FunctionOnBoundary_Scalar.dx",reference)
1109 def test_hex_contact_3D_order2_FunctionOnBoundary_Vector_dx(self):
1110 reference="hex_3D_o1_boundary_v.dx"
1111 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1112 x=FunctionOnBoundary(dom).getX()
1113 saveDX("hex_contact_3D_order2_FunctionOnBoundary_Vector.dx",data=x[0]*[1.,2.,3.])
1114 self.check_dx("data.hex_contact_3D_order2_FunctionOnBoundary_Vector.dx",reference)
1115 def test_hex_contact_3D_order2_FunctionOnBoundary_Tensor_dx(self):
1116 reference="hex_3D_o1_boundary_t.dx"
1117 dom=ReadMesh(REERFENCE_FILE_DIR+"/hex_contact_3D_order2.msh")
1118 x=FunctionOnBoundary(dom).getX()
1119 saveDX("hex_contact_3D_order2_FunctionOnBoundary_Tensor.dx",data=x[0]*[[11.,12.,13.],[21.,22.,23],[31.,32.,33.]])
1120 self.check_dx("data.hex_contact_3D_order2_FunctionOnBoundary_Tensor.dx",reference)
1121
1122 if __name__ == '__main__':
1123 suite = unittest.TestSuite()
1124 suite.addTest(unittest.makeSuite(Test_VTKFiles))
1125 suite.addTest(unittest.makeSuite(Test_DXFiles))
1126 s=unittest.TextTestRunner(verbosity=2).run(suite)
1127 if s.wasSuccessful():
1128 sys.exit(0)
1129 else:
1130 sys.exit(1)

Properties

Name Value
svn:eol-style native
svn:keywords Author Date Id Revision

  ViewVC Help
Powered by ViewVC 1.1.26