2802   PetscReal             vsum      = 0.0, vtmp, coordsTmp[3 * 3], origin[3];  in DMPlexComputeGeometryFVM_3D_Internal()  local
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()
2845       Volume_Tetrahedron_Origin_Internal(&vtmp, coordsTmp);  in DMPlexComputeGeometryFVM_3D_Internal()
2850           for (p = 0; p < 3; ++p) centroid[d] += coordsTmp[p * dim + d] * vtmp;  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()
2870       Volume_Tetrahedron_Origin_Internal(&vtmp, coordsTmp);  in DMPlexComputeGeometryFVM_3D_Internal()
2875           for (p = 0; p < 3; ++p) centroid[d] += coordsTmp[p * dim + d] * vtmp;  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()
2884       Volume_Tetrahedron_Origin_Internal(&vtmp, coordsTmp);  in DMPlexComputeGeometryFVM_3D_Internal()
2889           for (p = 0; p < 3; ++p) centroid[d] += coordsTmp[p * dim + d] * vtmp;  in DMPlexComputeGeometryFVM_3D_Internal()