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

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