Different behaviour of a mesh coming from generate_mesh() resp. UnitSquareMesh()?

mesh
21

Great, thank you for that!

I am going through the thread and I have a few questions:

  1. How can I tell if my mesh has triangles or tetrahedrons? I assume it is something on the gmsh creation?
  2. I use the create_mesh function presented here. Is there any way I can get the mesh node coordinates?
  3. Why do we need two xdmf filesL the triangle and the facet (cell type=lines)?
  4. It is not clear to me how I can read the mesh and use the markers I made in gmsh in order to add boundary conditions. Later, I will need to add a sudbomain in my mesh and problem, so I think that will answer my future question too.

Many thanks!