Fix up places where we use Point<dim> for directions when we should be

using Tensor<1,dim>.

diff --git a/source/base/quadrature_lib.cc b/source/base/quadrature_lib.cc
index dde1f1b9fa96f4fff5ecd9beb8279321a9cfc8d1..102d7b7fbb8fdc5b876b0649255a446641e0de09 100644 (file)
--- a/source/base/quadrature_lib.cc
+++ b/source/base/quadrature_lib.cc
@@ -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);