XDMFFile cannot read hex meshes
Issue #1018
resolved
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
Comments (6)
-
-
- edited description
Update MWE for dolfin master.
-
- changed status to resolved
Fix issue
#1018→ <<cset 5e5483a6727b>>
-
- changed status to open
Fix not in master yet.
-
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()
-
- changed status to resolved
Fixed in 5e5483a67
- Log in to comment