6 static PetscErrorCode testOrthogonality(PetscInt dim, PetscInt deg)  in testOrthogonality()  argument
14   PetscCall(PetscDTStroudConicalQuadrature(dim, 1, deg + 1, -1., 1., &q));  in testOrthogonality()
16   PetscCall(PetscDTBinomialInt(dim + deg, dim, &Npoly));  in testOrthogonality()
18   PetscCall(PetscDTPKDEvalJet(dim, npoints, points, deg, 0, p));  in testOrthogonality()
33 static PetscErrorCode testDerivativesLegendre(PetscInt dim, PetscInt deg, PetscInt k)  in testDerivativesLegendre()  argument
53   PetscCall(PetscDTBinomialInt(dim + deg, dim, &Np));  in testDerivativesLegendre()
57   PetscCall(PetscDTStroudConicalQuadrature(dim, 1, deg + 1, -1., 1., &q));  in testDerivativesLegendre()
60   PetscCall(PetscDTPKDEvalJet(dim, npoints, points, deg, 0, proj));  in testDerivativesLegendre()
72   PetscCall(PetscMalloc1(deg + 1, &degrees));  in testDerivativesLegendre()
73   for (i = 0; i < deg + 1; i++) degrees[i] = i;  in testDerivativesLegendre()
80     PetscCall(PetscMalloc1((deg + 1) * npoints, &B[d]));  in testDerivativesLegendre()
83     PetscCall(PetscDTLegendreEval(npoints, points1d, deg + 1, degrees, B[d], NULL, NULL));  in testDerivativesLegendre()
94       for (l = 0; l < dim; l++) valj *= B[l][i * (deg + 1) + degtup[l]];  in testDerivativesLegendre()
114   PetscCall(PetscDTPKDEvalJet(dim, 1, point, deg, k, pkd_jet_basis));  in testDerivativesLegendre()
124     PetscCall(PetscMalloc1((deg + 1) * (k + 1), &D[i]));  in testDerivativesLegendre()
125     PetscCall(PetscDTJacobiEvalJet(0., 0., 1, &point[i], deg, k, D[i]));  in testDerivativesLegendre()
165   PetscInt dim, deg, k;  in main()  local
170   deg = 4;  in main()
174 …PetscCall(PetscOptionsInt("-degree", "The degree of the polynomial space", "ex9.c", deg, &deg, NUL…  in main()
177   PetscCall(testOrthogonality(dim, deg));  in main()
178   PetscCall(testDerivativesLegendre(dim, deg, k));  in main()