test/python/ContactMeshCheck.py

# $Id$
"""checks the contact mesh generator for a single element by comparing the generated mesh file with a reference"""
import sys
import os
import unittest
                                                                                                                                                          
esys_root=os.getenv('ESYS_ROOT')
sys.path.append(esys_root+'/finley/lib')
sys.path.append(esys_root+'/escript/lib')
sys.path.append(esys_root+'/escript/py_src')
                                                                                                                                                          
from escript import *
from util import *
from linearPDE import *

import finley
import math
TMPFILE="tmp.msh"
numElm=1

Contact_2D_order1="""+Rectangular 2 x 2 mesh+Rectangular 2 x 2 mesh
2D-Nodes 8
0 0 11 0.000000000000000e+00 0.000000000000000e+00
1 1 12 1.000000000000000e+00 0.000000000000000e+00
2 2 21 0.000000000000000e+00 5.000000000000000e-01
3 3 22 1.000000000000000e+00 5.000000000000000e-01
4 4 11 0.000000000000000e+00 5.000000000000000e-01
5 5 12 1.000000000000000e+00 5.000000000000000e-01
6 6 21 0.000000000000000e+00 1.000000000000000e+00
7 7 22 1.000000000000000e+00 1.000000000000000e+00
Rec4 2
0 0 0 1 3 2
5 0 4 5 7 6
Line2 6
1 10 0 1
3 1 2 0
4 2 1 3
8 1 6 4
9 2 5 7
7 20 7 6
Line2_Contact 1
2 10 3 2 5 4
Point1 0
"""

Contact_2D_order1_onFace="""+Rectangular 2 x 2 mesh+Rectangular 2 x 2 mesh
2D-Nodes 8
0 0 11 0.000000000000000e+00 0.000000000000000e+00
1 1 12 1.000000000000000e+00 0.000000000000000e+00
2 2 21 0.000000000000000e+00 5.000000000000000e-01
3 3 22 1.000000000000000e+00 5.000000000000000e-01
4 4 11 0.000000000000000e+00 5.000000000000000e-01
5 5 12 1.000000000000000e+00 5.000000000000000e-01
6 6 21 0.000000000000000e+00 1.000000000000000e+00
7 7 22 1.000000000000000e+00 1.000000000000000e+00
Rec4 2
0 0 0 1 3 2
5 0 4 5 7 6
Rec4Face 6
1 10 0 1 3 2
3 1 2 0 1 3
4 2 1 3 2 0
7 20 7 6 4 5
8 1 6 4 5 7
9 2 5 7 6 4
Rec4Face_Contact 1
2 10 3 2 0 1 5 4 6 7
Point1 0
"""

Contact_2D_order2="""+Rectangular 3 x 3 mesh+Rectangular 3 x 3 mesh
2D-Nodes 16
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 2.500000000000000e-01
5 4 2 1.000000000000000e+00 2.500000000000000e-01
6 5 21 0.000000000000000e+00 5.000000000000000e-01
7 6 20 5.000000000000000e-01 5.000000000000000e-01
8 7 22 1.000000000000000e+00 5.000000000000000e-01
9 8 11 0.000000000000000e+00 5.000000000000000e-01
10 9 10 5.000000000000000e-01 5.000000000000000e-01
11 10 12 1.000000000000000e+00 5.000000000000000e-01
12 11 1 0.000000000000000e+00 7.500000000000000e-01
14 12 2 1.000000000000000e+00 7.500000000000000e-01
15 13 21 0.000000000000000e+00 1.000000000000000e+00
16 14 20 5.000000000000000e-01 1.000000000000000e+00
17 15 22 1.000000000000000e+00 1.000000000000000e+00
Rec8 2
0 0 0 2 8 6 1 5 7 3
5 0 9 11 17 15 10 14 16 12
Line3 6
1 10 0 2 1
3 1 6 0 3
4 2 2 8 5
8 1 15 9 12
9 2 11 17 14
7 20 17 15 16
Line3_Contact 1
2 10 8 6 7 11 9 10
Point1 0
"""

