1
15
16 #include "GmshUtil.hpp"
17
18 namespace moab
19 {
20
21
22
23 const int hex_27_node_order[] = { 0, 1, 2, 3, 4, 5, 6, 7,
24 8, 11, 12, 9, 13, 10, 14, 15, 16, 19, 17, 18,
25 24, 20, 23, 21, 22, 25,
26 26 };
27
28 const int tet_10_node_order[] = { 0, 1, 2, 3, 4, 5, 6, 7, 9, 8 };
29
30
31
32 const int pri_15_node_order[] = {
33 0, 1, 2, 3, 4, 5,
34 6, 8, 9, 7, 10, 11, 12, 14, 13
35 };
36 const int pyr_13_node_order[] = {
37 0, 1, 2, 3, 4,
38 5, 8, 9, 6, 10, 7, 11, 12
39 };
40
41
42 const GmshElemType GmshUtil::gmshElemTypes[] = {
43 { 0, 0, MBMAXTYPE, 0, 0 },
44 { "line", 1, MBEDGE, 2, 0 },
45 { "triangle", 2, MBTRI, 3, 0 },
46 { "quadrangle", 3, MBQUAD, 4, 0 },
47 { "tetrahedron", 4, MBTET, 4, 0 },
48 { "hexahedron", 5, MBHEX, 8, 0 },
49 { "prism", 6, MBPRISM, 6, 0 },
50 { "pyramid", 7, MBPYRAMID, 5, 0 },
51 { "2nd order line", 8, MBEDGE, 3, 0 },
52 { "2nd order triangle", 9, MBTRI, 6, 0 },
53 { "2nd order quadrangle", 10, MBQUAD, 9, 0 },
54 { "2nd order tetrahedron", 11, MBTET, 10, tet_10_node_order },
55 { "2nd order hexahedron", 12, MBHEX, 27, hex_27_node_order },
56 { "2nd order prism", 13, MBMAXTYPE, 0, 0 },
57 { "2nd order pyramid", 14, MBMAXTYPE, 0, 0 },
58 { "point", 15, MBVERTEX, 1, 0 },
59 { "2nd order quadrangle", 16, MBQUAD, 8, 0 },
60 { "2nd order hexahedron", 17, MBHEX, 20, hex_27_node_order },
61 { "2nd order prism", 18, MBPRISM, 15, pri_15_node_order },
62 { "2nd order pyramid", 19, MBPYRAMID, 13, pyr_13_node_order },
63 { "3rd order triangle", 20, MBMAXTYPE, 0, 0 },
64 { "3rd order triangle", 21, MBMAXTYPE, 0, 0 },
65 { "4th order triangle", 22, MBMAXTYPE, 0, 0 },
66 { "4th order triangle", 23, MBMAXTYPE, 0, 0 },
67 { "5th order triangle", 24, MBMAXTYPE, 0, 0 },
68 { "5th order triangle", 25, MBMAXTYPE, 0, 0 },
69 { "3rd order edge", 26, MBMAXTYPE, 0, 0 },
70 { "4th order edge", 27, MBMAXTYPE, 0, 0 },
71 { "5th order edge", 28, MBMAXTYPE, 0, 0 },
72 { "3rd order tetrahedron", 29, MBMAXTYPE, 0, 0 },
73 { "4th order tetrahedron", 30, MBMAXTYPE, 0, 0 },
74 { "5th order tetrahedron", 31, MBMAXTYPE, 0, 0 },
75 { 0, 32, MBMAXTYPE, 0, 0 } };
76
77 const unsigned GmshUtil::numGmshElemType = sizeof( GmshUtil::gmshElemTypes ) / sizeof( GmshUtil::gmshElemTypes[0] );
78
79
80
81 #define TWENTYEIGHT_ZEROS \
82 { \
83 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 \
84 }
85 const int MAX_NODES = 28;
86 const int mb_to_gmsh_type[][MAX_NODES] = {
87
88 TWENTYEIGHT_ZEROS,
89 { 0, 0, 1, 8, 26, 27, 28, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
90 { 0, 0, 0, 2, 0, 0, 9, 0, 0, 20, 21, 0, 22, 0, 0, 23, 0, 0, 0, 0, 0, 25, 0, 0, 0, 0, 0, 0 },
91 { 0, 0, 0, 0, 3, 0, 0, 0, 16, 10, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
92 TWENTYEIGHT_ZEROS,
93 { 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 11, 0, 0, 0, 0, 0, 0, 0, 0, 0, 29, 0, 0, 0, 0, 0, 0, 0 },
94 { 0, 0, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 0, 19, 14, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
95 { 0, 0, 0, 0, 0, 0, 6, 0, 0, 0, 0, 0, 0, 0, 0, 18, 0, 0, 13, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
96 TWENTYEIGHT_ZEROS,
97 { 0, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 17, 0, 0, 0, 0, 0, 0, 12 },
98 TWENTYEIGHT_ZEROS,
99 TWENTYEIGHT_ZEROS,
100 TWENTYEIGHT_ZEROS };
101
102 int GmshUtil::get_gmsh_type( EntityType type, unsigned num_nodes )
103 {
104 if( num_nodes >= (unsigned)MAX_NODES ) return -1;
105
106 int idx = mb_to_gmsh_type[type][num_nodes];
107 if( !idx ) return -1;
108
109 return gmshElemTypes[idx].mb_type == MBMAXTYPE ? -1 : idx;
110 }
111
112 }