FunctionSpace on BoundaryMesh does not work in parallel

[Øyvind Evju]

from dolfin import *

mesh = UnitCubeMesh(3,3,3)

bdry = BoundaryMesh(mesh, "exterior")

V = FunctionSpace(bdry, "DG", 0)
print V.dim()

When running this code in parallel (mpirun -n 2 python foo.py) this has V.dim() == 0.

It works fine on a single process, and it also works when writing the BoundaryMesh to file in serial, and then reading it in parallel.

[Øyvind Evju]

I have found that the global indices are the same as the local indices for the boundary mesh, so the BoundaryMesh doesn't seem to recognize that it is distributed?

Any ideas on this? Any workarounds maybe?

[Jan Blechta]

With master this even fails with Signal code: Integer divide-by-zero somewhere in DofMapBuilder::compute_node_ownership.

[Prof Garth Wells]

Fix issue 14: BoundaryMesh/BoundaryComputation now works in parallel.

→ <<cset a2630ce57a0a>>

[Prof Garth Wells]

Fix issue 14: BoundaryMesh/BoundaryComputation now works in parallel.

→ <<cset e04d7beb0443>>

[Prof Garth Wells]

Issue #177 was marked as a duplicate of this issue.