[BUG]: Reading XMDF with prism elements

[fagnelli]

Summarize the issue

The mesh export and import using XDMF file format does not appear to work for meshes with prismatic elements. It works for all other elements.

How to reproduce the bug

I prepared a small script to illustrate the error

Minimal Example (Python)

import numpy as np

# FeniCSx
from dolfinx import mesh, io
from mpi4py import MPI

# Define mesh parameters
comm = MPI.COMM_WORLD
nel = 10                                         # number of elements per domain

Dhi = mesh.create_box(comm, [np.array([0., 0., 0.]), np.array([1., 1., 1.])],
                     [nel,nel,nel], cell_type=mesh.CellType.prism)
with io.XDMFFile(comm, "box-pri6.xdmf", "w") as xdmf:
   xdmf.write_mesh(Dhi)
   
with io.XDMFFile(comm, "box-pri6.xdmf", "r") as xdmf:
    Dhf = xdmf.read_mesh()

Output (Python)

RuntimeError: Cannot recognise cell type. Unknown value: wedge

Version

0.8.0

DOLFINx git commit

No response

Installation

I used a Docker image on a Windows 11

Additional information

No response