1640   PetscInt pdim; // Dimension of approximation space P  in PetscFVCreateTabulation()  local
1655   pdim = Nc;  in PetscFVCreateTabulation()
1660   (*T)->Nb   = pdim;  in PetscFVCreateTabulation()
1664 …for (k = 0; k <= (*T)->K; ++k) PetscCall(PetscMalloc1(nrepl * npoints * pdim * Nc * PetscPowInt(cd…  in PetscFVCreateTabulation()
1667       for (d = 0; d < pdim; ++d)  in PetscFVCreateTabulation()
1668         for (c = 0; c < Nc; ++c) (*T)->T[0][(p * pdim + d) * Nc + c] = 1.;  in PetscFVCreateTabulation()
1672       for (d = 0; d < pdim; ++d)  in PetscFVCreateTabulation()
1674           for (e = 0; e < cdim; ++e) (*T)->T[1][((p * pdim + d) * Nc + c) * cdim + e] = 0.0;  in PetscFVCreateTabulation()
1678       for (d = 0; d < pdim; ++d)  in PetscFVCreateTabulation()
1680 …   for (e = 0; e < cdim * cdim; ++e) (*T)->T[2][((p * pdim + d) * Nc + c) * cdim * cdim + e] = 0.0;  in PetscFVCreateTabulation()
1913   PetscInt           dim, numConstants, pdim, totDim, Nc, off, f, d;  in PetscFVIntegrateRHSFunction_Upwind()  local
1923   PetscCall(PetscFVGetNumComponents(fvm, &pdim));  in PetscFVIntegrateRHSFunction_Upwind()
1925 …(*riemann)(dim, pdim, fgeom[f].centroid, fgeom[f].normal, &uL[f * Nc], &uR[f * Nc], numConstants, …  in PetscFVIntegrateRHSFunction_Upwind()
1926     for (d = 0; d < pdim; ++d) {  in PetscFVIntegrateRHSFunction_Upwind()
2195   PetscInt           dim, numConstants, pdim, Nc, totDim, off, f, d;  in PetscFVIntegrateRHSFunction_LeastSquares()  local
2205   PetscCall(PetscFVGetNumComponents(fvm, &pdim));  in PetscFVIntegrateRHSFunction_LeastSquares()
2207 …(*riemann)(dim, pdim, fgeom[f].centroid, fgeom[f].normal, &uL[f * Nc], &uR[f * Nc], numConstants, …  in PetscFVIntegrateRHSFunction_LeastSquares()
2208     for (d = 0; d < pdim; ++d) {  in PetscFVIntegrateRHSFunction_LeastSquares()