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 3ca857bad3baece573966842735dacc34c8b8e60..dc3ef1282a4413aa70598d7fa9adabb8648bb003 100644 (file)
--- 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)) };