Loading [MathJax]/extensions/tex2jax.js
Mesh Oriented datABase  (version 5.5.1)
An array-based unstructured mesh library
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
GmshUtil.cpp
Go to the documentation of this file.
1 /** 2  * MOAB, a Mesh-Oriented datABase, is a software component for creating, 3  * storing and accessing finite element mesh data. 4  * 5  * Copyright 2004 Sandia Corporation. Under the terms of Contract 6  * DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government 7  * retains certain rights in this software. 8  * 9  * This library is free software; you can redistribute it and/or 10  * modify it under the terms of the GNU Lesser General Public 11  * License as published by the Free Software Foundation; either 12  * version 2.1 of the License, or (at your option) any later version. 13  * 14  */ 15  16 #include "GmshUtil.hpp" 17  18 namespace moab 19 { 20  21 // Indexed by position in Gmsh order, containing cooresponding 22 // position in MOAB order. 23 const int hex_27_node_order[] = { 0, 1, 2, 3, 4, 5, 6, 7, // corners 24  8, 11, 12, 9, 13, 10, 14, 15, 16, 19, 17, 18, // edges 25  24, 20, 23, 21, 22, 25, // faces 26  26 }; // volume 27  28 const int tet_10_node_order[] = { 0, 1, 2, 3, 4, 5, 6, 7, 9, 8 }; 29  30 // Indexed by position in MOAB order, containing cooresponding 31 // position in Gmsh order. 32 const int pri_15_node_order[] = { 33  0, 1, 2, 3, 4, 5, // corners 34  6, 8, 9, 7, 10, 11, 12, 14, 13 // edges 35 }; 36 const int pyr_13_node_order[] = { 37  0, 1, 2, 3, 4, // corners 38  5, 8, 9, 6, 10, 7, 11, 12 // edges 39 }; 40  41 // List of GmshElemType structs, indexed by the VTK type number. 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 }, // prism w/ mid-face nodes on quads but not tris 57  { "2nd order pyramid", 14, MBMAXTYPE, 0, 0 }, // pyramid w/ mid-face nodes on quad but not tris 58  { "point", 15, MBVERTEX, 1, 0 }, // point element (0-rad sphere element?) // will be skipped 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 }, // triangle w/ 2 nodes per edge 64  { "3rd order triangle", 21, MBMAXTYPE, 0, 0 }, // " " " " " " and mid-face node 65  { "4th order triangle", 22, MBMAXTYPE, 0, 0 }, // triangle w/ 3 nodes per edge 66  { "4th order triangle", 23, MBMAXTYPE, 0, 0 }, // " " " " " " and 3 mid-face nodes 67  { "5th order triangle", 24, MBMAXTYPE, 0, 0 }, // triangle w/ 4 nodes per edge 68  { "5th order triangle", 25, MBMAXTYPE, 0, 0 }, // " " " " " " and 6 mid-face nodes 69  { "3rd order edge", 26, MBMAXTYPE, 0, 0 }, // 4-node edge 70  { "4th order edge", 27, MBMAXTYPE, 0, 0 }, // 5-node edge 71  { "5th order edge", 28, MBMAXTYPE, 0, 0 }, // 6-node edge 72  { "3rd order tetrahedron", 29, MBMAXTYPE, 0, 0 }, // tet w/ 2 nodes per edge and 1 per face 73  { "4th order tetrahedron", 30, MBMAXTYPE, 0, 0 }, // tet w/ 3 nodes per edge, 3 per face, and 1 mid-voluem 74  { "5th order tetrahedron", 31, MBMAXTYPE, 0, 0 }, // tet w/ 4 nodes per edge, 6 per face, and 4 mid-voluem 75  { 0, 32, MBMAXTYPE, 0, 0 } }; 76  77 const unsigned GmshUtil::numGmshElemType = sizeof( GmshUtil::gmshElemTypes ) / sizeof( GmshUtil::gmshElemTypes[0] ); 78  79 // Define an array, indexed by EntityType and number of nodes, 80 // containing the corresponding Gmsh element type. 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  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 88  TWENTYEIGHT_ZEROS, // MBVERTEX 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 }, // MBEDGE 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 }, // MBTRI 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 }, // MBQUAD 92  TWENTYEIGHT_ZEROS, // MBPOLYGON 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 }, // MBTET 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 }, // MBPYRAMID 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 }, // MBWEDGE 96  TWENTYEIGHT_ZEROS, // MBKNIFE 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 }, // MBHEX 98  TWENTYEIGHT_ZEROS, // MBPOLYHEDRON 99  TWENTYEIGHT_ZEROS, // MBENTITYSET 100  TWENTYEIGHT_ZEROS }; // MBMAXTYPE 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 } // namespace moab