Contact_2D_order2_onFace="""+Rectangular 3 x 3 mesh+Rectangular 3 x 3 mesh
2D-Nodes 16
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 2.500000000000000e-01
5 4 2 1.000000000000000e+00 2.500000000000000e-01
6 5 21 0.000000000000000e+00 5.000000000000000e-01
7 6 20 5.000000000000000e-01 5.000000000000000e-01
8 7 22 1.000000000000000e+00 5.000000000000000e-01
9 8 11 0.000000000000000e+00 5.000000000000000e-01
10 9 10 5.000000000000000e-01 5.000000000000000e-01
11 10 12 1.000000000000000e+00 5.000000000000000e-01
12 11 1 0.000000000000000e+00 7.500000000000000e-01
14 12 2 1.000000000000000e+00 7.500000000000000e-01
15 13 21 0.000000000000000e+00 1.000000000000000e+00
16 14 20 5.000000000000000e-01 1.000000000000000e+00
17 15 22 1.000000000000000e+00 1.000000000000000e+00
Rec8 2
0 0 0 2 8 6 1 5 7 3
5 0 9 11 17 15 10 14 16 12
Rec8Face 6
1 10 0 2 8 6 1 5 7 3
3 1 6 0 2 8 3 1 5 7
4 2 2 8 6 0 5 7 3 1
7 20 17 15 9 11 16 12 10 14
8 1 15 9 11 17 12 10 14 16
9 2 11 17 15 9 14 16 12 10
Rec8Face_Contact 1
2 10 8 6 0 2 7 3 1 5 11 9 15 17 10 12 16 14
Point1 0
"""

Contact_3D_order1="""+Rectangular 2 x 2 x 2 mesh+Rectangular 2 x 2 x 2 mesh
3D-Nodes 16
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 5.000000000000000e-01
5 5 212 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
6 6 221 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
7 7 222 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
8 8 111 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
9 9 112 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
10 10 121 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
11 11 122 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
12 12 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
13 13 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
14 14 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
15 15 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
Hex8 2
0 0 0 1 3 2 4 5 7 6
7 0 8 9 11 10 12 13 15 14
Rec4 10
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
10 1 8 12 14 10
12 10 8 9 13 12
11 2 9 11 15 13
13 20 10 14 15 11
9 200 12 13 15 14
Rec4_Contact 1
2 100 4 5 7 6 8 9 11 10
Point1 0
"""

Contact_3D_order1_onFace="""+Rectangular 2 x 2 x 2 mesh+Rectangular 2 x 2 x 2 mesh
3D-Nodes 16
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 5.000000000000000e-01
5 5 212 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
6 6 221 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
7 7 222 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
8 8 111 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
9 9 112 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
10 10 121 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
11 11 122 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
12 12 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
13 13 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
14 14 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
15 15 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
Hex8 2
0 0 0 1 3 2 4 5 7 6
7 0 8 9 11 10 12 13 15 14
Hex8Face 10
1 100 0 2 3 1 4 6 7 5
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
9 200 12 13 15 14 8 9 11 10
10 1 8 12 14 10 9 13 15 11
11 2 9 11 15 13 8 10 14 12
12 10 8 9 13 12 10 11 15 14
13 20 10 14 15 11 8 12 13 9
Hex8Face_Contact 1
2 100 4 5 7 6 0 1 3 2 8 9 11 10 12 13 15 14
Point1 0
"""

Contact_3D_order2="""+Rectangular 3 x 3 x 3 mesh+Rectangular 3 x 3 x 3 mesh
3D-Nodes 40
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 2.500000000000000e-01
11 9 12 1.000000000000000e+00 0.000000000000000e+00 2.500000000000000e-01
15 10 21 0.000000000000000e+00 1.000000000000000e+00 2.500000000000000e-01
17 11 22 1.000000000000000e+00 1.000000000000000e+00 2.500000000000000e-01
18 12 211 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
19 13 210 5.000000000000000e-01 0.000000000000000e+00 5.000000000000000e-01
20 14 212 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
21 15 201 0.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
23 16 202 1.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
24 17 221 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
25 18 220 5.000000000000000e-01 1.000000000000000e+00 5.000000000000000e-01
26 19 222 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
27 20 111 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
28 21 110 5.000000000000000e-01 0.000000000000000e+00 5.000000000000000e-01
29 22 112 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
30 23 101 0.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
32 24 102 1.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
33 25 121 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
34 26 120 5.000000000000000e-01 1.000000000000000e+00 5.000000000000000e-01
35 27 122 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
36 28 11 0.000000000000000e+00 0.000000000000000e+00 7.500000000000000e-01
38 29 12 1.000000000000000e+00 0.000000000000000e+00 7.500000000000000e-01
42 30 21 0.000000000000000e+00 1.000000000000000e+00 7.500000000000000e-01
44 31 22 1.000000000000000e+00 1.000000000000000e+00 7.500000000000000e-01
45 32 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
46 33 210 5.000000000000000e-01 0.000000000000000e+00 1.000000000000000e+00
47 34 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
48 35 201 0.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
50 36 202 1.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
51 37 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
52 38 220 5.000000000000000e-01 1.000000000000000e+00 1.000000000000000e+00
53 39 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
Hex20 2
0 0 0 2 8 6 18 20 26 24 1 5 7 3 9 11 17 15 19 23 25 21
7 0 27 29 35 33 45 47 53 51 28 32 34 30 36 38 44 42 46 50 52 48
Rec8 10
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
10 1 27 45 51 33 36 48 42 30
12 10 27 29 47 45 28 38 46 36
11 2 29 35 53 47 32 44 50 38
13 20 33 51 53 35 42 52 44 34
9 200 45 47 53 51 46 50 52 48
Rec8_Contact 1
2 100 18 20 26 24 19 23 25 21 27 29 35 33 28 32 34 30
Point1 0
"""

