I am trying to load a mesh from an xdmf file generating from the following code.
import pygmsh with pygmsh.occ.Geometry() as geom: geom.characteristic_length_min = 0.1 geom.characteristic_length_max = 0.1 rectangle1 = geom.add_rectangle([0.0, 0.0, 0.0], 1.0, 0.5) rectangle2 = geom.add_rectangle([0.0, 0.5, 0.0], 1.0, 0.5) geom.boolean_fragments(rectangle1, rectangle2) mesh = geom.generate_mesh() mesh.write("out.xdmf")
Then I use “infile” to read the mesh into fenics by below code.
from fenics import * mesh = Mesh() with XDMFFile("out.xdmf") as infile: infile.read(mesh)
But I encounter an error below:
Error: Unable to recognise cell type.
Reason: Unknown value “mixed”.
Where: This error was encountered inside XDMFFILE.cpp.
May I enquire that how can I resolve this error? Thank you.