Mesh 2D gmsh to fenics legacy (2019.2)

[dokken]

A general advice is to set physical group after all synchronisation is done.

Please also note that you could have presented a way smaller example.

Also note that you cannot write both triangles and lines to the same xdmf file and be compatible with legacy dolfin, this is for instance shown here:

Error defining surface and measures for a 3D mesh imported from Gmsh

Edge integrals are not supported in 3D meshes without using MeshView. You have corrected the integral measure for the surfaces (ds, not dS) Consider the following: from dolfin import * import meshio def create_mesh(mesh, cell_type, prune_z=False): cells = mesh.get_cells_type(cell_type) cell_data = mesh.get_cell_data("gmsh:physical", cell_type) out_mesh = meshio.Mesh(points=mesh.points, cells={ cell_type: cells}, cell_data={"name_to_read": [cell_data]})…

And

Mesh generation and conversion with GMSH and PYGMSH