test/python/RecMeshTest.py

# $Id$
"""checks the rectangular mesh generator for a single element by comparing the generated mesh file with a reference"""

import sys
import os
esys_root=os.getenv('ESYS_ROOT')
sys.path.append(esys_root+'/finley/lib')
sys.path.append(esys_root+'/escript/lib')
import escript
import finley

TMPFILE="tmp.msh"

ref_line1="""Rectangular mesh with 2 nodes
1D-Nodes 2
0 0 1 0.000000000000000e+00
1 1 2 1.000000000000000e+00
Line2 1
0 0 0 1
Point1 2
1 1 0
2 2 1
Point1_Contact 0
Point1 0
"""

ref_line2="""Rectangular mesh with 3 nodes
1D-Nodes 3
0 0 1 0.000000000000000e+00
1 1 0 5.000000000000000e-01
2 2 2 1.000000000000000e+00
Line3 1
0 0 0 2 1
Point1 2
1 1 0
2 2 2
Point1_Contact 0
Point1 0
"""

ref_rec1="""Rectangular 2 x 2 mesh
2D-Nodes 4
0 0 11 0.000000000000000e+00 0.000000000000000e+00
1 1 12 1.000000000000000e+00 0.000000000000000e+00
2 2 21 0.000000000000000e+00 1.000000000000000e+00
3 3 22 1.000000000000000e+00 1.000000000000000e+00
Rec4 1
0 0 0 1 3 2
Line2 4
1 10 0 1
3 1 2 0
4 2 1 3
2 20 3 2
Line2_Contact 0
Point1 0
"""

ref_rec2="""Rectangular 3 x 3 mesh
2D-Nodes 8
0 0 11 0.000000000000000e+00 0.000000000000000e+00
1 1 10 5.000000000000000e-01 0.000000000000000e+00
2 2 12 1.000000000000000e+00 0.000000000000000e+00
3 3 1 0.000000000000000e+00 5.000000000000000e-01
5 4 2 1.000000000000000e+00 5.000000000000000e-01
6 5 21 0.000000000000000e+00 1.000000000000000e+00
7 6 20 5.000000000000000e-01 1.000000000000000e+00
8 7 22 1.000000000000000e+00 1.000000000000000e+00
Rec8 1
0 0 0 2 8 6 1 5 7 3
Line3 4
1 10 0 2 1
3 1 6 0 3
4 2 2 8 5
2 20 8 6 7
Line3_Contact 0
Point1 0
"""

ref_brick1="""Rectangular 2 x 2 x 2 mesh
3D-Nodes 8
0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
1 1 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
2 2 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
3 3 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
4 4 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
5 5 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
6 6 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
7 7 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
Hex8 1
0 0 0 1 3 2 4 5 7 6
Rec4 6
1 100 0 2 3 1
3 1 0 4 6 2
5 10 0 1 5 4
4 2 1 3 7 5
6 20 2 6 7 3
2 200 4 5 7 6
Rec4_Contact 0
Point1 0
"""

ref_brick2="""Rectangular 3 x 3 x 3 mesh
3D-Nodes 20
0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
1 1 110 5.000000000000000e-01 0.000000000000000e+00 0.000000000000000e+00
2 2 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
3 3 101 0.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
5 4 102 1.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
6 5 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
7 6 120 5.000000000000000e-01 1.000000000000000e+00 0.000000000000000e+00
8 7 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
9 8 11 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
11 9 12 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
15 10 21 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
17 11 22 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
18 12 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
19 13 210 5.000000000000000e-01 0.000000000000000e+00 1.000000000000000e+00
20 14 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
21 15 201 0.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
23 16 202 1.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
24 17 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
25 18 220 5.000000000000000e-01 1.000000000000000e+00 1.000000000000000e+00
26 19 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
Hex20 1
0 0 0 2 8 6 18 20 26 24 1 5 7 3 9 11 17 15 19 23 25 21
Rec8 6
1 100 0 6 8 2 3 7 5 1
3 1 0 18 24 6 9 21 15 3
5 10 0 2 20 18 1 11 19 9
4 2 2 8 26 20 5 17 23 11
6 20 6 24 26 8 15 25 17 7
2 200 18 20 26 24 19 23 25 21
Rec8_Contact 0
Point1 0
"""
ref_line1_onElement="""Rectangular mesh with 2 nodes
1D-Nodes 2
0 0 1 0.000000000000000e+00
1 1 2 1.000000000000000e+00
Line2 1
0 0 0 1
Line2Face 2
1 1 0 1
2 2 1 0
Line2Face_Contact 0
Point1 0
"""