Contact_3D_order2_onFace="""+Rectangular 3 x 3 x 3 mesh+Rectangular 3 x 3 x 3 mesh
3D-Nodes 40
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 2.500000000000000e-01
11 9 12 1.000000000000000e+00 0.000000000000000e+00 2.500000000000000e-01
15 10 21 0.000000000000000e+00 1.000000000000000e+00 2.500000000000000e-01
17 11 22 1.000000000000000e+00 1.000000000000000e+00 2.500000000000000e-01
18 12 211 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
19 13 210 5.000000000000000e-01 0.000000000000000e+00 5.000000000000000e-01
20 14 212 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
21 15 201 0.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
23 16 202 1.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
24 17 221 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
25 18 220 5.000000000000000e-01 1.000000000000000e+00 5.000000000000000e-01
26 19 222 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
27 20 111 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
28 21 110 5.000000000000000e-01 0.000000000000000e+00 5.000000000000000e-01
29 22 112 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
30 23 101 0.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
32 24 102 1.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
33 25 121 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
34 26 120 5.000000000000000e-01 1.000000000000000e+00 5.000000000000000e-01
35 27 122 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
36 28 11 0.000000000000000e+00 0.000000000000000e+00 7.500000000000000e-01
38 29 12 1.000000000000000e+00 0.000000000000000e+00 7.500000000000000e-01
42 30 21 0.000000000000000e+00 1.000000000000000e+00 7.500000000000000e-01
44 31 22 1.000000000000000e+00 1.000000000000000e+00 7.500000000000000e-01
45 32 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
46 33 210 5.000000000000000e-01 0.000000000000000e+00 1.000000000000000e+00
47 34 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
48 35 201 0.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
50 36 202 1.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
51 37 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
52 38 220 5.000000000000000e-01 1.000000000000000e+00 1.000000000000000e+00
53 39 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
Hex20 2
0 0 0 2 8 6 18 20 26 24 1 5 7 3 9 11 17 15 19 23 25 21
7 0 27 29 35 33 45 47 53 51 28 32 34 30 36 38 44 42 46 50 52 48
Hex20Face 10
1 100 0 6 8 2 18 24 26 20 3 7 5 1 9 15 17 11 21 25 23 19
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
9 200 45 47 53 51 27 29 35 33 46 50 52 48 36 38 44 42 28 32 34 30
10 1 27 45 51 33 29 47 53 35 36 48 42 30 28 46 52 34 38 50 44 32
11 2 29 35 53 47 27 33 51 45 32 44 50 38 28 34 52 46 30 42 48 36
12 10 27 29 47 45 33 35 53 51 28 38 46 36 30 32 50 48 34 44 52 42
13 20 33 51 53 35 27 45 47 29 42 52 44 34 30 48 50 32 36 46 38 28
Hex20Face_Contact 1
2 100 18 20 26 24 0 2 8 6 19 23 25 21 9 11 17 15 1 5 7 3 27 29 35 33 45 47 53 51 28 32 34 30 36 38 44 42 46 50 52 48
Point1 0
"""

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

def mkMesh(dim,order,onElem):
    if dim==2:
         ms1=finley.Rectangle(numElm,numElm,order,l1=0.5,useElementsOnFace=onElem)
         ms2=finley.Rectangle(numElm,numElm,order,l1=0.5,useElementsOnFace=onElem)
         ms2.setX(ms2.getX()+[0,0.5])
    else:
         ms1=finley.Brick(numElm,numElm,numElm,order,l2=0.5,useElementsOnFace=onElem)
         ms2=finley.Brick(numElm,numElm,numElm,order,l2=0.5,useElementsOnFace=onElem)
         ms2.setX(ms2.getX()+[0,0,0.5])
    return finley.JoinFaces([ms1,ms2])

