XDMF output of cell-centred Functions not working properly

[Chris Richardson]

The following code produces weird output, in serial or parallel. Scalar FunctionSpaces also crash when ghost_mode is enabled.

mesh = UnitSquareMesh(20,20)
Q = VectorFunctionSpace(mesh, "DG", 0)
F = Function(Q)

E = Expression(("x[0]", "x[1]"))
F.interpolate(E)

xdmf = File("a.xdmf")
xdmf << F

[Martin Sandve Alnæs]

46798da7e9

[Prof Garth Wells]

Removing milestone: 1.5 (automated comment)