From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Mon, 24 Jul 2000 08:25:03 +0000 (+0000)
Subject: Remove QGauss8, since it is broken.
X-Git-Tag: v8.0.0~20227
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=89b075970e896207424dc595740ec1ed6a5c3147;p=dealii.git

Remove QGauss8, since it is broken.


git-svn-id: https://svn.dealii.org/trunk@3199 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/base/include/base/quadrature_lib.h b/deal.II/base/include/base/quadrature_lib.h
index a2868f65e3..65dd7a7129 100644
--- a/deal.II/base/include/base/quadrature_lib.h
+++ b/deal.II/base/include/base/quadrature_lib.h
@@ -117,26 +117,6 @@ class QGauss7 : public Quadrature<dim>
 };
 
 
-/**
- *  8-Point-Gauss quadrature formula, exact for polynomials of degree 15.
- *  I have not found explicite
- *  representations of the zeros of the Legendre functions of sixth
- *  and higher degree. If anyone finds them, please replace the existing
- *  numbers by these expressions.
- *
- *  Reference: J. E. Akin: "Application and Implementation of Finite
- *  Element Methods"
- *  For a comprehensive list of Gaussian quadrature formulae, see also:
- *  A. H. Strout, D. Secrest: "Gaussian Quadrature Formulas" 
- */
-template <int dim>
-class QGauss8 : public Quadrature<dim>
-{
-  public:
-    QGauss8 ();
-};
-
-
 /**
  * Midpoint quadrature rule, exact for linear polynomials.
  */
diff --git a/deal.II/base/source/quadrature_lib.cc b/deal.II/base/source/quadrature_lib.cc
index afe9e24f30..0bcec6fe37 100644
--- a/deal.II/base/source/quadrature_lib.cc
+++ b/deal.II/base/source/quadrature_lib.cc
@@ -228,54 +228,6 @@ QGauss7<1>::QGauss7 () :
 };
 
 
-template <>
-QGauss8<1>::QGauss8 () :
-		Quadrature<1> (8)
-{
-				   // points on [-1,1]
-  static const double xpts_normal[] = { -0.960289856497536,
-					-0.796666477413627,
-					-0.525532409916329,
-					-0.183434642495650,
-					+0.183434642495650,
-					+0.525532409916329,
-					+0.796666477413627,
-					+0.960289856497536 };
-				   // weights on [-1,1]
-  static const double wts_normal[]  = { 0.101228536200376,
-					0.222381034453374,
-					0.313706645877887,
-					0.362683783378362,
-					0.362683783378362,
-					0.313706645877887,
-					0.222381034453374,
-					0.101228536200376  };
-
-				   // points and weights on [0,1]
-  static const double xpts[] = { (xpts_normal[0]+1)/2.,
-				 (xpts_normal[1]+1)/2.,
-				 (xpts_normal[2]+1)/2.,
-				 (xpts_normal[3]+1)/2.,
-				 (xpts_normal[4]+1)/2.,
-				 (xpts_normal[5]+1)/2.,
-				 (xpts_normal[6]+1)/2.,
-				 (xpts_normal[7]+1)/2. };
-  static const double wts[]  = { wts_normal[0]/2.,
-				 wts_normal[1]/2.,
-				 wts_normal[2]/2.,
-				 wts_normal[3]/2.,
-				 wts_normal[4]/2.,
-				 wts_normal[5]/2.,
-				 wts_normal[6]/2.,
-				 wts_normal[7]/2. };
-
-  for (unsigned int i=0; i<n_quadrature_points; ++i) 
-    {
-      quadrature_points[i] = Point<1>(xpts[i]);
-      weights[i] = wts[i];
-    };
-};
-
 
 template <>
 QMidpoint<1>::QMidpoint () :
@@ -374,10 +326,6 @@ template <int dim>
 QGauss7<dim>::QGauss7 () :
 		Quadrature<dim> (QGauss7<dim-1>(), QGauss7<1>()){};
 
-template <int dim>
-QGauss8<dim>::QGauss8 () :
-		Quadrature<dim> (QGauss8<dim-1>(), QGauss8<1>()){};
-
 
 template <int dim>
 QMidpoint<dim>::QMidpoint () :
@@ -408,7 +356,6 @@ template class QGauss4<2>;
 template class QGauss5<2>;
 template class QGauss6<2>;
 template class QGauss7<2>;
-template class QGauss8<2>;
 template class QMidpoint<2>;
 template class QTrapez<2>;
 template class QSimpson<2>;
@@ -421,7 +368,6 @@ template class QGauss4<3>;
 template class QGauss5<3>;
 template class QGauss6<3>;
 template class QGauss7<3>;
-template class QGauss8<3>;
 template class QMidpoint<3>;
 template class QTrapez<3>;
 template class QSimpson<3>;