From: Wolfgang Bangerth Date: Thu, 5 Feb 2015 02:57:45 +0000 (-0600) Subject: Fix up places where we use Point for directions when we should be X-Git-Tag: v8.3.0-rc1~494^2~12 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=44ac8875d5f41d654eb077708a2c36286fd887f6;p=dealii.git Fix up places where we use Point for directions when we should be using Tensor<1,dim>. --- 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