4098   PetscInt           dim, Nf, f, Nc, numFaces = fEnd - fStart, iface, face;  in DMPlexGetFaceFields()  local
4142   for (face = fStart, iface = 0; face < fEnd; ++face) {  in DMPlexGetFaceFields()
4182 …      PetscCall(PetscFEEvaluateFaceFields_Internal(prob, f, faceLocL, xL, &uLl[iface * Nc + off]));  in DMPlexGetFaceFields()
4183 …izeR) PetscCall(PetscFEEvaluateFaceFields_Internal(prob, f, faceLocR, xR, &uRl[iface * Nc + off]));  in DMPlexGetFaceFields()
4185             for (d = 0; d < comp; ++d) uRl[iface * Nc + off + d] = uLl[iface * Nc + off + d];  in DMPlexGetFaceFields()
4188 …      PetscCall(PetscFEEvaluateFaceFields_Internal(prob, f, faceLocR, xR, &uRl[iface * Nc + off]));  in DMPlexGetFaceFields()
4190           for (d = 0; d < comp; ++d) uLl[iface * Nc + off + d] = uRl[iface * Nc + off + d];  in DMPlexGetFaceFields()
4210             uLl[iface * Nc + off + c] = xL[c] + DMPlex_DotD_Internal(dim, &gL[c * dim], dxL);  in DMPlexGetFaceFields()
4211             uRl[iface * Nc + off + c] = xR[c] + DMPlex_DotD_Internal(dim, &gR[c * dim], dxR);  in DMPlexGetFaceFields()
4215             uLl[iface * Nc + off + c] = xL[c];  in DMPlexGetFaceFields()
4216             uRl[iface * Nc + off + c] = xR[c];  in DMPlexGetFaceFields()
4221     ++iface;  in DMPlexGetFaceFields()
4223   *Nface = iface;  in DMPlexGetFaceFields()
4286   PetscInt           dim, numFaces = fEnd - fStart, iface, face;  in DMPlexGetFaceGeometry()  local
4302   for (face = fStart, iface = 0; face < fEnd; ++face) {  in DMPlexGetFaceGeometry()
4319       fgeoml[iface].centroid[d] = fg->centroid[d];  in DMPlexGetFaceGeometry()
4320       fgeoml[iface].normal[d]   = fg->normal[d];  in DMPlexGetFaceGeometry()
4322     voll[iface * 2 + 0] = cgL->volume;  in DMPlexGetFaceGeometry()
4323     voll[iface * 2 + 1] = cgR->volume;  in DMPlexGetFaceGeometry()
4324     ++iface;  in DMPlexGetFaceGeometry()
4326   *Nface = iface;  in DMPlexGetFaceGeometry()
5366       PetscInt     iface;  in DMPlexComputeResidualByKey()  local
5383         for (face = fS, iface = 0; face < fE; ++face) {  in DMPlexComputeResidualByKey()
5398             PetscCall(DMPrintCellVectorReal(face, "Residual: normal", cdim, fgeom[iface].normal));  in DMPlexComputeResidualByKey()
5399 …      PetscCall(DMPrintCellVector(face, "Residual: left state", pdim, &uL[iface * totDim + foff]));  in DMPlexComputeResidualByKey()
5400 …     PetscCall(DMPrintCellVector(face, "Residual: right state", pdim, &uR[iface * totDim + foff]));  in DMPlexComputeResidualByKey()
5401 …    PetscCall(DMPrintCellVector(face, "Residual: left flux", pdim, &fluxL[iface * totDim + foff]));  in DMPlexComputeResidualByKey()
5402 …   PetscCall(DMPrintCellVector(face, "Residual: right flux", pdim, &fluxR[iface * totDim + foff]));  in DMPlexComputeResidualByKey()
5405             if (fL) fL[d] -= fluxL[iface * totDim + foff + d];  in DMPlexComputeResidualByKey()
5406             if (fR) fR[d] += fluxR[iface * totDim + foff + d];  in DMPlexComputeResidualByKey()
5408           ++iface;  in DMPlexComputeResidualByKey()