Boundary mesh creation fails in parallel
Issue #381
duplicate
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)
-
-
Issue
#394was marked as a duplicate of this issue. -
- changed milestone to 1.6
-
- changed status to duplicate
Duplicate of #313.
-
reporter - removed milestone
Removing milestone: 1.6 (automated comment)
- Log in to comment
BoundaryMesh creates all sorts of problems in parallel - we've had a lot of trouble with it for ghosted meshes.
We could remove BoundaryMesh once we have MeshViews, and then implement it properly in parallel.