I would consider moving to gmsh instead of mshr, as mshr is no longer maintained (Different behaviour of a mesh coming from generate_mesh() resp. UnitSquareMesh()? - #5 by dokken)
I think then your geometry is much cleaner, especially for coarser meshes.
1 Like