1023     PetscReal          lp[3], up[3], *tmp;  in DMPlexComputeGridHash_Internal()  local
1047       up[d] = lp[d] + h[d];  in DMPlexComputeGridHash_Internal()
1058     for (PetscInt k = dlim[2 * 2 + 0]; k <= dlim[2 * 2 + 1]; ++k, lp[2] = up[2], up[2] += h[2]) {  in DMPlexComputeGridHash_Internal()
1059 …if (cdim > 2) PetscCall(DMPlexGetPlaneCellIntersection_Internal(dm, c, up, &normal[3 * 2], &upper[…  in DMPlexComputeGridHash_Internal()
1061 …on 2 (%g, %g, %g) does not intersect %s\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up…  in DMPlexComputeGridHash_Internal()
1062 … %g) intersects %" PetscInt_FMT " times\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up…  in DMPlexComputeGridHash_Internal()
1064       for (PetscInt j = dlim[1 * 2 + 0]; j <= dlim[1 * 2 + 1]; ++j, lp[1] = up[1], up[1] += h[1]) {  in DMPlexComputeGridHash_Internal()
1065 …if (cdim > 1) PetscCall(DMPlexGetPlaneCellIntersection_Internal(dm, c, up, &normal[3 * 1], &upper[…  in DMPlexComputeGridHash_Internal()
1068 …on 1 (%g, %g, %g) does not intersect %s\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up…  in DMPlexComputeGridHash_Internal()
1069 … %g) intersects %" PetscInt_FMT " times\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up…  in DMPlexComputeGridHash_Internal()
1071 …      for (PetscInt i = dlim[0 * 2 + 0]; i <= dlim[0 * 2 + 1]; ++i, lp[0] = up[0], up[0] += h[0]) {  in DMPlexComputeGridHash_Internal()
1078 …PetscCall(DMPlexGetPlaneCellIntersection_Internal(dm, c, up, &normal[3 * 0], &upper[0], &upperInt[…  in DMPlexComputeGridHash_Internal()
1081 …on 0 (%g, %g, %g) does not intersect %s\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up…  in DMPlexComputeGridHash_Internal()
1082 … %g) intersects %" PetscInt_FMT " times\n", c, (double)up[0], (double)up[1], cdim > 2 ? (double)up…  in DMPlexComputeGridHash_Internal()
1127 … + d] < (lp[d] - PETSC_SMALL)) || (lowerIntPoints[plane][ip * cdim + d] > (up[d] + PETSC_SMALL))) {  in DMPlexComputeGridHash_Internal()
1142 … + d] < (lp[d] - PETSC_SMALL)) || (upperIntPoints[plane][ip * cdim + d] > (up[d] + PETSC_SMALL))) {  in DMPlexComputeGridHash_Internal()
1143 …if (upperIntPoints[plane][ip * cdim + d] > (up[d] + PETSC_SMALL)) outsideCount[cdim + d]++; // The…  in DMPlexComputeGridHash_Internal()
1199 … - PETSC_SMALL || lowerIntPoints[plane][ip * cdim + nonIntersectingFace] < up[nonIntersectingFace]…  in DMPlexComputeGridHash_Internal()
1202 … - PETSC_SMALL || upperIntPoints[plane][ip * cdim + nonIntersectingFace] < up[nonIntersectingFace]…  in DMPlexComputeGridHash_Internal()
1220         up[0]             = lp[0] + h[0];  in DMPlexComputeGridHash_Internal()
1228       up[1]             = lp[1] + h[1];  in DMPlexComputeGridHash_Internal()