10   PetscInt   Nbpt; // number of trimmed polynomials  in constructTabulationAndMass()  local
16   PetscCall(PetscDTPTrimmedSize(dim, deg, form, &Nbpt));  in constructTabulationAndMass()
18   PetscCall(PetscMalloc1(Nbpt * Nf * Nk * npoints, &p_trimmed));  in constructTabulationAndMass()
23   PetscCall(PetscCalloc1(Nbpt * Nbpt, &M_trimmed));  in constructTabulationAndMass()
24   for (PetscInt i = 0; i < Nbpt; i++) {  in constructTabulationAndMass()
25     for (PetscInt j = 0; j < Nbpt; j++) {  in constructTabulationAndMass()
34       M_trimmed[i * Nbpt + j] += v;  in constructTabulationAndMass()
37   *_Nb = Nbpt;  in constructTabulationAndMass()
53   PetscInt         Nbpt; // number of trimmed polynomials  in test()  local
74 …PetscCall(constructTabulationAndMass(dim, deg, form, jetDegree, npoints, points, weights, &Nbpt, &…  in test()
80   PetscCall(PetscMalloc1(Nbpt * Nbpt, &Mcopy));  in test()
88     PetscCall(PetscBLASIntCast(Nbpt, &n));  in test()
89     PetscCall(PetscBLASIntCast(5 * Nbpt, &lwork));  in test()
90     PetscCall(PetscMalloc1(Nbpt, &S));  in test()
91     PetscCall(PetscMalloc1(5 * Nbpt, &work));  in test()
92     PetscCall(PetscArraycpy(Mcopy, M_trimmed, Nbpt * Nbpt));  in test()
95     PetscReal cond = S[0] / S[Nbpt - 1];  in test()
103   PetscCall(PetscCalloc1(Nbpt * Nbp * Nf, &M_moments));  in test()
105     for (PetscInt j = 0; j < Nbpt; j++) {  in test()
112         M_moments[(i * Nf + f) * Nbpt + j] += v;  in test()
119   PetscCall(PetscArraycpy(Mcopy, M_trimmed, Nbpt * Nbpt));  in test()
125     PetscCall(PetscBLASIntCast(Nbpt, &m));  in test()
126     PetscCall(PetscBLASIntCast(Nbpt, &n));  in test()
132 …for (PetscInt i = 0; i < Nbpt * Nbpt; i++) frob_err += PetscRealPart(Mcopy[i]) * PetscRealPart(Mco…  in test()
138   PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, Nbpt, Nbpt, M_trimmed, &mat_trimmed));  in test()
140   PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, Nbpt, Nbm1 * Nf, M_moments, &mat_moments_T));  in test()
211     PetscCall(PetscMalloc1(Nbpt * Nbpt1, &M_k_moment));  in test()
213       for (PetscInt j = 0; j < Nbpt; j++) {  in test()
222         M_k_moment[i * Nbpt + j] = val;  in test()
228     PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, Nbpt, Nbpt1, M_k_moment, &mat_k_moment_T));  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()
261     PetscCall(PetscBLASIntCast(Nbpt, &n));  in test()
264     PetscCall(PetscBLASIntCast(Nf * Nbpt, &ldb));  in test()
266 …EALgemm_("N", "T", &m, &n, &k, &alpha, p_scalar, &lda, &M_moment_real[f * Nbpt], &ldb, &beta, &p_t…  in test()
269 …for (PetscInt i = 0; i < Nbpt * Nf * Nk * npoints; i++) frob_err += (p_trimmed_copy[i] - p_trimmed…  in test()