]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Added assertion in QGuassOneOverR
authorLuca Heltai <luca.heltai@sissa.it>
Sat, 11 Sep 2010 16:14:40 +0000 (16:14 +0000)
committerLuca Heltai <luca.heltai@sissa.it>
Sat, 11 Sep 2010 16:14:40 +0000 (16:14 +0000)
git-svn-id: https://svn.dealii.org/trunk@21930 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/base/source/quadrature_lib.cc

index 0ee73d7fbf9ae74b3f622d0c3dfe0f35f5e40724..caf20b2350034e70bd008553ae231d32ee49b6b1 100644 (file)
@@ -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);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.