Home
last modified time | relevance | path

Searched refs:I_exact (Results 1 – 2 of 2) sorted by relevance

/petsc/src/dm/dt/tests/
H A Dex1.c62 PetscReal I_exact = 0.; in CheckQuadrature() local
70 I_exact = PetscPowReal(2.0, alpha + beta + 1.) / (2. * i + alpha + beta + 1.); in CheckQuadrature()
72I_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()
83 norm2diff = PetscAbsReal(norm * norm - I_exact); in CheckQuadrature()
84 …PetscCheck(norm2diff <= eps * I_exact, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Jacobi norm error %g", (d… in CheckQuadrature()
86 tol = eps * I_exact; in CheckQuadrature()
89 err = PetscAbsReal(I_exact - I_quad); in CheckQuadrature()
90 …", exact %g, err %g\n", npoints, (double)alpha, (double)beta, i, j, (double)I_exact, (double)err)); in CheckQuadrature()
91 …= %g: exact %g, err %g", i, j, npoints, (double)alpha, (double)beta, (double)I_exact, (double)err); in CheckQuadrature()
H A Dex15.c27 PetscReal I_exact = (i == j) ? 1.0 : 0.0; in testQuadrature() local
32 err = PetscAbsReal(I_exact - I_quad); in testQuadrature()