XDMFFile cannot read hex meshes

[David]

Hex meshes can be written but not read with XDMFFile.

MWE:

from dolfin import *

mesh = UnitCubeMesh.create(4, 4, 4, CellType.Type.hexahedron)
XDMFFile('mesh_hex.xdmf').write(mesh)​

mesh2 = Mesh()
XDMFFile('mesh_hex.xdmf').read(mesh2)​

Fails with the following error:

*** Error: Unable to recognise cell type.
*** Reason: Unknown value "hexahedron".
*** Where: This error was encountered inside XDMFFile.cpp.

Best wishes, David

[Jan Blechta]

Update MWE for dolfin master.

[Jan Blechta]

Fix issue #1018

→ <<cset 5e5483a6727b>>

[Jan Blechta]

Fix not in master yet.

[Jan Blechta]

Seems that the fix works:

from dolfin import *

mesh = UnitCubeMesh.create(4, 4, 4, CellType.Type.hexahedron)
XDMFFile('mesh_hex.xdmf').write(mesh)

mesh2 = Mesh()
XDMFFile('mesh_hex.xdmf').read(mesh2)

assert mesh.hash() == mesh2.hash()

[Jan Blechta]

Fixed in 5e5483a67