Solving PDEs on 3D meshes generated with gmsh can quickly result in faulty results, when the mesh is not perfectly contructed.