2021         PetscReal extPoint[4];  in DMPlexComputeRectangleGeometry_Internal()  local
2023         extPoint[0] = 1.;  in DMPlexComputeRectangleGeometry_Internal()
2024         extPoint[1] = xi;  in DMPlexComputeRectangleGeometry_Internal()
2025         extPoint[2] = eta;  in DMPlexComputeRectangleGeometry_Internal()
2026         extPoint[3] = xi * eta;  in DMPlexComputeRectangleGeometry_Internal()
2030           for (k = 0; k < Nv; k++) val += extPoint[k] * zCoeff[dim * k + j];  in DMPlexComputeRectangleGeometry_Internal()
2162         PetscReal extPoint[8];  in DMPlexComputeHexahedronGeometry_Internal()  local
2164         extPoint[0] = 1.;  in DMPlexComputeHexahedronGeometry_Internal()
2165         extPoint[1] = xi;  in DMPlexComputeHexahedronGeometry_Internal()
2166         extPoint[2] = eta;  in DMPlexComputeHexahedronGeometry_Internal()
2167         extPoint[3] = xi * eta;  in DMPlexComputeHexahedronGeometry_Internal()
2168         extPoint[4] = theta;  in DMPlexComputeHexahedronGeometry_Internal()
2169         extPoint[5] = theta * xi;  in DMPlexComputeHexahedronGeometry_Internal()
2170         extPoint[6] = theta * eta;  in DMPlexComputeHexahedronGeometry_Internal()
2171         extPoint[7] = theta * eta * xi;  in DMPlexComputeHexahedronGeometry_Internal()
2175           for (k = 0; k < Nv; k++) val += extPoint[k] * zCoeff[dim * k + j];  in DMPlexComputeHexahedronGeometry_Internal()
2310         PetscReal extPoint[6];  in DMPlexComputeTriangularPrismGeometry_Internal()  local
2313         extPoint[0] = 1.;  in DMPlexComputeTriangularPrismGeometry_Internal()
2314         extPoint[1] = eta;  in DMPlexComputeTriangularPrismGeometry_Internal()
2315         extPoint[2] = xi;  in DMPlexComputeTriangularPrismGeometry_Internal()
2316         extPoint[3] = zeta;  in DMPlexComputeTriangularPrismGeometry_Internal()
2317         extPoint[4] = xi * zeta;  in DMPlexComputeTriangularPrismGeometry_Internal()
2318         extPoint[5] = eta * zeta;  in DMPlexComputeTriangularPrismGeometry_Internal()
2322           for (k = 0; k < Nv; ++k) val += extPoint[k] * coeff[k * dim + c];  in DMPlexComputeTriangularPrismGeometry_Internal()