From 44ac8875d5f41d654eb077708a2c36286fd887f6 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Wed, 4 Feb 2015 20:57:45 -0600 Subject: [PATCH] Fix up places where we use Point for directions when we should be using Tensor<1,dim>. --- source/grid/tria_boundary_lib.cc | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/source/grid/tria_boundary_lib.cc b/source/grid/tria_boundary_lib.cc index f91640abb4..5eae3f2bce 100644 --- a/source/grid/tria_boundary_lib.cc +++ b/source/grid/tria_boundary_lib.cc @@ -1,6 +1,6 @@ // --------------------------------------------------------------------- // -// Copyright (C) 1999 - 2014 by the deal.II authors +// Copyright (C) 1999 - 2015 by the deal.II authors // // This file is part of the deal.II library. // @@ -74,8 +74,8 @@ get_new_point_on_line (const typename Triangulation::line_iterator // have to take into account the // offset point_on_axis and the // direction of the axis - const Point vector_from_axis = (middle-point_on_axis) - - ((middle-point_on_axis) * direction) * direction; + const Tensor<1,spacedim> vector_from_axis = (middle-point_on_axis) - + ((middle-point_on_axis) * direction) * direction; // scale it to the desired length // and put everything back // together, unless we have a point @@ -100,8 +100,8 @@ get_new_point_on_quad (const Triangulation<3>::quad_iterator &quad) const // same algorithm as above const unsigned int spacedim = 3; - const Point vector_from_axis = (middle-point_on_axis) - - ((middle-point_on_axis) * direction) * direction; + const Tensor<1,spacedim> vector_from_axis = (middle-point_on_axis) - + ((middle-point_on_axis) * direction) * direction; if (vector_from_axis.norm() <= 1e-10 * middle.norm()) return middle; else @@ -119,8 +119,8 @@ get_new_point_on_quad (const Triangulation<2,3>::quad_iterator &quad) const // same algorithm as above const unsigned int spacedim = 3; - const Point vector_from_axis = (middle-point_on_axis) - - ((middle-point_on_axis) * direction) * direction; + const Tensor<1,spacedim> vector_from_axis = (middle-point_on_axis) - + ((middle-point_on_axis) * direction) * direction; if (vector_from_axis.norm() <= 1e-10 * middle.norm()) return middle; else -- 2.39.5