Which is the preferred way to import into a FEniCS (python) program a mesh which was generated in an external generator which provides sub-domain and boundary markers? I was using MeshData->MeshFunctions for that purpose, but switching to parallel mpi program I've found that those MeshFunctions are not transferred to partial meshes used by mpi processes. MeshFunctions seem to be no longer supported - what is the alternative?
Using geometry information to recover sub-domains is extremely inconvenient for me.
I'm using my own (diffpack)grid to xml converter which is currently creating two MeshFunctions with subdomain and boundary information. Reading such xml Mesh file from a single-process code I'am able to get and use those MeshFunctions. However, when reading the same mesh from mpi code it is divided into sub-meshes and MeshFunctions are lost at this stage. My question is if may I create an xml mesh file with data which will not get lost during mesh sub-division.