]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Test project_to_surface without boundary objects.
authorDavid Wells <wellsd2@rpi.edu>
Sat, 12 Aug 2017 15:50:37 +0000 (11:50 -0400)
committerDavid Wells <wellsd2@rpi.edu>
Sun, 8 Oct 2017 23:22:35 +0000 (19:22 -0400)
Since GridTools::project_to_object now calls our implementation of
project_to_surface directly we can fix up these tests to just call the relevant
function.

tests/numerics/project_to_surface_01.cc
tests/numerics/project_to_surface_02.cc
tests/numerics/project_to_surface_03.cc

index 8b25a5cab1c82c6728eb5306bb8c9ef6d82ecbc9..39859c429d0d7bed76ee2ceabfb60a69d0cb4869 100644 (file)
 
 
 
-// test StraightBoundary::project_to_surface for lines
+// test GridTools::project_to_object for lines
 
 
 
 #include "../tests.h"
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/grid_generator.h>
@@ -82,7 +81,6 @@ void test ()
   deallog << "dim=" << dim << std::endl;
 
   Triangulation<dim> tria;
-  StraightBoundary<dim> boundary;
 
   for (unsigned int case_no=0; case_no<2; ++case_no)
     {
@@ -98,9 +96,9 @@ void test ()
         {
           deallog << "    Line " << e << ", projected point=";
           if (dim > 1)
-            deallog << boundary.project_to_surface (cell->line(e), trial_point);
+            deallog << GridTools::project_to_object (cell->line(e), trial_point);
           else
-            deallog << boundary.project_to_surface (cell, trial_point);
+            deallog << GridTools::project_to_object (cell, trial_point);
 
           deallog << "  (line is from ";
           if (dim > 1)
index 129b463ffc649cf16389c755106f7ec172c77afb..0a3739a8efbc9aa62a3d6e3293f85d0990d75707 100644 (file)
 
 
 
-// test StraightBoundary::project_to_surface for quads
+// test GridTools::project_to_object for quads
 
 
 
 #include "../tests.h"
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/grid_generator.h>
@@ -82,7 +81,6 @@ void test ()
   deallog << "dim=" << dim << std::endl;
 
   Triangulation<dim> tria;
-  StraightBoundary<dim> boundary;
 
   for (unsigned int case_no=0; case_no<2; ++case_no)
     {
@@ -102,7 +100,7 @@ void test ()
 
           deallog << "    Quad " << e << ", projected point=";
 
-          const Point<dim> p = boundary.project_to_surface (quad, trial_point);
+          const Point<dim> p = GridTools::project_to_object (quad, trial_point);
           deallog << p;
           deallog << "  (quad is from ";
           deallog << quad->vertex(0);
index 7907da9fb733053f5b2cd71ba76ce1844ff7bf20..ad96b69468d4a1981eddb3fd7908ac6b0ca995c3 100644 (file)
@@ -22,7 +22,6 @@
 #include "../tests.h"
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/grid_generator.h>
@@ -36,7 +35,6 @@ void test ()
   deallog << "dim=" << dim << std::endl;
 
   Triangulation<dim> tria;
-  StraightBoundary<dim> boundary;
 
   GridGenerator::hyper_cube(tria, 0, 1);
 
@@ -73,7 +71,7 @@ void test ()
 
           deallog << "    Quad " << e << ", projected point=";
 
-          const Point<dim> p = boundary.project_to_surface (quad, trial_point);
+          const Point<dim> p = GridTools::project_to_object (quad, trial_point);
           deallog << p;
           deallog << "  (quad is from ";
           deallog << quad->vertex(0);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.