2755     PetscReal c[3] = {0., 0., 0.}, n[3] = {0., 0., 0.}, origin[3] = {0., 0., 0.}, norm;  in DMPlexComputeGeometryFVM_2D_Internal()  local
2757     for (d = 0; d < cdim; d++) origin[d] = PetscRealPart(coords[d]);  in DMPlexComputeGeometryFVM_2D_Internal()
2761         e0[d] = PetscRealPart(coords[cdim * fv[p + 1] + d]) - origin[d];  in DMPlexComputeGeometryFVM_2D_Internal()
2762         e1[d] = PetscRealPart(coords[cdim * fv[p + 2] + d]) - origin[d];  in DMPlexComputeGeometryFVM_2D_Internal()
2772 …for (d = 0; d < cdim; d++) c[d] += a * PetscRealPart(origin[d] + coords[cdim * fv[p + 1] + d] + co…  in DMPlexComputeGeometryFVM_2D_Internal()
2802   PetscReal             vsum      = 0.0, vtmp, coordsTmp[3 * 3], origin[3];  in DMPlexComputeGeometryFVM_3D_Internal()  local
2837       for (d = 0; d < dim; d++) origin[d] = PetscRealPart(coords[d]);  in DMPlexComputeGeometryFVM_3D_Internal()
2841         coordsTmp[0 * dim + d] = PetscRealPart(coords[faces[fOff + 0] * dim + d]) - origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()
2842         coordsTmp[1 * dim + d] = PetscRealPart(coords[faces[fOff + 1] * dim + d]) - origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()
2843         coordsTmp[2 * dim + d] = PetscRealPart(coords[faces[fOff + 2] * dim + d]) - origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()
2866         coordsTmp[0 * dim + d] = PetscRealPart(coords[faces[fOff + fv[0]] * dim + d]) - origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()
2867         coordsTmp[1 * dim + d] = PetscRealPart(coords[faces[fOff + fv[1]] * dim + d]) - origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()
2868         coordsTmp[2 * dim + d] = PetscRealPart(coords[faces[fOff + fv[3]] * dim + d]) - origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()
2880         coordsTmp[0 * dim + d] = PetscRealPart(coords[faces[fOff + fv[1]] * dim + d]) - origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()
2881         coordsTmp[1 * dim + d] = PetscRealPart(coords[faces[fOff + fv[2]] * dim + d]) - origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()
2882         coordsTmp[2 * dim + d] = PetscRealPart(coords[faces[fOff + fv[3]] * dim + d]) - origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()
2905     for (d = 0; d < dim; ++d) centroid[d] = centroid[d] / (vsum * 4) + origin[d];  in DMPlexComputeGeometryFVM_3D_Internal()