154     PetscInt     Nf1, Nbpt1, Nk1;  in test()  local
166 …s(dim, deg, form < 0 ? form - 1 : form + 1, 0, npoints, points, weights, &Nbpt1, &Nf1, &Nk1, &p_tr…  in test()
172     PetscCall(PetscCalloc1(Nbpt1 * Nf * npoints, &p_koszul));  in test()
173     for (PetscInt b = 0; b < Nbpt1; b++) {  in test()
195     PetscCall(PetscMalloc1(Nbpt1 * Nbpt1, &M_koszul));  in test()
196     for (PetscInt i = 0; i < Nbpt1; i++) {  in test()
197       for (PetscInt j = 0; j < Nbpt1; j++) {  in test()
206         M_koszul[i * Nbpt1 + j] = val;  in test()
211     PetscCall(PetscMalloc1(Nbpt * Nbpt1, &M_k_moment));  in test()
212     for (PetscInt i = 0; i < Nbpt1; i++) {  in test()
227     PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, Nbpt1, Nbpt1, M_koszul, &mat_koszul));  in test()
228     PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, Nbpt, Nbpt1, M_k_moment, &mat_k_moment_T));  in test()