1023     PetscReal          lp[3], up[3], *tmp;  in DMPlexComputeGridHash_Internal()  local
1046       lp[d] = lbox->lower[d] + dlim[d * 2 + 0] * h[d];  in DMPlexComputeGridHash_Internal()
1047       up[d] = lp[d] + h[d];  in DMPlexComputeGridHash_Internal()
1050 …PetscCall(DMPlexGetPlaneCellIntersection_Internal(dm, c, lp, &normal[d * 3], &lower[d], &lowerInt[…  in DMPlexComputeGridHash_Internal()
1053 …" (%g, %g, %g) does not intersect %s\n", c, d, (double)lp[0], (double)lp[1], cdim > 2 ? (double)lp…  in DMPlexComputeGridHash_Internal()
1054 …) intersects %" PetscInt_FMT " times\n", c, d, (double)lp[0], (double)lp[1], cdim > 2 ? (double)lp…  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()
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()
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()
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()
1142 …if ((upperIntPoints[plane][ip * cdim + d] < (lp[d] - PETSC_SMALL)) || (upperIntPoints[plane][ip * …  in DMPlexComputeGridHash_Internal()
1199 …if (lowerIntPoints[plane][ip * cdim + nonIntersectingFace] > lp[nonIntersectingFace] - PETSC_SMALL…  in DMPlexComputeGridHash_Internal()
1202 …if (upperIntPoints[plane][ip * cdim + nonIntersectingFace] > lp[nonIntersectingFace] - PETSC_SMALL…  in DMPlexComputeGridHash_Internal()
1219         lp[0]             = lbox->lower[0] + dlim[0 * 2 + 0] * h[0];  in DMPlexComputeGridHash_Internal()
1220         up[0]             = lp[0] + h[0];  in DMPlexComputeGridHash_Internal()
1227       lp[1]             = lbox->lower[1] + dlim[1 * 2 + 0] * h[1];  in DMPlexComputeGridHash_Internal()
1228       up[1]             = lp[1] + h[1];  in DMPlexComputeGridHash_Internal()