Fail to read in parallel a mesh with marked vertices

Apparently mesh/MeshPartitioning.h:213, forbids to read mesh with a mesh value collection for vertices (dolfin_not_implemented()).

diff --git a/dolfin/mesh/MeshPartitioning.h b/dolfin/mesh/MeshPartitioning.h
index e0d2ce6..41eec8c 100644
--- a/dolfin/mesh/MeshPartitioning.h
+++ b/dolfin/mesh/MeshPartitioning.h
@@ -206,13 +206,6 @@ namespace dolfin
     // Initialise global entity numbering
     DistributedMeshTools::number_entities(mesh, dim);

-    if (dim == 0)
-    {
-      // MeshPartitioning::build_mesh_value_collection needs updating
-      // for vertices
-      dolfin_not_implemented();
-    }
-
     // Get mesh value collection used for marking
     MeshValueCollection& markers = mesh_values;

Removing the statement everything works: data structures are correct, and I'm able to build a MeshFunction (in parallel) from MeshDomains.

Comments (1)