Merge pull request #12 from francesco-ballarin/francesco/dolfinx-nightly

Update FEniCSx utils due to nanobind changes
NGSolve · Nov 16, 2023 · cbf27b2 · cbf27b2
2 parents fe93e80 + e07b605
commit cbf27b2
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 3 deletions.
diff --git a/ngsPETSc/utils/fenicsx.py b/ngsPETSc/utils/fenicsx.py
@@ -34,7 +34,7 @@ def model_to_mesh(self, hmax: float, gdim: int = 2,
         dolfinx.cpp.mesh.MeshTags_int32,dolfinx.cpp.mesh.MeshTags_int32]:
         """Given a NetGen model, take all physical entities of the highest
         topological dimension and create the corresponding DOLFINx mesh.
-        
+
         This function only works in serial, at the moment.
 
         Args:
@@ -76,6 +76,6 @@ def model_to_mesh(self, hmax: float, gdim: int = 2,
         ufl_domain = dolfinx.io.gmshio.ufl_mesh(_ngs_to_cells[(gdim,T.shape[1])],gdim)
         cell_perm = dolfinx.cpp.io.perm_gmsh(dolfinx.cpp.mesh.to_type(str(ufl_domain.ufl_cell())),
                                              T.shape[1])
-        T = T[:, cell_perm]
+        T = np.ascontiguousarray(T[:, cell_perm])
         mesh = dolfinx.mesh.create_mesh(self.comm, T, V, ufl_domain, partitioner)
         return mesh
diff --git a/tests/.gitignore b/tests/.gitignore
@@ -0,0 +1 @@
+XDMF
diff --git a/tests/test_fenicsx.py b/tests/test_fenicsx.py
@@ -43,7 +43,7 @@ def test_poisson_netgen():
     geo.AddRectangle((0,0),(np.pi,np.pi))
     geoModel = ngfx.GeometricModel(geo, MPI.COMM_WORLD)
     msh  = geoModel.model_to_mesh(hmax=0.1)
-    V = fem.FunctionSpace(msh, ("Lagrange", 2)) #pylint: disable=E1120
+    V = fem.functionspace(msh, ("Lagrange", 2)) #pylint: disable=E1120
     facetsLR = mesh.locate_entities_boundary(msh, dim=(msh.topology.dim - 1),
              marker=lambda x: np.logical_or(np.isclose(x[0], 0.0),
              np.isclose(x[0], np.pi)))