Is there a bug in `getElementFaceNodes`?
Dear developers,
If I understand the documentation correctly, getElementFaceNodes
works like this:
node_tags = gmsh.model.mesh.getElementFaceNodes(element_type, face_type)
where face_type
is 3 for triangular faces and 4 for quadrangular faces.
A prism element (element_type=6
) have both triangular and quadrilateral faces. Thus, I thought I would be able to do something like this for a mesh containing prisms:
entities = gmsh.model.getEntities(3)
for entity in entities:
element_types, _, _ = gmsh.model.mesh.getElements(entity[0], entity[1])
for element_type in element_types:
# Only consider prisms
if element_type == 6:
for face_type in [3, 4]:
print(face_type)
node_tags = gmsh.model.mesh.getElementFaceNodes(element_type, face_type)
print("ok")
This gives me a segmentation fault using both gmsh 4.9.3 (from pip) and latest dev (pip install ... gmsh-dev).
I've attached an example with the grid from tutorials/python/t3.py
.
Thanks!