Boundary mesh creation fails in parallel

The boundary of the boundary of the unit cube should be empty but instead contains the boundary edges between processes.

#!/usr/bin/env mpirun -n 3 python
from dolfin import *
mesh = UnitCubeMesh(2, 2, 2)
b0 = BoundaryMesh(mesh, "exterior")
b1 = BoundaryMesh(b0, "exterior")
File("b1.pvd") << b1
assert b1.num_vertices() == 0
assert b1.num_cells() == 0

Comments (5)