From: Ralf Hartmann Date: Mon, 22 Jan 2001 12:27:16 +0000 (+0000) Subject: Sort the quadrature points in ascending order. X-Git-Tag: v8.0.0~19796 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d5bace63fc28c04024b8faba55f900d6ad3f12f5;p=dealii.git Sort the quadrature points in ascending order. git-svn-id: https://svn.dealii.org/trunk@3752 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/base/source/quadrature_lib.cc b/deal.II/base/source/quadrature_lib.cc index 3ca857bad3..dc3ef1282a 100644 --- a/deal.II/base/source/quadrature_lib.cc +++ b/deal.II/base/source/quadrature_lib.cc @@ -77,13 +77,13 @@ QGauss4<1>::QGauss4 () : Quadrature<1> (4) { // points on [-1,1] - static const double xpts_normal[] = { -sqrt(1./7.*(3-4*sqrt(0.3))), - -sqrt(1./7.*(3+4*sqrt(0.3))), + static const double xpts_normal[] = { -sqrt(1./7.*(3+4*sqrt(0.3))), + -sqrt(1./7.*(3-4*sqrt(0.3))), +sqrt(1./7.*(3-4*sqrt(0.3))), +sqrt(1./7.*(3+4*sqrt(0.3))) }; // weights on [-1,1] - static const double wts_normal[] = { 1./2. + 1./12.*sqrt(10./3.), - 1./2. - 1./12.*sqrt(10./3.), + static const double wts_normal[] = { 1./2. - 1./12.*sqrt(10./3.), + 1./2. + 1./12.*sqrt(10./3.), 1./2. + 1./12.*sqrt(10./3.), 1./2. - 1./12.*sqrt(10./3.) }; @@ -110,14 +110,14 @@ QGauss5<1>::QGauss5 () : Quadrature<1> (5) { // points on [-1,1] - static const double xpts_normal[] = { -sqrt(1./9.*(5.-2*sqrt(10./7.))), - -sqrt(1./9.*(5.+2*sqrt(10./7.))), + static const double xpts_normal[] = { -sqrt(1./9.*(5.+2*sqrt(10./7.))), + -sqrt(1./9.*(5.-2*sqrt(10./7.))), 0, +sqrt(1./9.*(5.-2*sqrt(10./7.))), +sqrt(1./9.*(5.+2*sqrt(10./7.))) }; // weights on [-1,1] - static const double wts_normal[] = { 0.3*(-0.7+5.*sqrt(0.7))/(-2.+5.*sqrt(0.7)), - 0.3*(+0.7+5.*sqrt(0.7))/(+2.+5.*sqrt(0.7)), + static const double wts_normal[] = { 0.3*(+0.7+5.*sqrt(0.7))/(+2.+5.*sqrt(0.7)), + 0.3*(-0.7+5.*sqrt(0.7))/(-2.+5.*sqrt(0.7)), 128./225., 0.3*(-0.7+5.*sqrt(0.7))/(-2.+5.*sqrt(0.7)), 0.3*(+0.7+5.*sqrt(0.7))/(+2.+5.*sqrt(0.7)) };