From 55d8f9b313ebe739bf3a754659b446cc5dc0a011 Mon Sep 17 00:00:00 2001 From: wolf Date: Mon, 24 Jul 2000 08:25:03 +0000 Subject: [PATCH] Remove QGauss8, since it is broken. git-svn-id: https://svn.dealii.org/trunk@3199 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/base/include/base/quadrature_lib.h | 20 -------- deal.II/base/source/quadrature_lib.cc | 54 ---------------------- 2 files changed, 74 deletions(-) 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 }; -/** - * 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 -class QGauss8 : public Quadrature -{ - 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(xpts[i]); - weights[i] = wts[i]; - }; -}; - template <> QMidpoint<1>::QMidpoint () : @@ -374,10 +326,6 @@ template QGauss7::QGauss7 () : Quadrature (QGauss7(), QGauss7<1>()){}; -template -QGauss8::QGauss8 () : - Quadrature (QGauss8(), QGauss8<1>()){}; - template QMidpoint::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>; -- 2.39.5