How to use MeshValueCollection?

[rkenko]

Hi,

1 - You have to follow this procedure :

Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio

@Iactor, the problem with your code is your mesh. It is not a connected mesh, and consists of several bits and pieces. In particular, you add extra surfaces over the surfaces defined in the surface mesh. I’ve rewritten your code, using boolean fragments in gmsh to generate what I think is the desired mesh. Note that in this mesh, elements align over boundaries, which was not the case in your original code. // Gmsh project created by Jørgen S. Dokken 2020 SetFactory("OpenCASCADE"); Box(1) = {0,…

2 - If you still face errors, try to see if your mesh is well defined. Typically, with your mesh you have to be able to run this example without no errors :

Google Colaboratory

Good luck and please tell me if it solved your problem