8 static PetscErrorCode CheckPoints(const char *name, PetscInt npoints, const PetscReal *points, Pets…  in CheckPoints()  argument
14   PetscCall(PetscMalloc3(npoints * ndegrees, &B, npoints * ndegrees, &D, npoints * ndegrees, &D2));  in CheckPoints()
15   PetscCall(PetscDTLegendreEval(npoints, points, ndegrees, degrees, B, D, D2));  in CheckPoints()
17   for (i = 0; i < npoints; i++) {  in CheckPoints()
35 static PetscErrorCode CheckQuadrature_Basics(PetscInt npoints, PetscReal alpha, PetscReal beta, con…  in CheckQuadrature_Basics()  argument
40   for (i = 1; i < npoints; i++) {  in CheckQuadrature_Basics()
41 …, alpha = %g, beta = %g, i = %" PetscInt_FMT ", x[i] = %g, x[i-1] = %g", npoints, (double)alpha, (…  in CheckQuadrature_Basics()
43   for (i = 0; i < npoints; i++) {  in CheckQuadrature_Basics()
44 …_FMT " points, alpha = %g, beta = %g, i = %" PetscInt_FMT ", w[i] = %g", npoints, (double)alpha, (…  in CheckQuadrature_Basics()
49 static PetscErrorCode CheckQuadrature(PetscInt npoints, PetscReal alpha, PetscReal beta, const Pets…  in CheckQuadrature()  argument
57   PetscCall(PetscMalloc2(npoints, &Pi, npoints, &Pj));  in CheckQuadrature()
59     PetscCall(PetscDTJacobiEval(npoints, alpha, beta, x, 1, &i, Pi, NULL, NULL));  in CheckQuadrature()
64       PetscCall(PetscDTJacobiEval(npoints, alpha, beta, x, 1, &j, Pj, NULL, NULL));  in CheckQuadrature()
88       for (k = 0; k < npoints; k++) I_quad += w[k] * (Pi[k] * Pj[k]);  in CheckQuadrature()
90 … beta %g, i %" PetscInt_FMT ", j %" PetscInt_FMT ", exact %g, err %g\n", npoints, (double)alpha, (…  in CheckQuadrature()
91 …nt_FMT " point rule with alpha = %g, beta = %g: exact %g, err %g", i, j, npoints, (double)alpha, (…  in CheckQuadrature()
98 static PetscErrorCode CheckJacobiQuadrature(PetscInt npoints, PetscReal alpha, PetscReal beta, quad…  in CheckJacobiQuadrature()  argument
103   PetscCall(PetscMalloc2(npoints, &x, npoints, &w));  in CheckJacobiQuadrature()
104   PetscCall((*func)(npoints, -1., 1., alpha, beta, x, w));  in CheckJacobiQuadrature()
105   PetscCall(CheckQuadrature_Basics(npoints, alpha, beta, x, w));  in CheckJacobiQuadrature()
106   PetscCall(CheckQuadrature(npoints, alpha, beta, x, w, nexact));  in CheckJacobiQuadrature()
116     PetscCall(PetscMalloc2(npoints, &x2, npoints, &w2));  in CheckJacobiQuadrature()
117     PetscCall((*func)(npoints, -1., 1., alpha, beta, x2, w2));  in CheckJacobiQuadrature()
118     PetscCall(CheckQuadrature_Basics(npoints, alpha, beta, x2, w2));  in CheckJacobiQuadrature()
119     PetscCall(CheckQuadrature(npoints, alpha, beta, x2, w2, nexact));  in CheckJacobiQuadrature()
120     for (i = 0; i < npoints; i++) {  in CheckJacobiQuadrature()
127 …alpha %g, beta %g, i %" PetscInt_FMT ", xdiff/xtol %g, wdiff/wtol %g\n", npoints, (double)alpha, (…  in CheckJacobiQuadrature()
128 …FMT " points, alpha = %g, beta = %g, i = %" PetscInt_FMT ", xdiff = %g", npoints, (double)alpha, (…  in CheckJacobiQuadrature()
129 …FMT " points, alpha = %g, beta = %g, i = %" PetscInt_FMT ", wdiff = %g", npoints, (double)alpha, (…  in CheckJacobiQuadrature()
142   PetscInt  degrees[1000], ndegrees, npoints, two;  in main()  local
158     npoints   = 1000;  in main()
162 …tscOptionsRealArray("-points", "list of points at which to evaluate", "", points, &npoints, &flg));  in main()
164     if (!flg) npoints = 3;  in main()
181   PetscCall(CheckPoints("User-provided points", npoints, points, ndegrees, degrees));  in main()
183   PetscCall(PetscDTGaussQuadrature(npoints, interval[0], interval[1], points, weights));  in main()
185   PetscCall(PetscRealView(npoints, weights, PETSC_VIEWER_STDOUT_WORLD));  in main()
192     for (i = 0; i < npoints; i++) {  in main()
202   PetscCall(CheckPoints("Gauss points", npoints, points, ndegrees, degrees));  in main()