[Getfem-users] Bug in getfem::import_mesh

[Konstantinos Poulios]

Trying to import the attached gmsh file in getfem results to the
following error:

Traceback (most recent call last):
  File "tmp.py", line 23, in <module>
    m_2 = Mesh('import', 'gmsh', './TransmissionC_3D2.msh')
  File "/opt/getfem-rev3507-mumps/lib/python2.6/site-packages/getfem/getfem.py",
line 1446, in __init__
    generic_constructor(self,'mesh',*args)
  File "/opt/getfem-rev3507-mumps/lib/python2.6/site-packages/getfem/getfem.py",
line 66, in generic_constructor
    self.id = getfem_from_constructor(clname,*args)
RuntimeError: (Getfem::InterfaceError) -- Error in getfem_import.cc, line 234 :
Invalid node ID 2845 in gmsh convex 13831

It seems that the problem is related to the following object used to
map the gmsh node ids to getfem node ids:

dal::dynamic_tree_sorted<size_type> msh_node_2_getfem_node;

After I replaced this object in getfem_import.cc with:

std::map<size_type, size_type> msh_node_2_getfem_node;

the mesh can be imported without errors.

Since I have not much experience with the dal::dynamic_tree_sorted
object, every hint could be helpful.

Regards.

Kostas

TransmissionC_3D2.msh
Description: Mesh model