XDMFFile::write(CellFunction) fails in parallel when ghost_mode = shared_facet

[Nathan Sime]

The following code fails when ghost mode is shared_facet and the number of MPI processes is greater than 1.

from dolfin import *

parameters["ghost_mode"] = "shared_facet"  # Toggle on/off
mesh = UnitSquareMesh(16, 16)
rank = MPI.rank(mesh.mpi_comm())
procs = CellFunction('size_t', mesh, rank)
XDMFFile("procs.xdmf").write(procs)

Sometimes the following error is produced:

HDF5-DIAG: Error detected in HDF5 (1.8.16) MPI-process 1:
  #000: ../../../src/H5Dio.c line 266 in H5Dwrite(): file selection+offset not within extent
    major: Dataspace
    minor: Out of range

Facet and Vertex-Functions work as expected.

[Jan Blechta]

Duplicate of #837.

[Nathan Sime]

whoops, thanks