Fixes for Cell.geometric_dimension

NGSolve · Nov 14, 2024 · 89e8adc · 89e8adc
1 parent 8ede8bd
commit 89e8adc
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/ngsPETSc/utils/firedrake/meshes.py b/ngsPETSc/utils/firedrake/meshes.py
@@ -289,8 +289,7 @@ def createFromTopology(self, topology, name, comm):
         '''
         cell = topology.ufl_cell()
         geometric_dim = topology.topology_dm.getCoordinateDim()
-        cell = cell.reconstruct(geometric_dimension=geometric_dim)
-        element = fd.VectorElement("Lagrange", cell, 1)
+        element = fd.VectorElement("Lagrange", cell, 1, dim=geometric_dim)
         # Create mesh object
         self.firedrakeMesh = fd.MeshGeometry.__new__(fd.MeshGeometry, element, comm)
         self.firedrakeMesh._init_topology(topology)