From e75eabc217b85c364f51f376a148ccc5b2f71abf Mon Sep 17 00:00:00 2001
From: heltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Tue, 4 Feb 2014 00:41:22 +0000
Subject: [PATCH] Fixed interfaces in project_to_surface.

git-svn-id: https://svn.dealii.org/branches/branch_manifold_id@32399 0785d39b-7218-0410-832d-ea1e28bc413d
---
 tests/deal.II/project_to_surface_01.cc | 19 ++++++++++++++-----
 tests/deal.II/project_to_surface_02.cc |  8 ++++++--
 tests/deal.II/project_to_surface_03.cc |  7 ++++++-
 3 files changed, 26 insertions(+), 8 deletions(-)

diff --git a/tests/deal.II/project_to_surface_01.cc b/tests/deal.II/project_to_surface_01.cc
index 7fda38c7ae..f6775cc67e 100644
--- a/tests/deal.II/project_to_surface_01.cc
+++ b/tests/deal.II/project_to_surface_01.cc
@@ -97,14 +97,23 @@ void test ()
       Point<dim> trial_point;
       for (unsigned int i=0; i<dim; ++i)
         trial_point[i] = 1.5;
-
+      
+      std::vector<Point<dim> > vs(2);
+      
       for (unsigned int e=0; e<GeometryInfo<dim>::lines_per_cell; ++e)
         {
           deallog << "    Line " << e << ", projected point=";
-          if (dim > 1)
-            deallog << boundary.project_to_surface (cell->line(e), trial_point);
-          else
-            deallog << boundary.project_to_surface (cell, trial_point);
+	  vs.resize(2);
+          if (dim > 1) {
+	    vs[0] = cell->line(e)->vertex(0);
+	    vs[1] = cell->line(e)->vertex(1);
+            deallog << boundary.project_to_manifold (vs, trial_point);
+	  }
+          else {
+	    vs[0] = cell->vertex(0);
+	    vs[1] = cell->vertex(1);
+            deallog << boundary.project_to_manifold (vs, trial_point);
+	  }
 
           deallog << "  (line is from ";
           if (dim > 1)
diff --git a/tests/deal.II/project_to_surface_02.cc b/tests/deal.II/project_to_surface_02.cc
index 9e42fc3a6c..b098e96d90 100644
--- a/tests/deal.II/project_to_surface_02.cc
+++ b/tests/deal.II/project_to_surface_02.cc
@@ -97,7 +97,9 @@ void test ()
       Point<dim> trial_point;
       for (unsigned int i=0; i<dim; ++i)
         trial_point[i] = 1.5;
-
+      
+      std::vector<Point<dim> > vs(GeometryInfo<dim>::vertices_per_cell);
+      
       for (unsigned int e=0; e<GeometryInfo<dim>::quads_per_cell; ++e)
         {
           const typename Triangulation<dim>::quad_iterator
@@ -105,8 +107,10 @@ void test ()
                   *reinterpret_cast<const typename Triangulation<dim>::quad_iterator *>(&cell));
 
           deallog << "    Quad " << e << ", projected point=";
+	  for(unsigned int i=0; i<GeometryInfo<dim>::vertices_per_cell; ++i)
+	    vs[i] = quad->vertex(i);
 
-          const Point<dim> p = boundary.project_to_surface (quad, trial_point);
+          const Point<dim> p = boundary.project_to_manifold (vs, trial_point);
           deallog << p;
           deallog << "  (quad is from ";
           deallog << quad->vertex(0);
diff --git a/tests/deal.II/project_to_surface_03.cc b/tests/deal.II/project_to_surface_03.cc
index f99a6de525..9b488cf2b9 100644
--- a/tests/deal.II/project_to_surface_03.cc
+++ b/tests/deal.II/project_to_surface_03.cc
@@ -69,6 +69,8 @@ void test ()
 
       deallog << "Trial point = " << trial_point << std::endl;
 
+      std::vector<Point<dim> > vs(GeometryInfo<dim>::vertices_per_cell);
+
       for (unsigned int e=0; e<GeometryInfo<dim>::quads_per_cell; ++e)
         {
           const typename Triangulation<dim>::quad_iterator
@@ -76,8 +78,11 @@ void test ()
                   *reinterpret_cast<const typename Triangulation<dim>::quad_iterator *>(&cell));
 
           deallog << "    Quad " << e << ", projected point=";
+	  
+	  for(unsigned int i=0; i<GeometryInfo<dim>::vertices_per_cell; ++i)
+	    vs[i] = quad->vertex(i);
 
-          const Point<dim> p = boundary.project_to_surface (quad, trial_point);
+          const Point<dim> p = boundary.project_to_surface (vs, trial_point);
           deallog << p;
           deallog << "  (quad is from ";
           deallog << quad->vertex(0);
-- 
2.39.5