From: Luca Heltai Date: Sat, 11 Sep 2010 16:14:40 +0000 (+0000) Subject: Added assertion in QGuassOneOverR X-Git-Tag: v8.0.0~5543 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=67a0ab55c532d05507ff138b640adc82b1d1c297;p=dealii.git Added assertion in QGuassOneOverR git-svn-id: https://svn.dealii.org/trunk@21930 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/base/source/quadrature_lib.cc b/deal.II/base/source/quadrature_lib.cc index 0ee73d7fbf..caf20b2350 100644 --- a/deal.II/base/source/quadrature_lib.cc +++ b/deal.II/base/source/quadrature_lib.cc @@ -1075,6 +1075,14 @@ QGaussOneOverR<2>::QGaussOneOverR(const unsigned int n, const bool factor_out_singularity) : Quadrature<2>(2*n*n) { + // This version of the constructor + // works only for the 4 + // vertices. If you need a more + // general one, you should use the + // one with the Point<2> in the + // constructor. + Assert(vertex_index <4, ExcIndexRange(vertex_index, 0, 4)); + // Start with the gauss quadrature formula on the (u,v) reference // element. QGauss<2> gauss(n);