Hello,
I have a problem with the script dolfin-convert when generating a relative complex geometry. I define my geometry with .geo files and gmsh. I do not have any error on gmsh when reading the geometry and the 3D meshing is OK (without error and evaluation of the mesh is OK).
On basic geometry everything is going well. However when I add surface with labels on a surface that is generated with an extrusion I have this kind of behaviour:
dolfin-convert geometry.msh geometry.xml
Converting from Gmsh format (.msh, .gmsh) to DOLFIN XML format
Expecting 13377 vertices
Found all vertices
Expecting 52853 cells
Found all cells
Traceback (most recent call last):
File "/usr/bin/dolfin-convert", line 131, in
main(sys.argv[1:])
File "/usr/bin/dolfin-convert", line 78, in main
meshconvert.convert2xml(ifilename, ofilename, iformat=iformat)
File "/usr/lib/python2.7/dist-packages/dolfin_utils/meshconvert/meshconvert.py", line 1303, in convert2xml
convert(ifilename, XmlHandler(ofilename), iformat=iformat)
File "/usr/lib/python2.7/dist-packages/dolfin_utils/meshconvert/meshconvert.py", line 1324, in convert
gmsh2xml(ifilename, handler)
File "/usr/lib/python2.7/dist-packages/dolfin_utils/meshconvert/meshconvert.py", line 494, in gmsh2xml
index = nodes_as_facets[tuple(nodes)]
KeyError: (174, 175, 2065)
I tried with less complex combinations, and everytime it fails I have the same kind of errors, all the vertices are found and all the cells are found, the .xml is produced but not finished, the physical_region.xml is produced and not finished and the facets_region.xml is not created at all.
Where this comes from? is it a bug (gmsh 2.11 produces the mesh without error nor warning) or something missing?
thank you in advance.
PS: I can send my source file if requested, I just didn't know how to include them in my post.
