Exception generated when trying to get the face indices on a specific element about pymfem HOT 5 CLOSED

@i3s93
In my environment, both work just fine. See below. Actually, I do recall that I encountered this issue before.
Reading mesh.cpp, this table is created from GetElementToFaceTable. And this method is, generally, called
when dim>2 (3D mesh). How are you creating the mesh object? Maybe there is a specific pattern in which the table
is not create.

Python 3.8.10 (default, Nov 14 2022, 12:59:47) 
[GCC 9.4.0] on linux
Type "help", "copyright", "credits" or "license" for more information.
>>> import mfem
>>> mfem.__version__
'4.5.0.2rc1'
>>> import mfem.ser as mfem
>>> mesh = mfem.Mesh('beam-tet.mesh')
>>> faces, orientations = mesh.GetElementFaces(2)
>>> faces
[7, 8, 6, 9]
>>> mesh.ElementToFaceTable()
<mfem._ser.table.Table; proxy of <Swig Object of type 'mfem::Table *' at 0x37a88b0> >

from pymfem.

Hi @sshiraiwa, thank you for your response. Actually, the mesh I used in my example is a 2-D Cartesian grid, which I create using mesh = mfem.Mesh(N_x, N_y, elem_type). As a workaround, I manually created this lookup table, which works fine, but I was hoping to just use the existing functionality. By any chance, do you recall what caused this issue before?

from pymfem.

@sshiraiwa Just to follow up, the issue only seems to happen when I use a Cartesian mesh. If I read in the mesh, as you have done, the error goes away. It looks like the table is not being generated for Cartesian meshes.

from pymfem.

Mesh::GetElementFaces is a 3D-only method. In 2D you probably want to call Mesh::GetElementEdges.

One has to be cautious -- sometimes mfem uses "faces" to refer to (dim-1)-dimensional mesh entities and sometimes to refer specifically to 2-dimensional mesh entities.

from pymfem.

@v-dobrev , Thank you for the clarification. @i3s93, closing this for now. Feel free to open if your issue is not resolved.

from pymfem.