From: Guido Kanschat Date: Mon, 17 Oct 2005 05:45:47 +0000 (+0000) Subject: allow Gauss formulae with zero points X-Git-Tag: v8.0.0~12997 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=8b945c3585b099335cb9daf29122fa543ec40d95;p=dealii.git allow Gauss formulae with zero points git-svn-id: https://svn.dealii.org/trunk@11606 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/base/source/quadrature.cc b/deal.II/base/source/quadrature.cc index 08edb028e1..7d194347d8 100644 --- a/deal.II/base/source/quadrature.cc +++ b/deal.II/base/source/quadrature.cc @@ -140,13 +140,16 @@ Quadrature::Quadrature (const SubQuadrature &q1, }; #ifdef DEBUG - double sum = 0; - for (unsigned int i=0; i0.999999) && (sum<1.000001), ExcInternalError()); + if (n_quadrature_points > 0) + { + double sum = 0; + for (unsigned int i=0; i0.999999) && (sum<1.000001), ExcInternalError()); + } #endif } diff --git a/deal.II/base/source/quadrature_lib.cc b/deal.II/base/source/quadrature_lib.cc index ab876dbc64..2e267d5c7d 100644 --- a/deal.II/base/source/quadrature_lib.cc +++ b/deal.II/base/source/quadrature_lib.cc @@ -51,6 +51,9 @@ QGauss<1>::QGauss (const unsigned int n) : Quadrature<1> (n) { + if (n == 0) + return; + const unsigned int m = (n+1)/2; // tolerance for the Newton