]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix up places where we use Point<dim> for directions when we should be
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 5 Feb 2015 03:38:44 +0000 (21:38 -0600)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 5 Feb 2015 03:38:44 +0000 (21:38 -0600)
using Tensor<1,dim>.

source/base/quadrature_lib.cc

index dde1f1b9fa96f4fff5ecd9beb8279321a9cfc8d1..102d7b7fbb8fdc5b876b0649255a446641e0de09 100644 (file)
@@ -786,7 +786,7 @@ unsigned int QGaussOneOverR<2>::quad_size(const Point<2> singularity,
     if ( ( std::abs(singularity[i]  ) < eps ) ||
          ( std::abs(singularity[i]-1) < eps ) )
       on_edge = true;
-  if (on_edge && (std::abs( (singularity-Point<2>(.5, .5)).square()-.5)
+  if (on_edge && (std::abs( (singularity-Point<2>(.5, .5)).norm_square()-.5)
                   < eps) )
     on_vertex = true;
   if (on_vertex) return (2*n*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.