3020     PetscFVFaceGeom *fg;  in DMPlexComputeGeometryFVM()  local
3031     PetscCall(DMPlexPointLocalRef(dmFace, f, fgeom, &fg));  in DMPlexComputeGeometryFVM()
3032     PetscCall(DMPlexComputeCellGeometryFVM(dm, f, &area, fg->centroid, fg->normal));  in DMPlexComputeGeometryFVM()
3033     for (d = 0; d < dim; ++d) fg->normal[d] *= area;  in DMPlexComputeGeometryFVM()
3041       lcentroid = cells[0] >= cEndInterior ? fg->centroid : cL->centroid;  in DMPlexComputeGeometryFVM()
3044         rcentroid = cells[1] >= cEndInterior ? fg->centroid : cR->centroid;  in DMPlexComputeGeometryFVM()
3046         rcentroid = fg->centroid;  in DMPlexComputeGeometryFVM()
3048       PetscCall(DMLocalizeCoordinateReal_Internal(dm, dim, fg->centroid, lcentroid, l));  in DMPlexComputeGeometryFVM()
3049       PetscCall(DMLocalizeCoordinateReal_Internal(dm, dim, fg->centroid, rcentroid, r));  in DMPlexComputeGeometryFVM()
3051       if (DMPlex_DotRealD_Internal(dim, fg->normal, v) < 0) {  in DMPlexComputeGeometryFVM()
3052         for (d = 0; d < dim; ++d) fg->normal[d] = -fg->normal[d];  in DMPlexComputeGeometryFVM()
3054       if (DMPlex_DotRealD_Internal(dim, fg->normal, v) <= 0) {  in DMPlexComputeGeometryFVM()
3055 …FMT " could not be fixed, normal (%g,%g) v (%g,%g)", f, (double)fg->normal[0], (double)fg->normal[…  in DMPlexComputeGeometryFVM()
3056 … fixed, normal (%g,%g,%g) v (%g,%g,%g)", f, (double)fg->normal[0], (double)fg->normal[1], (double)…  in DMPlexComputeGeometryFVM()
3060         DMPlex_WaxpyD_Internal(dim, -1, fg->centroid, cL->centroid, v);  in DMPlexComputeGeometryFVM()
3064         DMPlex_WaxpyD_Internal(dim, -1, fg->centroid, cR->centroid, v);  in DMPlexComputeGeometryFVM()
3073     PetscFVFaceGeom *fg;  in DMPlexComputeGeometryFVM()  local
3083     PetscCall(DMPlexPointLocalRef(dmFace, cone[0], fgeom, &fg));  in DMPlexComputeGeometryFVM()
3092 …      DMPlex_WaxpyD_Internal(dim, -1, ci->centroid, fg->centroid, c2f); /* cell to face centroid */  in DMPlexComputeGeometryFVM()
3093 …a = DMPlex_DotRealD_Internal(dim, c2f, fg->normal) / DMPlex_DotRealD_Internal(dim, fg->normal, fg-…  in DMPlexComputeGeometryFVM()
3095         DMPlex_WaxpyD_Internal(dim, 2 * a, fg->normal, ci->centroid, cg->centroid);  in DMPlexComputeGeometryFVM()
3235       PetscFVFaceGeom *fg;  in BuildGradientReconstruction_Internal()  local
3245       PetscCall(DMPlexPointLocalRef(dmFace, faces[f], fgeom, &fg));  in BuildGradientReconstruction_Internal()
3248       gref[usedFaces++] = fg->grad[side]; /* Gradient reconstruction term will go here */  in BuildGradientReconstruction_Internal()
3351       PetscFVFaceGeom *fg;  in BuildGradientReconstruction_Internal_Tree()  local
3359       PetscCall(DMPlexPointLocalRef(dmFace, nface, fgeom, &fg));  in BuildGradientReconstruction_Internal_Tree()
3362       gref[f] = fg->grad[side]; /* Gradient reconstruction term will go here */  in BuildGradientReconstruction_Internal_Tree()