From b8c5a939075b35c03fcda671549ebd66e690a945 Mon Sep 17 00:00:00 2001 From: guido Date: Mon, 12 Feb 2001 13:35:50 +0000 Subject: [PATCH] new Gauss quadratures, larger numerical errors on SUN? git-svn-id: https://svn.dealii.org/trunk@3912 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/base/quadrature_test.cc | 28 +++++++++++++++------------- 1 file changed, 15 insertions(+), 13 deletions(-) diff --git a/tests/base/quadrature_test.cc b/tests/base/quadrature_test.cc index 2af6236bd7..83cabb0ac7 100644 --- a/tests/base/quadrature_test.cc +++ b/tests/base/quadrature_test.cc @@ -135,7 +135,8 @@ check_faces (vector*>& quadratures, bool sub) ++i; quadrature_int=0; - // Check the polynomial x^i*y^i + // Check the polynomial + // x^i*y^i*z^i for (unsigned int x=0; x*>& quadratures, bool sub) quadrature_int+=f*weights[x]; } - // the exact integral is 1/(i+1) - switch (dim) - { - case 2: - exact_int = 2 * (sub ? 2:1) / pow(i+1,dim-1); - break; - case 3: - exact_int = 3 * (sub ? 4:1) / pow(i+1,dim-1); - break; - } + // the exact integral is + // 1/(i+1)^(dim-1) + switch (dim) + { + case 2: + exact_int = 2 * (sub ? 2:1) / (double) (i+1); + break; + case 3: + exact_int = 3 * (sub ? 4:1) / (double) (i+1)/(i+1); + break; + } err = fabs(quadrature_int-exact_int); } - while (err<1e-16); + while (err<5e-15); // Uncomment here for testing -// deallog << " (Int " << quadrature_int << ',' << exact_int << ")"; + // deallog << " (Int " << quadrature_int << '-' << exact_int << '=' << err << ")"; deallog << " is exact for polynomials of degree " << i-1 << endl; } deallog.pop(); -- 2.39.5