35 static PetscErrorCode CheckQuadrature_Basics(PetscInt npoints, PetscReal alpha, PetscReal beta, con…  in CheckQuadrature_Basics()  argument
41 …ta = %g, i = %" PetscInt_FMT ", x[i] = %g, x[i-1] = %g", npoints, (double)alpha, (double)beta, i, …  in CheckQuadrature_Basics()
44 …lpha = %g, beta = %g, i = %" PetscInt_FMT ", w[i] = %g", npoints, (double)alpha, (double)beta, i, …  in CheckQuadrature_Basics()
49 static PetscErrorCode CheckQuadrature(PetscInt npoints, PetscReal alpha, PetscReal beta, const Pets…  in CheckQuadrature()  argument
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()
70         I_exact = PetscPowReal(2.0, alpha + beta + 1.) / (2. * i + alpha + beta + 1.);  in CheckQuadrature()
72 …I_exact *= PetscExpReal(PetscLGamma(i + alpha + 1.) + PetscLGamma(i + beta + 1.) - (PetscLGamma(i …  in CheckQuadrature()
78           for (k = 0; k < ibeta; k++) I_exact *= (i + 1. + k) / (i + alpha + 1. + k);  in CheckQuadrature()
82         PetscCall(PetscDTJacobiNorm(alpha, beta, i, &norm));  in CheckQuadrature()
90 …etscInt_FMT ", j %" PetscInt_FMT ", exact %g, err %g\n", npoints, (double)alpha, (double)beta, i, …  in CheckQuadrature()
91 …ule with alpha = %g, beta = %g: exact %g, err %g", i, j, npoints, (double)alpha, (double)beta, (do…  in CheckQuadrature()
98 static PetscErrorCode CheckJacobiQuadrature(PetscInt npoints, PetscReal alpha, PetscReal beta, quad…  in CheckJacobiQuadrature()  argument
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()
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()
127 …g, i %" PetscInt_FMT ", xdiff/xtol %g, wdiff/wtol %g\n", npoints, (double)alpha, (double)beta, i, …  in CheckJacobiQuadrature()
128 …pha = %g, beta = %g, i = %" PetscInt_FMT ", xdiff = %g", npoints, (double)alpha, (double)beta, i, …  in CheckJacobiQuadrature()
129 …pha = %g, beta = %g, i = %" PetscInt_FMT ", wdiff = %g", npoints, (double)alpha, (double)beta, i, …  in CheckJacobiQuadrature()