ref_line2_onElement="""Rectangular mesh with 3 nodes
1D-Nodes 3
0 0 1 0.000000000000000e+00
1 1 0 5.000000000000000e-01
2 2 2 1.000000000000000e+00
Line3 1
0 0 0 2 1
Line3Face 2
1 1 0 2 1
2 2 2 0 1
Line3Face_Contact 0
Point1 0
"""

ref_rec1_onElement="""Rectangular 2 x 2 mesh
2D-Nodes 4
0 0 11 0.000000000000000e+00 0.000000000000000e+00
1 1 12 1.000000000000000e+00 0.000000000000000e+00
2 2 21 0.000000000000000e+00 1.000000000000000e+00
3 3 22 1.000000000000000e+00 1.000000000000000e+00
Rec4 1
0 0 0 1 3 2
Rec4Face 4
1 10 0 1 3 2
2 20 3 2 0 1
3 1 2 0 1 3
4 2 1 3 2 0
Rec4Face_Contact 0
Point1 0
"""

ref_rec2_onElement="""Rectangular 3 x 3 mesh
2D-Nodes 8
0 0 11 0.000000000000000e+00 0.000000000000000e+00
1 1 10 5.000000000000000e-01 0.000000000000000e+00
2 2 12 1.000000000000000e+00 0.000000000000000e+00
3 3 1 0.000000000000000e+00 5.000000000000000e-01
5 4 2 1.000000000000000e+00 5.000000000000000e-01
6 5 21 0.000000000000000e+00 1.000000000000000e+00
7 6 20 5.000000000000000e-01 1.000000000000000e+00
8 7 22 1.000000000000000e+00 1.000000000000000e+00
Rec8 1
0 0 0 2 8 6 1 5 7 3
Rec8Face 4
1 10 0 2 8 6 1 5 7 3
2 20 8 6 0 2 7 3 1 5
3 1 6 0 2 8 3 1 5 7
4 2 2 8 6 0 5 7 3 1
Rec8Face_Contact 0
Point1 0
"""

ref_brick1_onElement="""Rectangular 2 x 2 x 2 mesh
3D-Nodes 8
0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
1 1 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
2 2 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
3 3 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
4 4 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
5 5 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
6 6 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
7 7 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
Hex8 1
0 0 0 1 3 2 4 5 7 6
Hex8Face 6
1 100 0 2 3 1 4 6 7 5
2 200 4 5 7 6 0 1 3 2
3 1 0 4 6 2 1 5 7 3
4 2 1 3 7 5 0 2 6 4
5 10 0 1 5 4 2 3 7 6
6 20 2 6 7 3 0 4 5 1
Hex8Face_Contact 0
Point1 0
"""

ref_brick2_onElement="""Rectangular 3 x 3 x 3 mesh
3D-Nodes 20
0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
1 1 110 5.000000000000000e-01 0.000000000000000e+00 0.000000000000000e+00
2 2 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
3 3 101 0.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
5 4 102 1.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
6 5 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
7 6 120 5.000000000000000e-01 1.000000000000000e+00 0.000000000000000e+00
8 7 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
9 8 11 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
11 9 12 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
15 10 21 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
17 11 22 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
18 12 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
19 13 210 5.000000000000000e-01 0.000000000000000e+00 1.000000000000000e+00
20 14 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
21 15 201 0.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
23 16 202 1.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
24 17 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
25 18 220 5.000000000000000e-01 1.000000000000000e+00 1.000000000000000e+00
26 19 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
Hex20 1
0 0 0 2 8 6 18 20 26 24 1 5 7 3 9 11 17 15 19 23 25 21
Hex20Face 6
1 100 0 6 8 2 18 24 26 20 3 7 5 1 9 15 17 11 21 25 23 19
2 200 18 20 26 24 0 2 8 6 19 23 25 21 9 11 17 15 1 5 7 3
3 1 0 18 24 6 2 20 26 8 9 21 15 3 1 19 25 7 11 23 17 5
4 2 2 8 26 20 0 6 24 18 5 17 23 11 1 7 25 19 3 15 21 9
5 10 0 2 20 18 6 8 26 24 1 11 19 9 3 5 23 21 7 17 25 15
6 20 6 24 26 8 0 18 20 2 15 25 17 7 3 21 23 5 9 19 11 1
Hex20Face_Contact 0
Point1 0
"""

