gmsh: add simplex and tensor cells up to order 10
@@ -51,6 +51,35 @@ num_nodes_per_cell = {
'quad36': 36,
'tetra56': 56,
'hexahedron216': 216,
+ #
+ 'line7': 7,
+ 'triangle28': 28,
+ 'quad49': 49,
+ 'tetra84': 84,
+ 'hexahedron343': 343,
+ #
+ 'line8': 8,
+ 'triangle36': 36,
+ 'quad64': 64,
+ 'tetra120': 120,
+ 'hexahedron512': 512,
+ #
+ 'line9': 9,
+ 'triangle45': 45,
+ 'quad81': 81,
+ 'tetra165': 165,
+ 'hexahedron729': 729,
+ #
+ 'line10': 10,
+ 'triangle55': 55,
+ 'quad100': 100,
+ 'tetra220': 220,
+ 'hexahedron1000': 1000,
+ #
+ 'line11': 11,
+ 'triangle66': 66,
+ 'quad121': 121,
+ 'tetra286': 286,
}
# Translate meshio types to gmsh codes
@@ -85,9 +114,33 @@ _gmsh_to_meshio_type = {
36: 'quad16',
37: 'quad25',
38: 'quad36',
+ 42: 'triangle28',
+ 43: 'triangle36',
+ 44: 'triangle45',
+ 45: 'triangle55',
+ 46: 'triangle66',
+ 47: 'quad49',
+ 48: 'quad64',
+ 49: 'quad81',
+ 50: 'quad100',
+ 51: 'quad121',
+ 62: 'line7',
+ 63: 'line8',
+ 64: 'line9',
+ 65: 'line10',
+ 66: 'line11',
+ 71: 'tetra84',
+ 72: 'tetra120',
+ 73: 'tetra165',
+ 74: 'tetra220',
+ 75: 'tetra286',
92: 'hexahedron64',
93: 'hexahedron125',
94: 'hexahedron216',
+ 95: 'hexahedron343',
+ 96: 'hexahedron512',
+ 97: 'hexahedron729',
+ 98: 'hexahedron1000',
}
_meshio_to_gmsh_type = {v: k for k, v in _gmsh_to_meshio_type.items()}