failed=[]

case="Contact: 2D, order 1"
my_dom=mkMesh(2,1,FALSE)
if checker(my_dom,Contact_2D_order1):
  print case," passed."
else:
  failed.append(case)

case="Contact: 2D, order 1,elements on face"
my_dom=mkMesh(2,1,TRUE)
if checker(my_dom,Contact_2D_order1_onFace):
  print case," passed."
else:
  failed.append(case)

case="Contact: 2D, order 2"
my_dom=mkMesh(2,2,FALSE)
if checker(my_dom,Contact_2D_order2):
  print case," passed."
else:
  failed.append(case)

case="Contact: 2D, order 2,elements on face"
my_dom=mkMesh(2,2,TRUE)
if checker(my_dom,Contact_2D_order2_onFace):
  print case," passed."
else:
  failed.append(case)

case="Contact: 3D, order 1"
my_dom=mkMesh(3,1,FALSE)
if checker(my_dom,Contact_3D_order1):
  print case," passed."
else:
  failed.append(case)

case="Contact: 3D, order 1,elements on face"
my_dom=mkMesh(3,1,TRUE)
if checker(my_dom,Contact_3D_order1_onFace):
  print case," passed."
else:
  failed.append(case)

case="Contact: 3D, order 2"
my_dom=mkMesh(3,2,FALSE)
if checker(my_dom,Contact_3D_order2):
  print case," passed."
else:
  failed.append(case)

case="Contact: 3D, order 2,elements on face"
my_dom=mkMesh(3,2,TRUE)
if checker(my_dom,Contact_3D_order2_onFace):
  print case," passed."
else:
  failed.append(case)

if len(failed) == 0:
   print "@@ contact mesh generation test passed"
else:
   print "@@ contact mesh generation failed for %s"%failed
