8   PetscInt   i, j, dimC, dimR;  in testIdentity()  local
13   PetscCall(DMGetCoordinateDim(dm, &dimC));  in testIdentity()
16   PetscCall(DMGetWorkArray(dm, dimC * numPoints, MPIU_REAL, &mapped));  in testIdentity()
45       for (j = 0; j < dimC; j++) {  in testIdentity()
46         PetscCall(PetscPrintf(PETSC_COMM_SELF, "%+f", (double)mapped[i * dimC + j]));  in testIdentity()
47         if (j < dimC - 1) PetscCall(PetscPrintf(PETSC_COMM_SELF, ","));  in testIdentity()
71       for (j = 0; j < dimC; j++) {  in testIdentity()
72 …PetscSNPrintfCount(strBuf + offset, BUFSIZ - offset, "%+f", &count, (double)mapped[i * dimC + j]));  in testIdentity()
74         if (j < dimC - 1) {  in testIdentity()
95   PetscCall(DMRestoreWorkArray(dm, dimC * numPoints, MPIU_REAL, &mapped));  in testIdentity()
111   PetscInt    dim, dimC, isSimplex, isFE, numTests = 10;  in main()  local
124     for (dimC = dim; dimC <= PetscMin(3, dim + 1); dimC++) {  in main()
159 … PetscInt_FMT "D mesh embedded in %" PetscInt_FMT "D\n", isSimplex ? "P" : "Q", order, dim, dimC));  in main()
162           if (dimC > dim) { /* reembed in higher dimension */  in main()
179               PetscCall(PetscSectionSetDof(newSec, p, (nDof / dim) * dimC));  in main()
180               PetscCall(PetscSectionSetFieldDof(newSec, p, 0, (nDof / dim) * dimC));  in main()
190               for (j = 0; j < dim; j++) newCoordArray[i * dimC + j] = coordArray[i * dim + j];  in main()
191               for (; j < dimC; j++) newCoordArray[i * dimC + j] = 0.;  in main()
195             PetscCall(DMSetCoordinateDim(dm, dimC));  in main()
196             PetscCall(DMSetCoordinateSection(dm, dimC, newSec));  in main()
205 …PetscCall(PetscFECreateDefault(PetscObjectComm((PetscObject)dm), dim, dimC, isSimplex ? PETSC_TRUE…  in main()
210             PetscCall(DMSetCoordinateDim(coordDM, dimC));  in main()