9   PetscInt   Nf;   // Number of form components  in constructTabulationAndMass()  local
15   PetscCall(PetscDTBinomialInt(dim, PetscAbsInt(form), &Nf));  in constructTabulationAndMass()
18   PetscCall(PetscMalloc1(Nbpt * Nf * Nk * npoints, &p_trimmed));  in constructTabulationAndMass()
28       for (PetscInt f = 0; f < Nf; f++) {  in constructTabulationAndMass()
29         const PetscReal *p_i = &p_trimmed[(i * Nf + f) * Nk * npoints];  in constructTabulationAndMass()
30         const PetscReal *p_j = &p_trimmed[(j * Nf + f) * Nk * npoints];  in constructTabulationAndMass()
38   *_Nf = Nf;  in constructTabulationAndMass()
51   PetscInt         Nf;   // Number of form components  in test()  local
74 …lationAndMass(dim, deg, form, jetDegree, npoints, points, weights, &Nbpt, &Nf, &Nk, &p_trimmed, &M…  in test()
103   PetscCall(PetscCalloc1(Nbpt * Nbp * Nf, &M_moments));  in test()
106       for (PetscInt f = 0; f < Nf; f++) {  in test()
109         const PetscReal *p_j = &p_trimmed[(j * Nf + f) * Nk * npoints];  in test()
112         M_moments[(i * Nf + f) * Nbpt + j] += v;  in test()
127     PetscCall(PetscBLASIntCast(Nbp * Nf, &k));  in test()
140   PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, Nbpt, Nbm1 * Nf, M_moments, &mat_moments_T));  in test()
172     PetscCall(PetscCalloc1(Nbpt1 * Nf * npoints, &p_koszul));  in test()
181         if (form < 0) i = Nf - 1 - i;  in test()
188         p_i = &p_koszul[(b * Nf + i) * npoints];  in test()
200         for (PetscInt v = 0; v < Nf; v++) {  in test()
201           const PetscReal *p_i = &p_koszul[(i * Nf + v) * npoints];  in test()
202           const PetscReal *p_j = &p_koszul[(j * Nf + v) * npoints];  in test()
216         for (PetscInt v = 0; v < Nf; v++) {  in test()
217           const PetscReal *p_i = &p_koszul[(i * Nf + v) * npoints];  in test()
218           const PetscReal *p_j = &p_trimmed[(j * Nf + v) * Nk * npoints];  in test()
252   PetscCall(PetscCalloc1(Nbpt * Nf * Nk * npoints, &p_trimmed_copy));  in test()
253   PetscCall(PetscMalloc1(Nbp * Nf * Nbpt, &M_moment_real));  in test()
254   for (PetscInt i = 0; i < Nbp * Nf * Nbpt; i++) M_moment_real[i] = PetscRealPart(M_moments[i]);  in test()
255   for (PetscInt f = 0; f < Nf; f++) {  in test()
264     PetscCall(PetscBLASIntCast(Nf * Nbpt, &ldb));  in test()
265     PetscCall(PetscBLASIntCast(Nf * Nk * npoints, &ldc));  in test()
269 …for (PetscInt i = 0; i < Nbpt * Nf * Nk * npoints; i++) frob_err += (p_trimmed_copy[i] - p_trimmed…  in test()