90   PetscInt        face, dim, indC[3], indS[3], seenA, flippedA, seenB, flippedB, mismatch;  in DMPlexCheckFace_Internal()  local
95   face = faceFIFO[(*fTop)++];  in DMPlexCheckFace_Internal()
99   PetscCall(DMPlexGetSupportSize(dm, face, &suppSize));  in DMPlexCheckFace_Internal()
100   PetscCall(DMPlexGetSupport(dm, face, &supp));  in DMPlexCheckFace_Internal()
131     if (coneA[c] == face) posA = c;  in DMPlexCheckFace_Internal()
134 …_WRONG, "Face %" PetscInt_FMT " could not be located in cell %" PetscInt_FMT, face, supp[indS[0]]);  in DMPlexCheckFace_Internal()
144     if (coneB[c] == face) posB = c;  in DMPlexCheckFace_Internal()
147 …_WRONG, "Face %" PetscInt_FMT " could not be located in cell %" PetscInt_FMT, face, supp[indS[1]]);  in DMPlexCheckFace_Internal()
172   PetscInt        face, dim, seenA, flippedA, seenB, flippedB, mismatch, c;  in DMPlexCheckFace_Old_Internal()  local
175   face = faceFIFO[(*fTop)++];  in DMPlexCheckFace_Old_Internal()
177   PetscCall(DMPlexGetSupportSize(dm, face, &supportSize));  in DMPlexCheckFace_Old_Internal()
178   PetscCall(DMPlexGetSupport(dm, face, &support));  in DMPlexCheckFace_Old_Internal()
197     if (coneA[c] == face) posA = c;  in DMPlexCheckFace_Old_Internal()
200 …ARG_WRONG, "Face %" PetscInt_FMT " could not be located in cell %" PetscInt_FMT, face, support[0]);  in DMPlexCheckFace_Old_Internal()
206     if (coneB[c] == face) posB = c;  in DMPlexCheckFace_Old_Internal()
209 …ARG_WRONG, "Face %" PetscInt_FMT " could not be located in cell %" PetscInt_FMT, face, support[1]);  in DMPlexCheckFace_Old_Internal()
350 …PetscInt           numLeaves, numRoots, dim, h, cStart, cEnd, c, cell, fStart, fEnd, face, off, to…  in DMPlexOrient()  local
429     for (face = 0; face < fEnd - fStart; ++face) {  in DMPlexOrient()
430       if (PetscBTLookup(seenFaces, face)) faceComp[face] = comp;  in DMPlexOrient()
455     for (face = fStart; face < fEnd; ++face) {  in DMPlexOrient()
459       PetscCall(DMPlexGetSupportSize(dm, face, &supportSize));  in DMPlexOrient()
461       PetscCall(DMPlexGetSupport(dm, face, &support));  in DMPlexOrient()
477         if (cone[c] == face) break;  in DMPlexOrient()
480         if (PetscBTLookup(flippedCells, neighbor - cStart)) rorntComp[face].rank = 1 - c * 2;  in DMPlexOrient()
481         else rorntComp[face].rank = c * 2 - 1;  in DMPlexOrient()
483 …   if (PetscBTLookup(flippedCells, neighbor - cStart)) rorntComp[face].rank = ornt[c] < 0 ? -1 : 1;  in DMPlexOrient()
484         else rorntComp[face].rank = ornt[c] < 0 ? 1 : -1;  in DMPlexOrient()
486       rorntComp[face].index = faceComp[face - fStart];  in DMPlexOrient()
504       const PetscInt face = lpoints ? lpoints[l] : l;  in DMPlexOrient()  local
507 …if ((face >= fStart) && (face < fEnd) && (faceComp[face - fStart] == comp) && rorntComp[face].rank…  in DMPlexOrient()
509         const PetscInt rcomp = lorntComp[face].index;  in DMPlexOrient()
516           PetscCall(DMPlexGetSupportSize(dm, face, &supportSize));  in DMPlexOrient()
520 … (%" PetscInt_FMT ", %" PetscInt_FMT ") with orientation %" PetscInt_FMT "\n", rank, comp, l, face,  in DMPlexOrient()
521                                              rpoints[l].index, rrank, rcomp, lorntComp[face].rank));  in DMPlexOrient()
535       const PetscInt face = lpoints ? lpoints[neighbors[comp][n]] : neighbors[comp][n];  in DMPlexOrient()  local
536       const PetscInt o    = rorntComp[face].rank * lorntComp[face].rank;  in DMPlexOrient()
540 …PetscInt_FMT ") neighbor: %" PetscInt_FMT " comp: %d", face, rorntComp[face].rank, lorntComp[face]…  in DMPlexOrient()
825       const PetscInt  face = faces ? faces[f] : f;  in DMPlexOrientCells_Internal()  local
829       PetscCall(DMPlexGetSupportSize(dm, face, &suppSize));  in DMPlexOrientCells_Internal()
830       PetscCall(DMPlexGetSupport(dm, face, &supp));  in DMPlexOrientCells_Internal()
850         if (cone[c] == face) break;  in DMPlexOrientCells_Internal()
853         if (PetscBTLookup(flippedCells, nind)) rorntComp[face].rank = 1 - c * 2;  in DMPlexOrientCells_Internal()
854         else rorntComp[face].rank = c * 2 - 1;  in DMPlexOrientCells_Internal()
856         if (PetscBTLookup(flippedCells, nind)) rorntComp[face].rank = ornt[c] < 0 ? -1 : 1;  in DMPlexOrientCells_Internal()
857         else rorntComp[face].rank = ornt[c] < 0 ? 1 : -1;  in DMPlexOrientCells_Internal()
859       rorntComp[face].index = faceComp[GetPointIndex(face, fStart, fEnd, faces)];  in DMPlexOrientCells_Internal()
877       const PetscInt face = lpoints[l];  in DMPlexOrientCells_Internal()  local
881       find = GetPointIndex(face, fStart, fEnd, faces);  in DMPlexOrientCells_Internal()
882       if ((find >= 0) && (faceComp[find] == comp) && rorntComp[face].rank) {  in DMPlexOrientCells_Internal()
884         const PetscInt rcomp = lorntComp[face].index;  in DMPlexOrientCells_Internal()
892           PetscCall(DMPlexGetSupport(dm, face, &supp));  in DMPlexOrientCells_Internal()
893           PetscCall(DMPlexGetSupportSize(dm, face, &suppSize));  in DMPlexOrientCells_Internal()
898 …RR_ARG_WRONG, "Boundary face %" PetscInt_FMT " should see one cell, not %" PetscInt_FMT, face, Ns);  in DMPlexOrientCells_Internal()
900 … (%" PetscInt_FMT ", %" PetscInt_FMT ") with orientation %" PetscInt_FMT "\n", rank, comp, l, face,  in DMPlexOrientCells_Internal()
901                                              rpoints[l].index, rrank, rcomp, lorntComp[face].rank));  in DMPlexOrientCells_Internal()
913       const PetscInt face = lpoints[neighbors[comp][n]];  in DMPlexOrientCells_Internal()  local
914       const PetscInt o    = rorntComp[face].rank * lorntComp[face].rank;  in DMPlexOrientCells_Internal()
919 … ") neighbor: %" PetscInt_FMT " comp: %" PetscInt_FMT, face, rorntComp[face].rank, lorntComp[face]…  in DMPlexOrientCells_Internal()