createCoordianteSpace() test

ksagiyam · ksagiyam · commit d43d152e04db · 2025-05-08T11:17:46.000+01:00
diff --git a/firedrake/cython/dmcommon.pyx b/firedrake/cython/dmcommon.pyx
@@ -1913,6 +1913,7 @@ def _set_dg_coordinates(PETSc.DM dm,
         PetscScalar *dg_coords
         const PetscScalar *firedrake_dg_coords
         PetscInt n, gdim, cStart, cEnd, c, offset, firedrake_offset, i, j, coord_size, ndof
+        PETSc.FE coord_fe
 
     gdim = firedrake_dg_coord_vec.getBlockSize()
     coord_dm = dm.getCoordinateDM()
@@ -1945,6 +1946,7 @@ def _set_dg_coordinates(PETSc.DM dm,
                                           dg_coord_vec,
                                           True)
     dm.setCellCoordinatesLocal(dg_coord_vec)
+    dm.createCoordinateSpace(1, True, False)
 
 
 @cython.boundscheck(False)
diff --git a/firedrake/utility_meshes.py b/firedrake/utility_meshes.py
@@ -734,6 +734,7 @@ def TensorRectangleMesh(
     plex = mesh.plex_from_cell_list(
         2, cells, coords, comm, mesh._generate_default_mesh_topology_name(name)
     )
+    plex.createCoordinateSpace(1, False, False)
 
     # mark boundary facets
     plex.createLabel(dmcommon.FACE_SETS_LABEL)

Original file line number	Diff line number	Diff line change
`@@ -734,6 +734,7 @@ def TensorRectangleMesh(`
`734`	`734`	`plex = mesh.plex_from_cell_list(`
`735`	`735`	`2, cells, coords, comm, mesh._generate_default_mesh_topology_name(name)`
`736`	`736`	`)`
	`737`	`+ plex.createCoordinateSpace(1, False, False)`
`737`	`738`
`738`	`739`	`# mark boundary facets`
`739`	`740`	`plex.createLabel(dmcommon.FACE_SETS_LABEL)`