1	jgs	82	# $Id$
2			"""checks the contact mesh generator for a single element by comparing the generated mesh file with a reference"""
3			import sys
4			import os
5			import unittest
6
7			esys_root=os.getenv('ESYS_ROOT')
8			sys.path.append(esys_root+'/finley/lib')
9			sys.path.append(esys_root+'/escript/lib')
10			sys.path.append(esys_root+'/escript/py_src')
11
12			from escript import *
13			from util import *
14			from linearPDE import *
15
16			import finley
17			import math
18			TMPFILE="tmp.msh"
19			numElm=1
20
21			Contact_2D_order1="""+Rectangular 2 x 2 mesh+Rectangular 2 x 2 mesh
22			2D-Nodes 8
23			0 0 11 0.000000000000000e+00 0.000000000000000e+00
24			1 1 12 1.000000000000000e+00 0.000000000000000e+00
25			2 2 21 0.000000000000000e+00 5.000000000000000e-01
26			3 3 22 1.000000000000000e+00 5.000000000000000e-01
27			4 4 11 0.000000000000000e+00 5.000000000000000e-01
28			5 5 12 1.000000000000000e+00 5.000000000000000e-01
29			6 6 21 0.000000000000000e+00 1.000000000000000e+00
30			7 7 22 1.000000000000000e+00 1.000000000000000e+00
31			Rec4 2
32			0 0 0 1 3 2
33			5 0 4 5 7 6
34			Line2 6
35			1 10 0 1
36			3 1 2 0
37			4 2 1 3
38			8 1 6 4
39			9 2 5 7
40			7 20 7 6
41			Line2_Contact 1
42			2 10 3 2 5 4
43			Point1 0
44			"""
45
46			Contact_2D_order1_onFace="""+Rectangular 2 x 2 mesh+Rectangular 2 x 2 mesh
47			2D-Nodes 8
48			0 0 11 0.000000000000000e+00 0.000000000000000e+00
49			1 1 12 1.000000000000000e+00 0.000000000000000e+00
50			2 2 21 0.000000000000000e+00 5.000000000000000e-01
51			3 3 22 1.000000000000000e+00 5.000000000000000e-01
52			4 4 11 0.000000000000000e+00 5.000000000000000e-01
53			5 5 12 1.000000000000000e+00 5.000000000000000e-01
54			6 6 21 0.000000000000000e+00 1.000000000000000e+00
55			7 7 22 1.000000000000000e+00 1.000000000000000e+00
56			Rec4 2
57			0 0 0 1 3 2
58			5 0 4 5 7 6
59			Rec4Face 6
60			1 10 0 1 3 2
61			3 1 2 0 1 3
62			4 2 1 3 2 0
63			7 20 7 6 4 5
64			8 1 6 4 5 7
65			9 2 5 7 6 4
66			Rec4Face_Contact 1
67			2 10 3 2 0 1 5 4 6 7
68			Point1 0
69			"""
70
71			Contact_2D_order2="""+Rectangular 3 x 3 mesh+Rectangular 3 x 3 mesh
72			2D-Nodes 16
73			0 0 11 0.000000000000000e+00 0.000000000000000e+00
74			1 1 10 5.000000000000000e-01 0.000000000000000e+00
75			2 2 12 1.000000000000000e+00 0.000000000000000e+00
76			3 3 1 0.000000000000000e+00 2.500000000000000e-01
77			5 4 2 1.000000000000000e+00 2.500000000000000e-01
78			6 5 21 0.000000000000000e+00 5.000000000000000e-01
79			7 6 20 5.000000000000000e-01 5.000000000000000e-01
80			8 7 22 1.000000000000000e+00 5.000000000000000e-01
81			9 8 11 0.000000000000000e+00 5.000000000000000e-01
82			10 9 10 5.000000000000000e-01 5.000000000000000e-01
83			11 10 12 1.000000000000000e+00 5.000000000000000e-01
84			12 11 1 0.000000000000000e+00 7.500000000000000e-01
85			14 12 2 1.000000000000000e+00 7.500000000000000e-01
86			15 13 21 0.000000000000000e+00 1.000000000000000e+00
87			16 14 20 5.000000000000000e-01 1.000000000000000e+00
88			17 15 22 1.000000000000000e+00 1.000000000000000e+00
89			Rec8 2
90			0 0 0 2 8 6 1 5 7 3
91			5 0 9 11 17 15 10 14 16 12
92			Line3 6
93			1 10 0 2 1
94			3 1 6 0 3
95			4 2 2 8 5
96			8 1 15 9 12
97			9 2 11 17 14
98			7 20 17 15 16
99			Line3_Contact 1
100			2 10 8 6 7 11 9 10
101			Point1 0
102			"""
103
104			Contact_2D_order2_onFace="""+Rectangular 3 x 3 mesh+Rectangular 3 x 3 mesh
105			2D-Nodes 16
106			0 0 11 0.000000000000000e+00 0.000000000000000e+00
107			1 1 10 5.000000000000000e-01 0.000000000000000e+00
108			2 2 12 1.000000000000000e+00 0.000000000000000e+00
109			3 3 1 0.000000000000000e+00 2.500000000000000e-01
110			5 4 2 1.000000000000000e+00 2.500000000000000e-01
111			6 5 21 0.000000000000000e+00 5.000000000000000e-01
112			7 6 20 5.000000000000000e-01 5.000000000000000e-01
113			8 7 22 1.000000000000000e+00 5.000000000000000e-01
114			9 8 11 0.000000000000000e+00 5.000000000000000e-01
115			10 9 10 5.000000000000000e-01 5.000000000000000e-01
116			11 10 12 1.000000000000000e+00 5.000000000000000e-01
117			12 11 1 0.000000000000000e+00 7.500000000000000e-01
118			14 12 2 1.000000000000000e+00 7.500000000000000e-01
119			15 13 21 0.000000000000000e+00 1.000000000000000e+00
120			16 14 20 5.000000000000000e-01 1.000000000000000e+00
121			17 15 22 1.000000000000000e+00 1.000000000000000e+00
122			Rec8 2
123			0 0 0 2 8 6 1 5 7 3
124			5 0 9 11 17 15 10 14 16 12
125			Rec8Face 6
126			1 10 0 2 8 6 1 5 7 3
127			3 1 6 0 2 8 3 1 5 7
128			4 2 2 8 6 0 5 7 3 1
129			7 20 17 15 9 11 16 12 10 14
130			8 1 15 9 11 17 12 10 14 16
131			9 2 11 17 15 9 14 16 12 10
132			Rec8Face_Contact 1
133			2 10 8 6 0 2 7 3 1 5 11 9 15 17 10 12 16 14
134			Point1 0
135			"""
136
137			Contact_3D_order1="""+Rectangular 2 x 2 x 2 mesh+Rectangular 2 x 2 x 2 mesh
138			3D-Nodes 16
139			0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
140			1 1 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
141			2 2 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
142			3 3 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
143			4 4 211 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
144			5 5 212 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
145			6 6 221 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
146			7 7 222 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
147			8 8 111 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
148			9 9 112 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
149			10 10 121 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
150			11 11 122 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
151			12 12 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
152			13 13 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
153			14 14 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
154			15 15 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
155			Hex8 2
156			0 0 0 1 3 2 4 5 7 6
157			7 0 8 9 11 10 12 13 15 14
158			Rec4 10
159			1 100 0 2 3 1
160			3 1 0 4 6 2
161			5 10 0 1 5 4
162			4 2 1 3 7 5
163			6 20 2 6 7 3
164			10 1 8 12 14 10
165			12 10 8 9 13 12
166			11 2 9 11 15 13
167			13 20 10 14 15 11
168			9 200 12 13 15 14
169			Rec4_Contact 1
170			2 100 4 5 7 6 8 9 11 10
171			Point1 0
172			"""
173
174			Contact_3D_order1_onFace="""+Rectangular 2 x 2 x 2 mesh+Rectangular 2 x 2 x 2 mesh
175			3D-Nodes 16
176			0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
177			1 1 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
178			2 2 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
179			3 3 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
180			4 4 211 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
181			5 5 212 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
182			6 6 221 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
183			7 7 222 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
184			8 8 111 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
185			9 9 112 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
186			10 10 121 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
187			11 11 122 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
188			12 12 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
189			13 13 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
190			14 14 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
191			15 15 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
192			Hex8 2
193			0 0 0 1 3 2 4 5 7 6
194			7 0 8 9 11 10 12 13 15 14
195			Hex8Face 10
196			1 100 0 2 3 1 4 6 7 5
197			3 1 0 4 6 2 1 5 7 3
198			4 2 1 3 7 5 0 2 6 4
199			5 10 0 1 5 4 2 3 7 6
200			6 20 2 6 7 3 0 4 5 1
201			9 200 12 13 15 14 8 9 11 10
202			10 1 8 12 14 10 9 13 15 11
203			11 2 9 11 15 13 8 10 14 12
204			12 10 8 9 13 12 10 11 15 14
205			13 20 10 14 15 11 8 12 13 9
206			Hex8Face_Contact 1
207			2 100 4 5 7 6 0 1 3 2 8 9 11 10 12 13 15 14
208			Point1 0
209			"""
210
211			Contact_3D_order2="""+Rectangular 3 x 3 x 3 mesh+Rectangular 3 x 3 x 3 mesh
212			3D-Nodes 40
213			0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
214			1 1 110 5.000000000000000e-01 0.000000000000000e+00 0.000000000000000e+00
215			2 2 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
216			3 3 101 0.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
217			5 4 102 1.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
218			6 5 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
219			7 6 120 5.000000000000000e-01 1.000000000000000e+00 0.000000000000000e+00
220			8 7 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
221			9 8 11 0.000000000000000e+00 0.000000000000000e+00 2.500000000000000e-01
222			11 9 12 1.000000000000000e+00 0.000000000000000e+00 2.500000000000000e-01
223			15 10 21 0.000000000000000e+00 1.000000000000000e+00 2.500000000000000e-01
224			17 11 22 1.000000000000000e+00 1.000000000000000e+00 2.500000000000000e-01
225			18 12 211 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
226			19 13 210 5.000000000000000e-01 0.000000000000000e+00 5.000000000000000e-01
227			20 14 212 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
228			21 15 201 0.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
229			23 16 202 1.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
230			24 17 221 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
231			25 18 220 5.000000000000000e-01 1.000000000000000e+00 5.000000000000000e-01
232			26 19 222 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
233			27 20 111 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
234			28 21 110 5.000000000000000e-01 0.000000000000000e+00 5.000000000000000e-01
235			29 22 112 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
236			30 23 101 0.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
237			32 24 102 1.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
238			33 25 121 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
239			34 26 120 5.000000000000000e-01 1.000000000000000e+00 5.000000000000000e-01
240			35 27 122 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
241			36 28 11 0.000000000000000e+00 0.000000000000000e+00 7.500000000000000e-01
242			38 29 12 1.000000000000000e+00 0.000000000000000e+00 7.500000000000000e-01
243			42 30 21 0.000000000000000e+00 1.000000000000000e+00 7.500000000000000e-01
244			44 31 22 1.000000000000000e+00 1.000000000000000e+00 7.500000000000000e-01
245			45 32 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
246			46 33 210 5.000000000000000e-01 0.000000000000000e+00 1.000000000000000e+00
247			47 34 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
248			48 35 201 0.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
249			50 36 202 1.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
250			51 37 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
251			52 38 220 5.000000000000000e-01 1.000000000000000e+00 1.000000000000000e+00
252			53 39 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
253			Hex20 2
254			0 0 0 2 8 6 18 20 26 24 1 5 7 3 9 11 17 15 19 23 25 21
255			7 0 27 29 35 33 45 47 53 51 28 32 34 30 36 38 44 42 46 50 52 48
256			Rec8 10
257			1 100 0 6 8 2 3 7 5 1
258			3 1 0 18 24 6 9 21 15 3
259			5 10 0 2 20 18 1 11 19 9
260			4 2 2 8 26 20 5 17 23 11
261			6 20 6 24 26 8 15 25 17 7
262			10 1 27 45 51 33 36 48 42 30
263			12 10 27 29 47 45 28 38 46 36
264			11 2 29 35 53 47 32 44 50 38
265			13 20 33 51 53 35 42 52 44 34
266			9 200 45 47 53 51 46 50 52 48
267			Rec8_Contact 1
268			2 100 18 20 26 24 19 23 25 21 27 29 35 33 28 32 34 30
269			Point1 0
270			"""
271
272			Contact_3D_order2_onFace="""+Rectangular 3 x 3 x 3 mesh+Rectangular 3 x 3 x 3 mesh
273			3D-Nodes 40
274			0 0 111 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
275			1 1 110 5.000000000000000e-01 0.000000000000000e+00 0.000000000000000e+00
276			2 2 112 1.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
277			3 3 101 0.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
278			5 4 102 1.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
279			6 5 121 0.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
280			7 6 120 5.000000000000000e-01 1.000000000000000e+00 0.000000000000000e+00
281			8 7 122 1.000000000000000e+00 1.000000000000000e+00 0.000000000000000e+00
282			9 8 11 0.000000000000000e+00 0.000000000000000e+00 2.500000000000000e-01
283			11 9 12 1.000000000000000e+00 0.000000000000000e+00 2.500000000000000e-01
284			15 10 21 0.000000000000000e+00 1.000000000000000e+00 2.500000000000000e-01
285			17 11 22 1.000000000000000e+00 1.000000000000000e+00 2.500000000000000e-01
286			18 12 211 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
287			19 13 210 5.000000000000000e-01 0.000000000000000e+00 5.000000000000000e-01
288			20 14 212 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
289			21 15 201 0.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
290			23 16 202 1.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
291			24 17 221 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
292			25 18 220 5.000000000000000e-01 1.000000000000000e+00 5.000000000000000e-01
293			26 19 222 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
294			27 20 111 0.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
295			28 21 110 5.000000000000000e-01 0.000000000000000e+00 5.000000000000000e-01
296			29 22 112 1.000000000000000e+00 0.000000000000000e+00 5.000000000000000e-01
297			30 23 101 0.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
298			32 24 102 1.000000000000000e+00 5.000000000000000e-01 5.000000000000000e-01
299			33 25 121 0.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
300			34 26 120 5.000000000000000e-01 1.000000000000000e+00 5.000000000000000e-01
301			35 27 122 1.000000000000000e+00 1.000000000000000e+00 5.000000000000000e-01
302			36 28 11 0.000000000000000e+00 0.000000000000000e+00 7.500000000000000e-01
303			38 29 12 1.000000000000000e+00 0.000000000000000e+00 7.500000000000000e-01
304			42 30 21 0.000000000000000e+00 1.000000000000000e+00 7.500000000000000e-01
305			44 31 22 1.000000000000000e+00 1.000000000000000e+00 7.500000000000000e-01
306			45 32 211 0.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
307			46 33 210 5.000000000000000e-01 0.000000000000000e+00 1.000000000000000e+00
308			47 34 212 1.000000000000000e+00 0.000000000000000e+00 1.000000000000000e+00
309			48 35 201 0.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
310			50 36 202 1.000000000000000e+00 5.000000000000000e-01 1.000000000000000e+00
311			51 37 221 0.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
312			52 38 220 5.000000000000000e-01 1.000000000000000e+00 1.000000000000000e+00
313			53 39 222 1.000000000000000e+00 1.000000000000000e+00 1.000000000000000e+00
314			Hex20 2
315			0 0 0 2 8 6 18 20 26 24 1 5 7 3 9 11 17 15 19 23 25 21
316			7 0 27 29 35 33 45 47 53 51 28 32 34 30 36 38 44 42 46 50 52 48
317			Hex20Face 10
318			1 100 0 6 8 2 18 24 26 20 3 7 5 1 9 15 17 11 21 25 23 19
319			3 1 0 18 24 6 2 20 26 8 9 21 15 3 1 19 25 7 11 23 17 5
320			4 2 2 8 26 20 0 6 24 18 5 17 23 11 1 7 25 19 3 15 21 9
321			5 10 0 2 20 18 6 8 26 24 1 11 19 9 3 5 23 21 7 17 25 15
322			6 20 6 24 26 8 0 18 20 2 15 25 17 7 3 21 23 5 9 19 11 1
323			9 200 45 47 53 51 27 29 35 33 46 50 52 48 36 38 44 42 28 32 34 30
324			10 1 27 45 51 33 29 47 53 35 36 48 42 30 28 46 52 34 38 50 44 32
325			11 2 29 35 53 47 27 33 51 45 32 44 50 38 28 34 52 46 30 42 48 36
326			12 10 27 29 47 45 33 35 53 51 28 38 46 36 30 32 50 48 34 44 52 42
327			13 20 33 51 53 35 27 45 47 29 42 52 44 34 30 48 50 32 36 46 38 28
328			Hex20Face_Contact 1
329			2 100 18 20 26 24 0 2 8 6 19 23 25 21 9 11 17 15 1 5 7 3 27 29 35 33 45 47 53 51 28 32 34 30 36 38 44 42 46 50 52 48
330			Point1 0
331			"""
332
333			def checker(dom,reference):
334			dom.write(TMPFILE)
335			if reference != open(TMPFILE).read():
336			return None
337			else:
338			return not None
339
340			def mkMesh(dim,order,onElem):
341			if dim==2:
342			ms1=finley.Rectangle(numElm,numElm,order,l1=0.5,useElementsOnFace=onElem)
343			ms2=finley.Rectangle(numElm,numElm,order,l1=0.5,useElementsOnFace=onElem)
344			ms2.setX(ms2.getX()+[0,0.5])
345			else:
346			ms1=finley.Brick(numElm,numElm,numElm,order,l2=0.5,useElementsOnFace=onElem)
347			ms2=finley.Brick(numElm,numElm,numElm,order,l2=0.5,useElementsOnFace=onElem)
348			ms2.setX(ms2.getX()+[0,0,0.5])
349			return finley.JoinFaces([ms1,ms2])
350
351			failed=[]
352
353			case="Contact: 2D, order 1"
354			my_dom=mkMesh(2,1,FALSE)
355			if checker(my_dom,Contact_2D_order1):
356			print case," passed."
357			else:
358			failed.append(case)
359
360			case="Contact: 2D, order 1,elements on face"
361			my_dom=mkMesh(2,1,TRUE)
362			if checker(my_dom,Contact_2D_order1_onFace):
363			print case," passed."
364			else:
365			failed.append(case)
366
367			case="Contact: 2D, order 2"
368			my_dom=mkMesh(2,2,FALSE)
369			if checker(my_dom,Contact_2D_order2):
370			print case," passed."
371			else:
372			failed.append(case)
373
374			case="Contact: 2D, order 2,elements on face"
375			my_dom=mkMesh(2,2,TRUE)
376			if checker(my_dom,Contact_2D_order2_onFace):
377			print case," passed."
378			else:
379			failed.append(case)
380
381			case="Contact: 3D, order 1"
382			my_dom=mkMesh(3,1,FALSE)
383			if checker(my_dom,Contact_3D_order1):
384			print case," passed."
385			else:
386			failed.append(case)
387
388			case="Contact: 3D, order 1,elements on face"
389			my_dom=mkMesh(3,1,TRUE)
390			if checker(my_dom,Contact_3D_order1_onFace):
391			print case," passed."
392			else:
393			failed.append(case)
394
395			case="Contact: 3D, order 2"
396			my_dom=mkMesh(3,2,FALSE)
397			if checker(my_dom,Contact_3D_order2):
398			print case," passed."
399			else:
400			failed.append(case)
401
402			case="Contact: 3D, order 2,elements on face"
403			my_dom=mkMesh(3,2,TRUE)
404			if checker(my_dom,Contact_3D_order2_onFace):
405			print case," passed."
406			else:
407			failed.append(case)
408
409			if len(failed) == 0:
410			print "@@ contact mesh generation test passed"
411			else:
412			print "@@ contact mesh generation failed for %s"%failed
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision