From ae9199e0ff1856b3c66e473749c4e70da2178ad1 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Thu, 30 Jul 2009 18:48:05 +0000 Subject: [PATCH] First step in implementing project_to_surface, git-svn-id: https://svn.dealii.org/trunk@19143 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/grid/tria_boundary.cc | 44 ++++++++++++++------ 1 file changed, 32 insertions(+), 12 deletions(-) diff --git a/deal.II/deal.II/source/grid/tria_boundary.cc b/deal.II/deal.II/source/grid/tria_boundary.cc index c89556a76e..dd6ba2ef5d 100644 --- a/deal.II/deal.II/source/grid/tria_boundary.cc +++ b/deal.II/deal.II/source/grid/tria_boundary.cc @@ -167,10 +167,15 @@ template Point Boundary:: project_to_surface (const typename Triangulation::line_iterator &, - const Point &) const + const Point &trial_point) const { - Assert (false, ExcPureFunctionCalled()); - return Point(); + if (spacedim <= 1) + return trial_point; + else + { + Assert (false, ExcPureFunctionCalled()); + return Point(); + } } @@ -179,10 +184,15 @@ template Point Boundary:: project_to_surface (const typename Triangulation::quad_iterator &, - const Point &) const + const Point &trial_point) const { - Assert (false, ExcPureFunctionCalled()); - return Point(); + if (spacedim <= 2) + return trial_point; + else + { + Assert (false, ExcPureFunctionCalled()); + return Point(); + } } @@ -489,10 +499,15 @@ template Point StraightBoundary:: project_to_surface (const typename Triangulation::line_iterator &, - const Point &) const + const Point &trial_point) const { - Assert (false, ExcPureFunctionCalled()); - return Point(); + if (spacedim <= 1) + return trial_point; + else + { + Assert (false, ExcPureFunctionCalled()); + return Point(); + } } @@ -501,10 +516,15 @@ template Point StraightBoundary:: project_to_surface (const typename Triangulation::quad_iterator &, - const Point &) const + const Point &trial_point) const { - Assert (false, ExcPureFunctionCalled()); - return Point(); + if (spacedim <= 2) + return trial_point; + else + { + Assert (false, ExcPureFunctionCalled()); + return Point(); + } } -- 2.39.5