@@ -233,7 +233,7 @@ This is because low level helper functions (e.g. getNodes) are doing same things
233233 point = Cases [allElements ,PointElement [__ ],2 ];
234234 line = Cases [allElements ,LineElement [__ ],2 ];
235235 surface = Cases [allElements ,TriangleElement [__ ]|QuadElement [__ ],2 ];
236- solid = Cases [allElements ,TetrahedronElement [__ ]|HexahedronElement [__ ],2 ];
236+ solid = Cases [allElements ,TetrahedronElement [__ ]|PrismElement [ __ ]| HexahedronElement [__ ],2 ];
237237
238238 (* If number of dimensions is 3 but no solid elements are specified,
239239 then we use ToBoundaryMesh to create ElementMesh. And similarly for 2 dimensions. *)
@@ -311,6 +311,8 @@ extract element topology (e.g. HexahedronElement) and its order from element typ
311311processVolume [type_ ,string_ ]:= Which [
312312 StringStartsQ [string ,"4" ],{TetrahedronElement ,4 },
313313 StringStartsQ [string ,"10" ],{TetrahedronElement ,10 },
314+ StringStartsQ [string ,"6" ],{PrismElement ,6 },
315+ StringStartsQ [string ,"15" ],{PrismElement ,15 },
314316 StringStartsQ [string ,"8" ],{HexahedronElement ,8 },
315317 StringStartsQ [string ,"20" ],{HexahedronElement ,20 },
316318 True ,Message [ImportMesh ::eltype ,type ]; Throw [$Failed ]
0 commit comments