80   PetscInt     spaceDim, p;  in CreatePoints_Centroid()  local
87   PetscCall(DMGetCoordinateDim(dm, &spaceDim));  in CreatePoints_Centroid()
89   PetscCall(PetscCalloc1(*Np * spaceDim, pcoords));  in CreatePoints_Centroid()
95     num = size / spaceDim;  in CreatePoints_Centroid()
97 …for (d = 0; d < spaceDim; ++d) (*pcoords)[p * spaceDim + d] += PetscRealPart(coords[n * spaceDim +…  in CreatePoints_Centroid()
100     for (d = 0; d < spaceDim; ++d) {  in CreatePoints_Centroid()
101 …  PetscCall(PetscSynchronizedPrintf(PETSC_COMM_WORLD, "%g", (double)(*pcoords)[p * spaceDim + d]));  in CreatePoints_Centroid()
102       if (d < spaceDim - 1) PetscCall(PetscSynchronizedPrintf(PETSC_COMM_WORLD, ", "));  in CreatePoints_Centroid()
116   PetscInt      N = 3, n = 0, dim, spaceDim, i, j, k, *ind, d;  in CreatePoints_Grid()  local
123   PetscCall(DMGetCoordinateDim(dm, &spaceDim));  in CreatePoints_Grid()
124   PetscCall(PetscCalloc1(spaceDim, &ind));  in CreatePoints_Grid()
125   PetscCall(PetscCalloc1(spaceDim, &h));  in CreatePoints_Grid()
138   PetscCall(PetscCalloc1(*Np * spaceDim, pcoords));  in CreatePoints_Grid()
146         for (d = 0; d < spaceDim; ++d) (*pcoords)[n * spaceDim + d] = ind[d] * h[d];  in CreatePoints_Grid()
148         for (d = 0; d < spaceDim; ++d) {  in CreatePoints_Grid()
149 …  PetscCall(PetscSynchronizedPrintf(PETSC_COMM_WORLD, "%g", (double)(*pcoords)[n * spaceDim + d]));  in CreatePoints_Grid()
150           if (d < spaceDim - 1) PetscCall(PetscSynchronizedPrintf(PETSC_COMM_WORLD, ", "));  in CreatePoints_Grid()
166   PetscInt    N = 3, n = 0, dim, spaceDim, i, j, k, *ind, d;  in CreatePoints_GridReplicated()  local
173   PetscCall(DMGetCoordinateDim(dm, &spaceDim));  in CreatePoints_GridReplicated()
174   PetscCall(PetscCalloc1(spaceDim, &ind));  in CreatePoints_GridReplicated()
175   PetscCall(PetscCalloc1(spaceDim, &h));  in CreatePoints_GridReplicated()
180   PetscCall(PetscCalloc1(*Np * spaceDim, pcoords));  in CreatePoints_GridReplicated()
188         for (d = 0; d < spaceDim; ++d) (*pcoords)[n * spaceDim + d] = ind[d] * h[d];  in CreatePoints_GridReplicated()
190         for (d = 0; d < spaceDim; ++d) {  in CreatePoints_GridReplicated()
191 …  PetscCall(PetscSynchronizedPrintf(PETSC_COMM_WORLD, "%g", (double)(*pcoords)[n * spaceDim + d]));  in CreatePoints_GridReplicated()
192           if (d < spaceDim - 1) PetscCall(PetscSynchronizedPrintf(PETSC_COMM_WORLD, ", "));  in CreatePoints_GridReplicated()
268   PetscInt            dim, spaceDim, Nc, c, Np, p;  in main()  local
277   PetscCall(DMGetCoordinateDim(dm, &spaceDim));  in main()
284   PetscCall(DMInterpolationSetDim(interpolator, spaceDim));  in main()
331 … for (PetscInt i = 0; i < spaceDim; i++) vcoordsReal[i] = PetscRealPart(vcoords[p * spaceDim + i]);  in main()
333       const PetscReal *vcoordsReal = &vcoords[p * spaceDim];  in main()