2256     PetscReal      coeff[18];  in DMPlexComputeTriangularPrismGeometry_Internal()  local
2291 …coeff[dim * 0 + j] = (1. / 4.) * (verts[dim * 1 + j] + verts[dim * 2 + j] + verts[dim * 4 + j] + v…  in DMPlexComputeTriangularPrismGeometry_Internal()
2292 …coeff[dim * 1 + j] = (1. / 4.) * (-verts[dim * 0 + j] + verts[dim * 1 + j] - verts[dim * 3 + j] + …  in DMPlexComputeTriangularPrismGeometry_Internal()
2293 …coeff[dim * 2 + j] = (1. / 4.) * (-verts[dim * 0 + j] + verts[dim * 2 + j] - verts[dim * 3 + j] + …  in DMPlexComputeTriangularPrismGeometry_Internal()
2294 …coeff[dim * 3 + j] = (1. / 4.) * (-verts[dim * 1 + j] - verts[dim * 2 + j] + verts[dim * 4 + j] + …  in DMPlexComputeTriangularPrismGeometry_Internal()
2295 …coeff[dim * 4 + j] = (1. / 4.) * (verts[dim * 0 + j] - verts[dim * 2 + j] - verts[dim * 3 + j] + v…  in DMPlexComputeTriangularPrismGeometry_Internal()
2296 …coeff[dim * 5 + j] = (1. / 4.) * (verts[dim * 0 + j] - verts[dim * 1 + j] - verts[dim * 3 + j] + v…  in DMPlexComputeTriangularPrismGeometry_Internal()
2322           for (k = 0; k < Nv; ++k) val += extPoint[k] * coeff[k * dim + c];  in DMPlexComputeTriangularPrismGeometry_Internal()
2352             for (l = 0; l < Nv; l++) val += coeff[dim * l + j] * extJ[dimR * l + k];  in DMPlexComputeTriangularPrismGeometry_Internal()
3631           PetscReal coeff = cellCoeffs[dimC * k + l];  in DMPlexCoordinatesToReference_Tensor()  local
3633           resNeg[l] -= coeff * extCoord;  in DMPlexCoordinatesToReference_Tensor()
3634           for (m = 0; m < dimR; m++) J[dimR * l + m] += coeff * extJ[m];  in DMPlexCoordinatesToReference_Tensor()
3716         PetscReal coeff = cellCoeffs[dimC * k + l];  in DMPlexReferenceToCoordinates_Tensor()  local
3718         mapped[l] += coeff * extCoord;  in DMPlexReferenceToCoordinates_Tensor()