def checker(dom,reference):
   dom.write(TMPFILE)
   mshfile=open(TMPFILE).read()
   if reference != mshfile:
      print reference
      print mshfile
      return None
   else:
      return not None

failed=[]

my_dom=finley.Interval(1,1)
if checker(my_dom,ref_line1):
  print "@@@@ line, order 1 passed"
else:
  failed.append("line, order 1")

my_dom=finley.Interval(1,2)
if checker(my_dom,ref_line2):
  print "@@@@ line, order 2 passed"
else:
  failed.append("line, order 2")

my_dom=finley.Rectangle(1,1,1)
if checker(my_dom,ref_rec1):
  print "@@@@ rec, order 1 passed"
else:
  failed.append("rec, order 1")

my_dom=finley.Rectangle(1,1,2)
if checker(my_dom,ref_rec2):
  print "@@@@ rec, order 2 passed"
else:
  failed.append("rec, order 2")

my_dom=finley.Brick(1,1,1,1)
if checker(my_dom,ref_brick1):
  print "@@@@ brick, order 1 passed"
else:
  failed.append("brick, order 1")

my_dom=finley.Brick(1,1,1,2)
if checker(my_dom,ref_brick2):
  print "@@@@ brick, order 2 passed"
else:
  failed.append("brick, order 2")

my_dom=finley.Interval(1,1,useElementsOnFace=1)
if checker(my_dom,ref_line1_onElement):
  print "@@@@ line, order 1, onElement passed"
else:
  failed.append("line, order 1, onElement")

my_dom=finley.Interval(1,2,useElementsOnFace=1)
if checker(my_dom,ref_line2_onElement):
  print "@@@@ line, order 2, onElement passed"
else:
  failed.append("line, order 2, onElement")

my_dom=finley.Rectangle(1,1,1,useElementsOnFace=1)
if checker(my_dom,ref_rec1_onElement):
  print "@@@@ rec, order 1, onElement passed"
else:
  failed.append("rec, order 1, onElement")

my_dom=finley.Rectangle(1,1,2,useElementsOnFace=1)
if checker(my_dom,ref_rec2_onElement):
  print "@@@@ rec, order 2, onElement passed"
else:
  failed.append("rec, order 2, onElement")

my_dom=finley.Brick(1,1,1,1,useElementsOnFace=1)
if checker(my_dom,ref_brick1_onElement):
  print "@@@@ brick, order 1, onElement passed"
else:
  failed.append("brick, order 1, onElement")

my_dom=finley.Brick(1,1,1,2,useElementsOnFace=1)
if checker(my_dom,ref_brick2_onElement):
  print "@@@@ brick, order 2, onElement passed"
else:
  failed.append("brick, order 2, onElement")

if len(failed) == 0:
   print "@@ mesh generation test passed"
else:
   print "@@ mesh generation failed for %s"%failed
  
