]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Sort the quadrature points in ascending order.
authorRalf Hartmann <Ralf.Hartmann@dlr.de>
Mon, 22 Jan 2001 12:27:16 +0000 (12:27 +0000)
committerRalf Hartmann <Ralf.Hartmann@dlr.de>
Mon, 22 Jan 2001 12:27:16 +0000 (12:27 +0000)
git-svn-id: https://svn.dealii.org/trunk@3752 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/base/source/quadrature_lib.cc

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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.