From: hartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 22 Jan 2001 12:27:16 +0000 (+0000)
Subject: Sort the quadrature points in ascending order.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f5a46d9cd4d93abaf66dad07bf29a9efa65bb2d4;p=dealii-svn.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)) };