7 static PetscErrorCode constructTabulationAndMass(PetscInt dim, PetscInt deg, PetscInt form, PetscIn…  in constructTabulationAndMass()  argument
15   PetscCall(PetscDTBinomialInt(dim, PetscAbsInt(form), &Nf));  in constructTabulationAndMass()
16   PetscCall(PetscDTPTrimmedSize(dim, deg, form, &Nbpt));  in constructTabulationAndMass()
19   PetscCall(PetscDTPTrimmedEvalJet(dim, npoints, points, deg, form, jetDegree, p_trimmed));  in constructTabulationAndMass()
45 static PetscErrorCode test(PetscInt dim, PetscInt deg, PetscInt form, PetscInt jetDegree, PetscBool…  in test()  argument
74 …PetscCall(constructTabulationAndMass(dim, deg, form, jetDegree, npoints, points, weights, &Nbpt, &…  in test()
96 … %" PetscInt_FMT ", form %" PetscInt_FMT ": condition number %g\n", dim, deg, form, (double)cond));  in test()
135 …Int_FMT ", form %" PetscInt_FMT ": trimmed projection error %g", dim, deg, form, (double)frob_err);  in test()
147 …", form %" PetscInt_FMT ": trimmed reverse projection error %g", dim, deg, form, (double)frob_err);  in test()
153   if (PetscAbsInt(form) < dim) {  in test()
166 …PetscCall(constructTabulationAndMass(dim, deg, form < 0 ? form - 1 : form + 1, 0, npoints, points,…  in test()
168     PetscCall(PetscMalloc1(Nf1 * (PetscAbsInt(form) + 1), &pattern));  in test()
169     PetscCall(PetscDTAltVInteriorPattern(dim, PetscAbsInt(form) + 1, pattern));  in test()
174       for (PetscInt a = 0; a < Nf1 * (PetscAbsInt(form) + 1); a++) {  in test()
181         if (form < 0) i = Nf - 1 - i;  in test()
183         if (form < 0) j = Nf1 - 1 - j;  in test()
186         if (form < 0 && (i & 1) ^ (j & 1)) sign = -sign;  in test()
234 …%" PetscInt_FMT "): koszul projection error %g", dim, deg, form, form < 0 ? (form - 1) : (form + 1…  in test()
272 …degree %" PetscInt_FMT ", form %" PetscInt_FMT ": jet error %g", dim, deg, form, (double)frob_err);  in test()
303 …for (PetscInt form = -dim + 1; form <= dim; form++) PetscCall(test(dim, deg, form, PetscMax(1, k),…  in main()  local