From 67a0ab55c532d05507ff138b640adc82b1d1c297 Mon Sep 17 00:00:00 2001 From: Luca Heltai Date: Sat, 11 Sep 2010 16:14:40 +0000 Subject: [PATCH] Added assertion in QGuassOneOverR git-svn-id: https://svn.dealii.org/trunk@21930 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/base/source/quadrature_lib.cc | 8 ++++++++ 1 file changed, 8 insertions(+) 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); -- 2.39.5