Lines Matching refs:npoints
10 PetscInt Npoly, npoints, i, j, k; in testOrthogonality() local
15 PetscCall(PetscQuadratureGetData(q, NULL, NULL, &npoints, &points, &weights)); in testOrthogonality()
17 PetscCall(PetscMalloc1(Npoly * npoints, &p)); in testOrthogonality()
18 PetscCall(PetscDTPKDEvalJet(dim, npoints, points, deg, 0, p)); in testOrthogonality()
24 … for (k = 0; k < npoints; k++) integral += weights[k] * p[i * npoints + k] * p[j * npoints + k]; in testOrthogonality()
35 PetscInt Np, Nk, i, j, l, d, npoints; in testDerivativesLegendre() local
58 PetscCall(PetscQuadratureGetData(q, NULL, NULL, &npoints, &points, &weights)); in testDerivativesLegendre()
59 PetscCall(PetscMalloc1(npoints * Np, &proj)); in testDerivativesLegendre()
60 PetscCall(PetscDTPKDEvalJet(dim, npoints, points, deg, 0, proj)); in testDerivativesLegendre()
62 for (j = 0; j < npoints; j++) proj[i * npoints + j] *= weights[j]; in testDerivativesLegendre()
77 PetscCall(PetscMalloc1(npoints, &points1d)); in testDerivativesLegendre()
80 PetscCall(PetscMalloc1((deg + 1) * npoints, &B[d])); in testDerivativesLegendre()
82 for (i = 0; i < npoints; i++) points1d[i] = points[i * dim + d]; in testDerivativesLegendre()
83 PetscCall(PetscDTLegendreEval(npoints, points1d, deg + 1, degrees, B[d], NULL, NULL)); in testDerivativesLegendre()
86 for (i = 0; i < npoints; i++) { in testDerivativesLegendre()
97 for (j = 0; j < Np; j++) pkd_coeffs[j] += proj[j * npoints + i] * val; in testDerivativesLegendre()