Fix for legacy readers

Author: finsberg

src/io4dolfinx/readers.py

@@ -6,6 +6,7 @@
 
 from __future__ import annotations
 
+import inspect
 import pathlib
 import typing
 from pathlib import Path
@@ -149,6 +150,7 @@ def read_mesh_from_legacy_h5(
     group: str,
     cell_type: str = "tetrahedron",
     backend: str = "adios2",
+    max_facet_to_cell_links: int = 2,
 ) -> dolfinx.mesh.Mesh:
     """
     Read mesh from `h5`-file generated by legacy DOLFIN `HDF5File.write` or `XDMF.write_checkpoint`.
@@ -158,6 +160,10 @@ def read_mesh_from_legacy_h5(
         filename: Path to `h5` or `xdmf` file
         group: Name of mesh in `h5`-file
         cell_type: What type of cell type, by default tetrahedron.
+        backend: The IO backend to use when reading the mesh (must
+            support legacy mesh reading, e.g., "adios2").
+        max_facet_to_cell_links: Maximum number of cells a facet
+            can be connected to.
     """
     # Make sure we use the HDF5File and check that the file is present
     check_file_exists(filename)
@@ -175,9 +181,18 @@ def read_mesh_from_legacy_h5(
         shape=(mesh_geometry.shape[1],),
     )
     domain = ufl.Mesh(element)
+    sig = inspect.signature(dolfinx.mesh.create_mesh)
+    kwargs: dict[str, int] = {}
+    if "max_facet_to_cell_links" in list(sig.parameters.keys()):
+        kwargs["max_facet_to_cell_links"] = max_facet_to_cell_links
 
     return dolfinx.mesh.create_mesh(
-        comm=MPI.COMM_WORLD, cells=mesh_topology, x=mesh_geometry, e=domain
+        comm=MPI.COMM_WORLD,
+        cells=mesh_topology,
+        x=mesh_geometry,
+        e=domain,
+        partitioner=None,
+        **kwargs,
     )