Lines Matching refs:cdim
40 PetscInt c, cdim, i, j, o, p, vStart, vEnd; in DMPlexFindVertices() local
49 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexFindVertices()
54 …cdim == 0, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Given coordinates Vec has local length %" Petsc… in DMPlexFindVertices()
55 npoints = n / cdim; in DMPlexFindVertices()
69 …k(ndof == cdim, PETSC_COMM_SELF, PETSC_ERR_PLIB, "point %" PetscInt_FMT ": ndof = %" PetscInt_FMT … in DMPlexFindVertices()
74 for (i = 0, j = 0; i < npoints; i++, j += cdim) { in DMPlexFindVertices()
76 for (p = vStart, o = 0; p < vEnd; p++, o += cdim) { in DMPlexFindVertices()
77 for (c = 0; c < cdim; c++) { in DMPlexFindVertices()
80 if (c == cdim) { in DMPlexFindVertices()
87 for (i = 0, j = 0; i < npoints; i++, j += cdim) { in DMPlexFindVertices()
91 for (p = vStart, o = 0; p < vEnd; p++, o += cdim) { in DMPlexFindVertices()
93 for (c = 0; c < cdim; c++) norm += PetscRealPart(PetscSqr(coord[j + c] - allCoords[o + c])); in DMPlexFindVertices()
205 …GetPlaneSimplexIntersection_Coords_Internal(DM dm, PetscInt dim, PetscInt cdim, const PetscScalar … in DMPlexGetPlaneSimplexIntersection_Coords_Internal() argument
215 PetscReal mag = DMPlex_NormD_Internal(cdim, normal); in DMPlexGetPlaneSimplexIntersection_Coords_Internal()
219 dp = DMPlex_DotRealD_Internal(cdim, normal, p); in DMPlexGetPlaneSimplexIntersection_Coords_Internal()
224 for (PetscInt i = 0; i < cdim; ++i) c[i] = PetscRealPart(coords[v * cdim + i]); in DMPlexGetPlaneSimplexIntersection_Coords_Internal()
225 d[v] = DMPlex_DotRealD_Internal(cdim, normal, c); in DMPlexGetPlaneSimplexIntersection_Coords_Internal()
227 d[v] = DMPlex_DotRealD_Internal(cdim, normal, &coords[v * cdim]); in DMPlexGetPlaneSimplexIntersection_Coords_Internal()
249 … for (PetscInt i = 0; i < cdim; ++i) intPoints[n * cdim + i] = PetscRealPart(coords[v * cdim + i]); in DMPlexGetPlaneSimplexIntersection_Coords_Internal()
258 …scInt i = 0; i < cdim; ++i) intPoints[n * cdim + i] = (d[v] * PetscRealPart(coords[w * cdim + i]) … in DMPlexGetPlaneSimplexIntersection_Coords_Internal()
275 PetscInt cdim; in DMPlexGetPlaneSimplexIntersection_Internal() local
278 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexGetPlaneSimplexIntersection_Internal()
279 …PetscCheck(cdim == dim, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "DM has coordinates in %" PetscInt_F… in DMPlexGetPlaneSimplexIntersection_Internal()
284 …PetscCall(DMPlexGetPlaneSimplexIntersection_Coords_Internal(dm, dim, cdim, coords, p, normal, pos,… in DMPlexGetPlaneSimplexIntersection_Internal()
296 PetscInt cdim; in DMPlexGetPlaneQuadIntersection_Internal() local
303 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexGetPlaneQuadIntersection_Internal()
304 …PetscCheck(cdim == dim, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "DM has coordinates in %" PetscInt_F… in DMPlexGetPlaneQuadIntersection_Internal()
310 for (PetscInt d = 0; d < cdim; ++d) tcoords[v * cdim + d] = coords[vertsA[v] * cdim + d]; in DMPlexGetPlaneQuadIntersection_Internal()
311 …PetscCall(DMPlexGetPlaneSimplexIntersection_Coords_Internal(dm, dim, cdim, tcoords, p, normal, pos… in DMPlexGetPlaneQuadIntersection_Internal()
313 for (PetscInt d = 0; d < cdim; ++d) tcoords[v * cdim + d] = coords[vertsB[v] * cdim + d]; in DMPlexGetPlaneQuadIntersection_Internal()
314 …SimplexIntersection_Coords_Internal(dm, dim, cdim, tcoords, p, normal, pos, &NintB, &intPoints[Nin… in DMPlexGetPlaneQuadIntersection_Internal()
327 PetscInt cdim; in DMPlexGetPlaneHexIntersection_Internal() local
339 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexGetPlaneHexIntersection_Internal()
340 …PetscCheck(cdim == dim, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "DM has coordinates in %" PetscInt_F… in DMPlexGetPlaneHexIntersection_Internal()
346 for (PetscInt d = 0; d < cdim; ++d) tcoords[v * cdim + d] = coords[vertsA[v] * cdim + d]; in DMPlexGetPlaneHexIntersection_Internal()
347 …eSimplexIntersection_Coords_Internal(dm, dim, cdim, tcoords, p, normal, pos, &NintA, &intPoints[Ns… in DMPlexGetPlaneHexIntersection_Internal()
350 for (PetscInt d = 0; d < cdim; ++d) tcoords[v * cdim + d] = coords[vertsB[v] * cdim + d]; in DMPlexGetPlaneHexIntersection_Internal()
351 …eSimplexIntersection_Coords_Internal(dm, dim, cdim, tcoords, p, normal, pos, &NintB, &intPoints[Ns… in DMPlexGetPlaneHexIntersection_Internal()
354 for (PetscInt d = 0; d < cdim; ++d) tcoords[v * cdim + d] = coords[vertsC[v] * cdim + d]; in DMPlexGetPlaneHexIntersection_Internal()
355 …eSimplexIntersection_Coords_Internal(dm, dim, cdim, tcoords, p, normal, pos, &NintC, &intPoints[Ns… in DMPlexGetPlaneHexIntersection_Internal()
358 for (PetscInt d = 0; d < cdim; ++d) tcoords[v * cdim + d] = coords[vertsD[v] * cdim + d]; in DMPlexGetPlaneHexIntersection_Internal()
359 …eSimplexIntersection_Coords_Internal(dm, dim, cdim, tcoords, p, normal, pos, &NintD, &intPoints[Ns… in DMPlexGetPlaneHexIntersection_Internal()
362 for (PetscInt d = 0; d < cdim; ++d) tcoords[v * cdim + d] = coords[vertsE[v] * cdim + d]; in DMPlexGetPlaneHexIntersection_Internal()
363 …eSimplexIntersection_Coords_Internal(dm, dim, cdim, tcoords, p, normal, pos, &NintE, &intPoints[Ns… in DMPlexGetPlaneHexIntersection_Internal()
366 for (PetscInt d = 0; d < cdim; ++d) tcoords[v * cdim + d] = coords[vertsF[v] * cdim + d]; in DMPlexGetPlaneHexIntersection_Internal()
367 …eSimplexIntersection_Coords_Internal(dm, dim, cdim, tcoords, p, normal, pos, &NintF, &intPoints[Ns… in DMPlexGetPlaneHexIntersection_Internal()
722 PetscInt cdim, cStart, cEnd; in DMPlexCreateGridHash() local
725 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexCreateGridHash()
730 PetscCall(PetscGridHashCreate(PetscObjectComm((PetscObject)dm), cdim, a, box)); in DMPlexCreateGridHash()
739 …for (PetscInt i = 0; i < numCoords / cdim; ++i) PetscCall(PetscGridHashEnlarge(*box, &coords[i * c… in DMPlexCreateGridHash()
986 PetscInt cdim, cStart, cEnd, Nl = -1; in DMPlexComputeGridHash_Internal() local
990 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexComputeGridHash_Internal()
998 for (PetscInt i = d; i < cdim; ++i) n[i] = n[d - 1]; in DMPlexComputeGridHash_Internal()
1000 …for (PetscInt i = 0; i < cdim; ++i) n[i] = PetscMax(2, PetscFloorReal(PetscPowReal((PetscReal)(cEn… in DMPlexComputeGridHash_Internal()
1004 …t_FMT "\n h %g %g %g\n", (double)lbox->lower[0], (double)lbox->lower[1], cdim > 2 ? (double)lbox-… in DMPlexComputeGridHash_Internal()
1005 …e)lbox->upper[1], cdim > 2 ? (double)lbox->upper[2] : 0, n[0], n[1], cdim > 2 ? n[2] : 0, (double)… in DMPlexComputeGridHash_Internal()
1011 PetscCall(PetscCalloc2(16 * cdim, &dboxes, 16, &boxes)); in DMPlexComputeGridHash_Internal()
1031 PetscCall(PetscGridHashGetEnclosingBox(lbox, numCoords / cdim, coords, dboxes, boxes)); in DMPlexComputeGridHash_Internal()
1033 for (PetscInt d = 0; d < cdim; ++d) dlim[d * 2 + 0] = dlim[d * 2 + 1] = dboxes[d]; in DMPlexComputeGridHash_Internal()
1034 for (PetscInt d = cdim; d < 3; ++d) dlim[d * 2 + 0] = dlim[d * 2 + 1] = 0; in DMPlexComputeGridHash_Internal()
1035 for (PetscInt e = 1; e < numCoords / cdim; ++e) { in DMPlexComputeGridHash_Internal()
1036 for (PetscInt d = 0; d < cdim; ++d) { in DMPlexComputeGridHash_Internal()
1037 dlim[d * 2 + 0] = PetscMin(dlim[d * 2 + 0], dboxes[e * cdim + d]); in DMPlexComputeGridHash_Internal()
1038 dlim[d * 2 + 1] = PetscMax(dlim[d * 2 + 1], dboxes[e * cdim + d]); in DMPlexComputeGridHash_Internal()
1042 …for (PetscInt d = 0; d < cdim; ++d) PetscCall(PetscPrintf(PETSC_COMM_SELF, " Cell %" PetscInt_FMT… in DMPlexComputeGridHash_Internal()
1045 for (PetscInt d = 0; d < cdim; ++d) { in DMPlexComputeGridHash_Internal()
1049 for (PetscInt d = 0; d < cdim; ++d) { in DMPlexComputeGridHash_Internal()
1053 …%g, %g, %g) does not intersect %s\n", c, d, (double)lp[0], (double)lp[1], cdim > 2 ? (double)lp[2]… in DMPlexComputeGridHash_Internal()
1054 …ntersects %" PetscInt_FMT " times\n", c, d, (double)lp[0], (double)lp[1], cdim > 2 ? (double)lp[2]… in DMPlexComputeGridHash_Internal()
1059 …if (cdim > 2) PetscCall(DMPlexGetPlaneCellIntersection_Internal(dm, c, up, &normal[3 * 2], &upper[… in DMPlexComputeGridHash_Internal()
1060 if (cdim > 2 && debug > 4) { in DMPlexComputeGridHash_Internal()
1061 …2 (%g, %g, %g) does not intersect %s\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up[2]… in DMPlexComputeGridHash_Internal()
1062 …) intersects %" PetscInt_FMT " times\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up[2]… in DMPlexComputeGridHash_Internal()
1065 …if (cdim > 1) PetscCall(DMPlexGetPlaneCellIntersection_Internal(dm, c, up, &normal[3 * 1], &upper[… in DMPlexComputeGridHash_Internal()
1066 if (cdim > 1 && debug > 4) { in DMPlexComputeGridHash_Internal()
1068 …1 (%g, %g, %g) does not intersect %s\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up[2]… in DMPlexComputeGridHash_Internal()
1069 …) intersects %" PetscInt_FMT " times\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up[2]… in DMPlexComputeGridHash_Internal()
1081 …0 (%g, %g, %g) does not intersect %s\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up[2]… in DMPlexComputeGridHash_Internal()
1082 …) intersects %" PetscInt_FMT " times\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up[2]… in DMPlexComputeGridHash_Internal()
1084 for (PetscInt d = 0; d < cdim; ++d) { in DMPlexComputeGridHash_Internal()
1091 for (PetscInt d = 0; d < cdim; ++d) in DMPlexComputeGridHash_Internal()
1097 for (PetscInt d = 0; d < cdim; ++d) in DMPlexComputeGridHash_Internal()
1122 for (PetscInt plane = 0; plane < cdim; ++plane) { in DMPlexComputeGridHash_Internal()
1126 for (d = 0; d < cdim; ++d) { in DMPlexComputeGridHash_Internal()
1127 …if ((lowerIntPoints[plane][ip * cdim + d] < (lp[d] - PETSC_SMALL)) || (lowerIntPoints[plane][ip * … in DMPlexComputeGridHash_Internal()
1128 …if (lowerIntPoints[plane][ip * cdim + d] < (lp[d] - PETSC_SMALL)) outsideCount[d]++; // The lower … in DMPlexComputeGridHash_Internal()
1132 if (d == cdim) { in DMPlexComputeGridHash_Internal()
1141 for (d = 0; d < cdim; ++d) { in DMPlexComputeGridHash_Internal()
1142 …if ((upperIntPoints[plane][ip * cdim + d] < (lp[d] - PETSC_SMALL)) || (upperIntPoints[plane][ip * … in DMPlexComputeGridHash_Internal()
1143 …if (upperIntPoints[plane][ip * cdim + d] > (up[d] + PETSC_SMALL)) outsideCount[cdim + d]++; // The… in DMPlexComputeGridHash_Internal()
1147 if (d == cdim) { in DMPlexComputeGridHash_Internal()
1161 if (cdim == 2) { in DMPlexComputeGridHash_Internal()
1163 … for (PetscInt d = 0; d < cdim; ++d) nIntersects += (outsideCount[d] + outsideCount[d + cdim]); in DMPlexComputeGridHash_Internal()
1175 if (cdim == 3) { in DMPlexComputeGridHash_Internal()
1179 for (PetscInt d = 0; d < cdim; ++d) in DMPlexComputeGridHash_Internal()
1180 if (outsideCount[d] >= 3 && outsideCount[cdim + d] >= 3) { in DMPlexComputeGridHash_Internal()
1195 for (PetscInt plane = 0; plane < cdim; ++plane) { in DMPlexComputeGridHash_Internal()
1199 …werIntPoints[plane][ip * cdim + nonIntersectingFace] > lp[nonIntersectingFace] - PETSC_SMALL || lo… in DMPlexComputeGridHash_Internal()
1202 …perIntPoints[plane][ip * cdim + nonIntersectingFace] > lp[nonIntersectingFace] - PETSC_SMALL || up… in DMPlexComputeGridHash_Internal()
1251 PetscInt cdim, Nl = 0, cStart, cEnd, numCells; in DMLocatePoints_Plex() local
1268 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMLocatePoints_Plex()
1273 …cdim, PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_WRONG, "Block size for point vector %" Petsc… in DMLocatePoints_Plex()
1319 for (PetscInt d = 0; d < cdim; d++) { in DMLocatePoints_Plex()
1341 PetscCall(DMPlexLocatePoint_Internal(dm, cdim, point, c, &cell)); in DMLocatePoints_Plex()
1353 …ank, p, (double)PetscRealPart(point[0]), (double)PetscRealPart(point[1]), cdim > 2 ? (double)Petsc… in DMLocatePoints_Plex()
1360 …", %" PetscInt_FMT ", %" PetscInt_FMT ")\n", rank, bin, dbin[0], dbin[1], cdim > 2 ? dbin[2] : 0)); in DMLocatePoints_Plex()
1366 PetscCall(DMPlexLocatePoint_Internal(dm, cdim, point, boxCells[c], &cell)); in DMLocatePoints_Plex()
1384 PetscCall(DMPlexLocatePoint_Internal(dm, cdim, point, c, &cell)); in DMLocatePoints_Plex()
1409 PetscCall(DMPlexClosestPoint_Internal(dm, cdim, point, boxCells[c], cpoint)); in DMLocatePoints_Plex()
1410 for (PetscInt d = 0; d < cdim; ++d) diff[d] = cpoint[d] - PetscRealPart(point[d]); in DMLocatePoints_Plex()
1411 dist = DMPlex_NormD_Internal(cdim, diff); in DMLocatePoints_Plex()
1413 for (PetscInt d = 0; d < cdim; ++d) best[d] = cpoint[d]; in DMLocatePoints_Plex()
1422 for (PetscInt d = 0; d < cdim; ++d) a[p * bs + d] = best[d]; in DMLocatePoints_Plex()
2498 PetscInt dim, cdim, pdim, qdim, Nq, q; in DMPlexComputeCellGeometryFEM_FE() local
2502 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexComputeCellGeometryFEM_FE()
2518 …cdim, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "There are %" PetscInt_FMT " coordinates for point %" Pe… in DMPlexComputeCellGeometryFEM_FE()
2530 …PetscAssert(cdim == T->Nc, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Nc %" PetscInt_FMT " != %" PetscIn… in DMPlexComputeCellGeometryFEM_FE()
2531 …PetscAssert(dim == T->cdim, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "cdim %" PetscInt_FMT " != %" Pets… in DMPlexComputeCellGeometryFEM_FE()
2533 PetscCall(PetscArrayzero(v, Nq * cdim)); in DMPlexComputeCellGeometryFEM_FE()
2538 const PetscInt vertex = k / cdim; in DMPlexComputeCellGeometryFEM_FE()
2539 …for (i = 0; i < cdim; ++i) v[q * cdim + i] += basis[(q * pdim + k) * cdim + i] * PetscRealPart(coo… in DMPlexComputeCellGeometryFEM_FE()
2541 PetscCall(PetscLogFlops(2.0 * pdim * cdim)); in DMPlexComputeCellGeometryFEM_FE()
2545 PetscCall(PetscArrayzero(J, Nq * cdim * cdim)); in DMPlexComputeCellGeometryFEM_FE()
2551 const PetscInt vertex = k / cdim; in DMPlexComputeCellGeometryFEM_FE()
2553 … (i = 0; i < cdim; ++i) J[(q * cdim + i) * cdim + j] += basisDer[((q * pdim + k) * cdim + i) * dim… in DMPlexComputeCellGeometryFEM_FE()
2556 PetscCall(PetscLogFlops(2.0 * pdim * dim * cdim)); in DMPlexComputeCellGeometryFEM_FE()
2557 if (cdim > dim) { in DMPlexComputeCellGeometryFEM_FE()
2558 for (c = dim; c < cdim; ++c) in DMPlexComputeCellGeometryFEM_FE()
2559 for (r = 0; r < cdim; ++r) J[r * cdim + c] = r == c ? 1.0 : 0.0; in DMPlexComputeCellGeometryFEM_FE()
2563 switch (cdim) { in DMPlexComputeCellGeometryFEM_FE()
2565 DMPlex_Det3D_Internal(&detJt, &J[q * cdim * dim]); in DMPlexComputeCellGeometryFEM_FE()
2566 if (invJ) DMPlex_Invert3D_Internal(&invJ[q * cdim * dim], &J[q * cdim * dim], detJt); in DMPlexComputeCellGeometryFEM_FE()
2569 DMPlex_Det2D_Internal(&detJt, &J[q * cdim * dim]); in DMPlexComputeCellGeometryFEM_FE()
2570 if (invJ) DMPlex_Invert2D_Internal(&invJ[q * cdim * dim], &J[q * cdim * dim], detJt); in DMPlexComputeCellGeometryFEM_FE()
2573 detJt = J[q * cdim * dim]; in DMPlexComputeCellGeometryFEM_FE()
2574 if (invJ) invJ[q * cdim * dim] = 1.0 / detJt; in DMPlexComputeCellGeometryFEM_FE()
2616 PetscInt dim, cdim; in DMPlexComputeCellGeometryFEM() local
2621 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexComputeCellGeometryFEM()
2637 …if (!fe || (dim != cdim)) PetscCall(DMPlexComputeCellGeometryFEM_Implicit(dm, cell, quad, v, J, in… in DMPlexComputeCellGeometryFEM()
2691 PetscInt cdim, coordSize, d; in DMPlexComputeGeometryFVM_1D_Internal() local
2695 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexComputeGeometryFVM_1D_Internal()
2697 …k(coordSize == cdim * 2, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Edge has %" PetscInt_FMT " coordin… in DMPlexComputeGeometryFVM_1D_Internal()
2699 for (d = 0; d < cdim; ++d) centroid[d] = 0.5 * PetscRealPart(coords[d] + coords[cdim + d]); in DMPlexComputeGeometryFVM_1D_Internal()
2704 switch (cdim) { in DMPlexComputeGeometryFVM_1D_Internal()
2708 normal[0] = -PetscRealPart(coords[1] - coords[cdim + 1]); in DMPlexComputeGeometryFVM_1D_Internal()
2709 normal[1] = PetscRealPart(coords[0] - coords[cdim + 0]); in DMPlexComputeGeometryFVM_1D_Internal()
2715 … SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Dimension %" PetscInt_FMT " not supported", cdim); in DMPlexComputeGeometryFVM_1D_Internal()
2717 norm = DMPlex_NormD_Internal(cdim, normal); in DMPlexComputeGeometryFVM_1D_Internal()
2718 for (d = 0; d < cdim; ++d) normal[d] /= norm; in DMPlexComputeGeometryFVM_1D_Internal()
2722 for (d = 0; d < cdim; ++d) *vol += PetscSqr(PetscRealPart(coords[d] - coords[cdim + d])); in DMPlexComputeGeometryFVM_1D_Internal()
2738 PetscInt cdim, numCorners, p, d; in DMPlexComputeGeometryFVM_2D_Internal() local
2751 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexComputeGeometryFVM_2D_Internal()
2757 for (d = 0; d < cdim; d++) origin[d] = PetscRealPart(coords[d]); in DMPlexComputeGeometryFVM_2D_Internal()
2760 for (d = 0; d < cdim; d++) { in DMPlexComputeGeometryFVM_2D_Internal()
2761 e0[d] = PetscRealPart(coords[cdim * fv[p + 1] + d]) - origin[d]; in DMPlexComputeGeometryFVM_2D_Internal()
2762 e1[d] = PetscRealPart(coords[cdim * fv[p + 2] + d]) - origin[d]; in DMPlexComputeGeometryFVM_2D_Internal()
2772 …for (d = 0; d < cdim; d++) c[d] += a * PetscRealPart(origin[d] + coords[cdim * fv[p + 1] + d] + co… in DMPlexComputeGeometryFVM_2D_Internal()
2786 for (d = 0; d < cdim; ++d) centroid[d] = c[d]; in DMPlexComputeGeometryFVM_2D_Internal()
2788 for (d = 0; d < cdim; ++d) normal[d] = n[d]; in DMPlexComputeGeometryFVM_2D_Internal()
4140 PetscInt dim, cdim, Nc, vStart, vEnd; in DMPlexRemapGeometry() local
4143 PetscCall(DMGetCoordinateDim(dm, &cdim)); in DMPlexRemapGeometry()
4144 …(cdim <= 16, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Affine version of DMPlexRemapGeometry is currently … in DMPlexRemapGeometry()
4150 PetscInt uOff[2] = {0, cdim}; in DMPlexRemapGeometry()
4155 for (c = 0; c < cdim; ++c) coords[off + c] = f[c]; in DMPlexRemapGeometry()