Empty cells refining in parallel with ghost_mode = shared_facet

[Nathan Sime]

When refining a mesh in parallel with the option

parameters['ghost_mode'] = 'shared_facet'

produces output as follows:

[Nathan Sime]

It seems the cells exist, this is an issue with XDMF output when the ghost_mode is shared_facet.

Run the following with more than one mpi process.

from dolfin import *
parameters['ghost_mode'] = 'shared_facet'

mesh = UnitSquareMesh(8, 8)
fi = XDMFFile(mesh.mpi_comm(), 'out.xdmf')
fi.write(mesh)
fi = File('out.pvd')
fi << mesh

[Anders Logg (Chalmers)]

Removing milestone: 1.7 (automated comment)

[Jan Blechta]

Does not seem to be an issue any more.