1259 …PetscInt        *pMaxNew;    /* The first replicated point at each depth in the new mesh, hybrids …  in DMPlexConstructCohesiveCells_Internal()  local
1272 …epthMax, depth + 1, &depthEnd, 2 * (depth + 1), &depthShift, depth + 1, &pMaxNew, depth + 1, &numH…  in DMPlexConstructCohesiveCells_Internal()
1275     PetscCall(DMPlexGetDepthStratum(dm, d, NULL, &pMaxNew[d]));  in DMPlexConstructCohesiveCells_Internal()
1277     depthEnd[d]           = pMaxNew[d];  in DMPlexConstructCohesiveCells_Internal()
1326 …for (d = 0; d <= depth; ++d) pMaxNew[d] = DMPlexShiftPoint_Internal(pMaxNew[d], depth, depthShift)…  in DMPlexConstructCohesiveCells_Internal()
1333       const PetscInt  splitp = p + pMaxNew[dep];  in DMPlexConstructCohesiveCells_Internal()
1345         const PetscInt hybcell = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1372         const PetscInt hybedge = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1392         const PetscInt hybface = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1425 …     const PetscInt hybedge = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1] + numSplitPoints[dep];  in DMPlexConstructCohesiveCells_Internal()
1446 …     const PetscInt hybface = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1] + numSplitPoints[dep];  in DMPlexConstructCohesiveCells_Internal()
1477       const PetscInt  splitp = p + pMaxNew[dep];  in DMPlexConstructCohesiveCells_Internal()
1490         const PetscInt  hybcell    = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1507             coneNew[2 + q] = v + pMaxNew[dep - 1];  in DMPlexConstructCohesiveCells_Internal()
1595             coneNew[2 + q] = v + pMaxNew[dep] + numSplitPoints[dep] + numSplitPoints[dep - 1];  in DMPlexConstructCohesiveCells_Internal()
1597             coneNew[2 + q] = v + pMaxNew[dep] + numSplitPoints[dep];  in DMPlexConstructCohesiveCells_Internal()
1606         const PetscInt hybedge = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1625             supportNew[qp++] = edge + pMaxNew[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1643             supportNew[qf++] = edge + pMaxNew[dep + 2] + numSplitPoints[dep + 2];  in DMPlexConstructCohesiveCells_Internal()
1648         const PetscInt hybface = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1659             coneNew[q] = v + pMaxNew[dep - 1];  in DMPlexConstructCohesiveCells_Internal()
1674         supportNew[q++] = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1684             supportNew[q++] = face + pMaxNew[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1689         supportNew[q++] = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1700             coneNew[2 + v] = vertex + pMaxNew[dep] + numSplitPoints[dep] + numSplitPoints[dep - 1];  in DMPlexConstructCohesiveCells_Internal()
1702             coneNew[2 + v] = vertex + pMaxNew[dep] + numSplitPoints[dep];  in DMPlexConstructCohesiveCells_Internal()
1713             supportNew[qf++] = face + pMaxNew[dep + 2] + numSplitPoints[dep + 2];  in DMPlexConstructCohesiveCells_Internal()
1732 …     const PetscInt hybedge = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1] + numSplitPoints[dep];  in DMPlexConstructCohesiveCells_Internal()
1739           if (e >= 0) supportNew[q++] = e + pMaxNew[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1756             supportNew[qf++] = edge + pMaxNew[dep + 2] + numSplitPoints[dep + 2];  in DMPlexConstructCohesiveCells_Internal()
1760 …    supportNew[qf++] = edge + pMaxNew[dep + 2] + numSplitPoints[dep + 2] + numSplitPoints[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1765 …     const PetscInt hybface = p + pMaxNew[dep + 1] + numSplitPoints[dep + 1] + numSplitPoints[dep];  in DMPlexConstructCohesiveCells_Internal()
1776             supportNew[qf++] = face + pMaxNew[dep + 1];  in DMPlexConstructCohesiveCells_Internal()
1791         coneNew[2] = v + pMaxNew[dep] + numSplitPoints[dep] + numSplitPoints[dep - 1];  in DMPlexConstructCohesiveCells_Internal()
1794         coneNew[3] = v + pMaxNew[dep] + numSplitPoints[dep] + numSplitPoints[dep - 1];  in DMPlexConstructCohesiveCells_Internal()
1802             supportNew[qf++] = face + pMaxNew[dep + 2] + numSplitPoints[dep + 2];  in DMPlexConstructCohesiveCells_Internal()
1843           csplitp = pMaxNew[dep - 1] + cp;  in DMPlexConstructCohesiveCells_Internal()
1862     const PetscInt splitp = pMaxNew[0] + v;  in DMPlexConstructCohesiveCells_Internal()
1881       const PetscInt splitp = pMaxNew[dep] + p;  in DMPlexConstructCohesiveCells_Internal()
1922     pMaxNew[d] = pEnd - numHybridPoints[d] - numHybridPointsOld[d];  in DMPlexConstructCohesiveCells_Internal()
1925   PetscCall(PetscFree5(depthMax, depthEnd, depthShift, pMaxNew, numHybridPointsOld));  in DMPlexConstructCohesiveCells_Internal()