1	# $Id$
2	"""checks the rectangular mesh generator for a single element by comparing the generated mesh file with a reference"""
3
4	import sys
5	import os
6	esys_root=os.getenv('ESYS_ROOT')
7	sys.path.append(esys_root+'/finley/lib')
8	sys.path.append(esys_root+'/escript/lib')
9	import escript
10	import finley
11
12	TMPFILE="tmp.msh"
13
14	ref_line1="""Rectangular mesh with 2 nodes
15	1D-Nodes 2
16	0 0 1 0.000000000000000e+00
17	1 1 2 1.000000000000000e+00
18	Line2 1
19	0 0 0 1
20	Point1 2
21	1 1 0
22	2 2 1
23	Point1_Contact 0
24	Point1 0
25	"""
26
27	ref_line2="""Rectangular mesh with 3 nodes
28	1D-Nodes 3
29	0 0 1 0.000000000000000e+00
30	1 1 0 5.000000000000000e-01
31	2 2 2 1.000000000000000e+00
32	Line3 1
33	0 0 0 2 1
34	Point1 2
35	1 1 0
36	2 2 2
37	Point1_Contact 0
38	Point1 0
39	"""
40
41	ref_rec1="""Rectangular 2 x 2 mesh
42	2D-Nodes 4
43	0 0 11 0.000000000000000e+00 0.000000000000000e+00
44	1 1 12 1.000000000000000e+00 0.000000000000000e+00
45	2 2 21 0.000000000000000e+00 1.000000000000000e+00
46	3 3 22 1.000000000000000e+00 1.000000000000000e+00
47	Rec4 1
48	0 0 0 1 3 2
49	Line2 4
50	1 10 0 1
51	3 1 2 0
52	4 2 1 3
53	2 20 3 2
54	Line2_Contact 0
55	Point1 0
56	"""
57
58	ref_rec2="""Rectangular 3 x 3 mesh
59	2D-Nodes 8
60	0 0 11 0.000000000000000e+00 0.000000000000000e+00
61	1 1 10 5.000000000000000e-01 0.000000000000000e+00
62	2 2 12 1.000000000000000e+00 0.000000000000000e+00
63	3 3 1 0.000000000000000e+00 5.000000000000000e-01
64	5 4 2 1.000000000000000e+00 5.000000000000000e-01
65	6 5 21 0.000000000000000e+00 1.000000000000000e+00
66	7 6 20 5.000000000000000e-01 1.000000000000000e+00
67	8 7 22 1.000000000000000e+00 1.000000000000000e+00
68	Rec8 1
69	0 0 0 2 8 6 1 5 7 3
70	Line3 4
71	1 10 0 2 1
72	3 1 6 0 3
73	4 2 2 8 5
74	2 20 8 6 7
75	Line3_Contact 0
76	Point1 0
77	"""
78
79	ref_brick1="""Rectangular 2 x 2 x 2 mesh
80	3D-Nodes 8
81	0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
82	1 1 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
83	2 2 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
84	3 3 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
85	4 4 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
86	5 5 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
87	6 6 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
88	7 7 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
89	Hex8 1
90	0 0 0 1 3 2 4 5 7 6
91	Rec4 6
92	1 100 0 2 3 1
93	3 1 0 4 6 2
94	5 10 0 1 5 4
95	4 2 1 3 7 5
96	6 20 2 6 7 3
97	2 200 4 5 7 6
98	Rec4_Contact 0
99	Point1 0
100	"""
101
102	ref_brick2="""Rectangular 3 x 3 x 3 mesh
103	3D-Nodes 20
104	0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
105	1 1 110 5.000000000000000e-01 0.000000000000000e+00 0.000000000000000e+00
106	2 2 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
107	3 3 101 0.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
108	5 4 102 1.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
109	6 5 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
110	7 6 120 5.000000000000000e-01 1.000000000000000e+00 0.000000000000000e+00
111	8 7 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
112	9 8 11 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
113	11 9 12 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
114	15 10 21 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
115	17 11 22 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
116	18 12 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
117	19 13 210 5.000000000000000e-01 0.000000000000000e+00 1.000000000000000e+00
118	20 14 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
119	21 15 201 0.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
120	23 16 202 1.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
121	24 17 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
122	25 18 220 5.000000000000000e-01 1.000000000000000e+00 1.000000000000000e+00
123	26 19 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
124	Hex20 1
125	0 0 0 2 8 6 18 20 26 24 1 5 7 3 9 11 17 15 19 23 25 21
126	Rec8 6
127	1 100 0 6 8 2 3 7 5 1
128	3 1 0 18 24 6 9 21 15 3
129	5 10 0 2 20 18 1 11 19 9
130	4 2 2 8 26 20 5 17 23 11
131	6 20 6 24 26 8 15 25 17 7
132	2 200 18 20 26 24 19 23 25 21
133	Rec8_Contact 0
134	Point1 0
135	"""
136	ref_line1_onElement="""Rectangular mesh with 2 nodes
137	1D-Nodes 2
138	0 0 1 0.000000000000000e+00
139	1 1 2 1.000000000000000e+00
140	Line2 1
141	0 0 0 1
142	Line2Face 2
143	1 1 0 1
144	2 2 1 0
145	Line2Face_Contact 0
146	Point1 0
147	"""
148
149	ref_line2_onElement="""Rectangular mesh with 3 nodes
150	1D-Nodes 3
151	0 0 1 0.000000000000000e+00
152	1 1 0 5.000000000000000e-01
153	2 2 2 1.000000000000000e+00
154	Line3 1
155	0 0 0 2 1
156	Line3Face 2
157	1 1 0 2 1
158	2 2 2 0 1
159	Line3Face_Contact 0
160	Point1 0
161	"""
162
163	ref_rec1_onElement="""Rectangular 2 x 2 mesh
164	2D-Nodes 4
165	0 0 11 0.000000000000000e+00 0.000000000000000e+00
166	1 1 12 1.000000000000000e+00 0.000000000000000e+00
167	2 2 21 0.000000000000000e+00 1.000000000000000e+00
168	3 3 22 1.000000000000000e+00 1.000000000000000e+00
169	Rec4 1
170	0 0 0 1 3 2
171	Rec4Face 4
172	1 10 0 1 3 2
173	2 20 3 2 0 1
174	3 1 2 0 1 3
175	4 2 1 3 2 0
176	Rec4Face_Contact 0
177	Point1 0
178	"""
179
180	ref_rec2_onElement="""Rectangular 3 x 3 mesh
181	2D-Nodes 8
182	0 0 11 0.000000000000000e+00 0.000000000000000e+00
183	1 1 10 5.000000000000000e-01 0.000000000000000e+00
184	2 2 12 1.000000000000000e+00 0.000000000000000e+00
185	3 3 1 0.000000000000000e+00 5.000000000000000e-01
186	5 4 2 1.000000000000000e+00 5.000000000000000e-01
187	6 5 21 0.000000000000000e+00 1.000000000000000e+00
188	7 6 20 5.000000000000000e-01 1.000000000000000e+00
189	8 7 22 1.000000000000000e+00 1.000000000000000e+00
190	Rec8 1
191	0 0 0 2 8 6 1 5 7 3
192	Rec8Face 4
193	1 10 0 2 8 6 1 5 7 3
194	2 20 8 6 0 2 7 3 1 5
195	3 1 6 0 2 8 3 1 5 7
196	4 2 2 8 6 0 5 7 3 1
197	Rec8Face_Contact 0
198	Point1 0
199	"""
200
201	ref_brick1_onElement="""Rectangular 2 x 2 x 2 mesh
202	3D-Nodes 8
203	0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
204	1 1 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
205	2 2 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
206	3 3 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
207	4 4 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
208	5 5 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
209	6 6 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
210	7 7 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
211	Hex8 1
212	0 0 0 1 3 2 4 5 7 6
213	Hex8Face 6
214	1 100 0 2 3 1 4 6 7 5
215	2 200 4 5 7 6 0 1 3 2
216	3 1 0 4 6 2 1 5 7 3
217	4 2 1 3 7 5 0 2 6 4
218	5 10 0 1 5 4 2 3 7 6
219	6 20 2 6 7 3 0 4 5 1
220	Hex8Face_Contact 0
221	Point1 0
222	"""
223
224	ref_brick2_onElement="""Rectangular 3 x 3 x 3 mesh
225	3D-Nodes 20
226	0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
227	1 1 110 5.000000000000000e-01 0.000000000000000e+00 0.000000000000000e+00
228	2 2 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
229	3 3 101 0.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
230	5 4 102 1.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
231	6 5 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
232	7 6 120 5.000000000000000e-01 1.000000000000000e+00 0.000000000000000e+00
233	8 7 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
234	9 8 11 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
235	11 9 12 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
236	15 10 21 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
237	17 11 22 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
238	18 12 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
239	19 13 210 5.000000000000000e-01 0.000000000000000e+00 1.000000000000000e+00
240	20 14 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
241	21 15 201 0.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
242	23 16 202 1.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
243	24 17 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
244	25 18 220 5.000000000000000e-01 1.000000000000000e+00 1.000000000000000e+00
245	26 19 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
246	Hex20 1
247	0 0 0 2 8 6 18 20 26 24 1 5 7 3 9 11 17 15 19 23 25 21
248	Hex20Face 6
249	1 100 0 6 8 2 18 24 26 20 3 7 5 1 9 15 17 11 21 25 23 19
250	2 200 18 20 26 24 0 2 8 6 19 23 25 21 9 11 17 15 1 5 7 3
251	3 1 0 18 24 6 2 20 26 8 9 21 15 3 1 19 25 7 11 23 17 5
252	4 2 2 8 26 20 0 6 24 18 5 17 23 11 1 7 25 19 3 15 21 9
253	5 10 0 2 20 18 6 8 26 24 1 11 19 9 3 5 23 21 7 17 25 15
254	6 20 6 24 26 8 0 18 20 2 15 25 17 7 3 21 23 5 9 19 11 1
255	Hex20Face_Contact 0
256	Point1 0
257	"""
258
259	def checker(dom,reference):
260	dom.write(TMPFILE)
261	mshfile=open(TMPFILE).read()
262	if reference != mshfile:
263	print reference
264	print mshfile
265	return None
266	else:
267	return not None
268
269	failed=[]
270
271	my_dom=finley.Interval(1,1)
272	if checker(my_dom,ref_line1):
273	print "@@@@ line, order 1 passed"
274	else:
275	failed.append("line, order 1")
276
277	my_dom=finley.Interval(1,2)
278	if checker(my_dom,ref_line2):
279	print "@@@@ line, order 2 passed"
280	else:
281	failed.append("line, order 2")
282
283	my_dom=finley.Rectangle(1,1,1)
284	if checker(my_dom,ref_rec1):
285	print "@@@@ rec, order 1 passed"
286	else:
287	failed.append("rec, order 1")
288
289	my_dom=finley.Rectangle(1,1,2)
290	if checker(my_dom,ref_rec2):
291	print "@@@@ rec, order 2 passed"
292	else:
293	failed.append("rec, order 2")
294
295	my_dom=finley.Brick(1,1,1,1)
296	if checker(my_dom,ref_brick1):
297	print "@@@@ brick, order 1 passed"
298	else:
299	failed.append("brick, order 1")
300
301	my_dom=finley.Brick(1,1,1,2)
302	if checker(my_dom,ref_brick2):
303	print "@@@@ brick, order 2 passed"
304	else:
305	failed.append("brick, order 2")
306
307	my_dom=finley.Interval(1,1,useElementsOnFace=1)
308	if checker(my_dom,ref_line1_onElement):
309	print "@@@@ line, order 1, onElement passed"
310	else:
311	failed.append("line, order 1, onElement")
312
313	my_dom=finley.Interval(1,2,useElementsOnFace=1)
314	if checker(my_dom,ref_line2_onElement):
315	print "@@@@ line, order 2, onElement passed"
316	else:
317	failed.append("line, order 2, onElement")
318
319	my_dom=finley.Rectangle(1,1,1,useElementsOnFace=1)
320	if checker(my_dom,ref_rec1_onElement):
321	print "@@@@ rec, order 1, onElement passed"
322	else:
323	failed.append("rec, order 1, onElement")
324
325	my_dom=finley.Rectangle(1,1,2,useElementsOnFace=1)
326	if checker(my_dom,ref_rec2_onElement):
327	print "@@@@ rec, order 2, onElement passed"
328	else:
329	failed.append("rec, order 2, onElement")
330
331	my_dom=finley.Brick(1,1,1,1,useElementsOnFace=1)
332	if checker(my_dom,ref_brick1_onElement):
333	print "@@@@ brick, order 1, onElement passed"
334	else:
335	failed.append("brick, order 1, onElement")
336
337	my_dom=finley.Brick(1,1,1,2,useElementsOnFace=1)
338	if checker(my_dom,ref_brick2_onElement):
339	print "@@@@ brick, order 2, onElement passed"
340	else:
341	failed.append("brick, order 2, onElement")
342
343	if len(failed) == 0:
344	print "@@ mesh generation test passed"
345	else:
346	print "@@ mesh generation failed for %s"%failed
347
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision