7 …(PetscInt dim, PetscInt deg, PetscInt form, PetscInt jetDegree, PetscInt npoints, const PetscReal …  in constructTabulationAndMass()  argument
18   PetscCall(PetscMalloc1(Nbpt * Nf * Nk * npoints, &p_trimmed));  in constructTabulationAndMass()
19   PetscCall(PetscDTPTrimmedEvalJet(dim, npoints, points, deg, form, jetDegree, p_trimmed));  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()
32         for (PetscInt pt = 0; pt < npoints; pt++) v += p_i[pt] * p_j[pt] * weights[pt];  in constructTabulationAndMass()
48   PetscInt         npoints;  in test()  local
72   PetscCall(PetscQuadratureGetData(q, NULL, NULL, &npoints, &points, &weights));  in test()
74 …PetscCall(constructTabulationAndMass(dim, deg, form, jetDegree, npoints, points, weights, &Nbpt, &…  in test()
77   PetscCall(PetscMalloc1(Nbp * Nk * npoints, &p_scalar));  in test()
78   PetscCall(PetscDTPKDEvalJet(dim, npoints, points, deg, jetDegree, p_scalar));  in test()
108         const PetscReal *p_i = &p_scalar[i * Nk * npoints];  in test()
109         const PetscReal *p_j = &p_trimmed[(j * Nf + f) * Nk * npoints];  in test()
111         for (PetscInt pt = 0; pt < npoints; pt++) v += p_i[pt] * p_j[pt] * weights[pt];  in test()
166 …PetscCall(constructTabulationAndMass(dim, deg, form < 0 ? form - 1 : form + 1, 0, npoints, points,…  in test()
172     PetscCall(PetscCalloc1(Nbpt1 * Nf * npoints, &p_koszul));  in test()
188         p_i = &p_koszul[(b * Nf + i) * npoints];  in test()
189         p_j = &p_trimmed1[(b * Nf1 + j) * npoints];  in test()
190         for (PetscInt pt = 0; pt < npoints; pt++) p_i[pt] += p_j[pt] * points[pt * dim + k] * sign;  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()
204           for (PetscInt pt = 0; pt < npoints; pt++) val += p_i[pt] * p_j[pt] * weights[pt];  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()
220           for (PetscInt pt = 0; pt < npoints; pt++) val += p_i[pt] * p_j[pt] * weights[pt];  in test()
252   PetscCall(PetscCalloc1(Nbpt * Nf * Nk * npoints, &p_trimmed_copy));  in test()
260     PetscCall(PetscBLASIntCast(Nk * npoints, &m));  in test()
263     PetscCall(PetscBLASIntCast(Nk * npoints, &lda));  in test()
265     PetscCall(PetscBLASIntCast(Nf * Nk * npoints, &ldc));  in test()
266 …, p_scalar, &lda, &M_moment_real[f * Nbpt], &ldb, &beta, &p_trimmed_copy[f * 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()