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

Contents of /trunk/finley/test/python/run_visualization_interface.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 757 - (show annotations)
Mon Jun 26 13:12:56 2006 UTC (13 years, 9 months ago) by woo409
File MIME type: text/x-python
File size: 82726 byte(s)
+ Merge of intelc_win32 branch (revision 741:755) with trunk. Tested on iVEC altix (run_tests and py_tests all pass)

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

Properties

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

  ViewVC Help
Powered by ViewVC 1.1.26