11   PetscInt   Nk;   // jet size  in constructTabulationAndMass()  local
17   PetscCall(PetscDTBinomialInt(dim + jetDegree, dim, &Nk));  in constructTabulationAndMass()
18   PetscCall(PetscMalloc1(Nbpt * Nf * Nk * npoints, &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()
39   *_Nk = Nk;  in constructTabulationAndMass()
52   PetscInt         Nk;   // jet size  in test()  local
74 …nAndMass(dim, deg, form, jetDegree, npoints, points, weights, &Nbpt, &Nf, &Nk, &p_trimmed, &M_trim…  in test()
77   PetscCall(PetscMalloc1(Nbp * Nk * npoints, &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()
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()
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()