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	# $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