From 8b945c3585b099335cb9daf29122fa543ec40d95 Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Mon, 17 Oct 2005 05:45:47 +0000 Subject: [PATCH] allow Gauss formulae with zero points git-svn-id: https://svn.dealii.org/trunk@11606 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/base/source/quadrature.cc | 17 ++++++++++------- deal.II/base/source/quadrature_lib.cc | 3 +++ 2 files changed, 13 insertions(+), 7 deletions(-) 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 -- 2.39.5