]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Replaced and regexped most tests from Boundary to Manifold.
authorheltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 14 Nov 2013 16:15:54 +0000 (16:15 +0000)
committerheltai <heltai@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 14 Nov 2013 16:15:54 +0000 (16:15 +0000)
git-svn-id: https://svn.dealii.org/branches/branch_manifold_id@31654 0785d39b-7218-0410-832d-ea1e28bc413d

153 files changed:
tests/bits/cone_01.cc
tests/bits/cylinder_01.cc
tests/bits/cylinder_02.cc
tests/bits/cylinder_03.cc
tests/bits/cylinder_04.cc
tests/bits/distorted_cells_03.cc
tests/bits/distorted_cells_04.cc
tests/bits/distorted_cells_05.cc
tests/bits/distorted_cells_06.cc
tests/bits/distorted_mapped_cells_02.cc
tests/bits/fe_field_function_04.cc
tests/bits/fe_field_function_04_vector.cc
tests/bits/fe_field_function_05.cc
tests/bits/fe_field_function_05_vector.cc
tests/bits/fe_field_function_06_vector.cc
tests/bits/fe_field_function_07_vector.cc
tests/bits/fe_field_function_08_vector.cc
tests/bits/joa_1.cc
tests/bits/step-10-high-order.cc
tests/bits/step-10.cc
tests/bits/step-11.cc
tests/bits/step-5.cc
tests/bits/step-6.cc
tests/codim_one/boundary_indicator_01.cc
tests/codim_one/extract_boundary_mesh_02.cc
tests/codim_one/extract_boundary_mesh_03.cc
tests/codim_one/extract_boundary_mesh_04.cc
tests/codim_one/grid_refinement.cc
tests/codim_one/mapping_01.cc
tests/codim_one/mapping_02.cc
tests/codim_one/torus_01.cc
tests/deal.II/boundaries.cc
tests/deal.II/dof_test.cc
tests/deal.II/fe_values_view_01.cc
tests/deal.II/fe_values_view_02.cc
tests/deal.II/fe_values_view_03.cc
tests/deal.II/fe_values_view_04.cc
tests/deal.II/fe_values_view_05.cc
tests/deal.II/fe_values_view_06.cc
tests/deal.II/fe_values_view_07.cc
tests/deal.II/fe_values_view_08.cc
tests/deal.II/fe_values_view_09.cc
tests/deal.II/fe_values_view_10.cc
tests/deal.II/fe_values_view_10_single_01.cc
tests/deal.II/fe_values_view_10_single_02.cc
tests/deal.II/fe_values_view_10_single_03.cc
tests/deal.II/fe_values_view_10_single_04.cc
tests/deal.II/fe_values_view_11.cc
tests/deal.II/fe_values_view_12.cc
tests/deal.II/fe_values_view_13.cc
tests/deal.II/fe_values_view_14.cc
tests/deal.II/fe_values_view_15.cc
tests/deal.II/fe_values_view_16.cc
tests/deal.II/fe_values_view_17.cc
tests/deal.II/fe_values_view_18.cc
tests/deal.II/fe_values_view_19.cc
tests/deal.II/fe_values_view_20.cc
tests/deal.II/fe_values_view_23.cc
tests/deal.II/fe_values_view_23_nonsymmetric.cc
tests/deal.II/fe_values_view_24.cc
tests/deal.II/fe_values_view_24_nonsymmetric.cc
tests/deal.II/grid_hyper_shell.cc
tests/deal.II/grid_hyper_shell_02.cc
tests/deal.II/grid_hyper_shell_03.cc
tests/deal.II/grid_hyper_shell_04.cc
tests/deal.II/grid_out.cc
tests/deal.II/grid_output_input.cc
tests/deal.II/grid_test.cc
tests/deal.II/grid_tools_03.cc
tests/deal.II/grid_transform.cc
tests/deal.II/grid_transform_3d.cc
tests/deal.II/memory_consumption_01.cc
tests/deal.II/no_flux_03.cc
tests/deal.II/no_flux_04.cc
tests/deal.II/no_flux_05.cc
tests/deal.II/no_flux_10.cc
tests/deal.II/no_flux_hp_03.cc
tests/deal.II/no_flux_hp_04.cc
tests/deal.II/no_flux_hp_05.cc
tests/deal.II/normal_flux_03.cc
tests/deal.II/project_01_curved_boundary.cc
tests/deal.II/project_to_surface_01.cc
tests/deal.II/project_to_surface_02.cc
tests/deal.II/project_to_surface_03.cc
tests/fe/cell_similarity_01.cc
tests/fe/cell_similarity_02.cc
tests/fe/cell_similarity_03.cc
tests/fe/cell_similarity_04.cc
tests/fe/cell_similarity_05.cc
tests/fe/cell_similarity_06.cc
tests/fe/cell_similarity_07.cc
tests/fe/cell_similarity_08.cc
tests/fe/cell_similarity_09.cc
tests/fe/cell_similarity_10.cc
tests/fe/cell_similarity_dgp_monomial_01.cc
tests/fe/cell_similarity_dgp_monomial_02.cc
tests/fe/cell_similarity_dgp_monomial_03.cc
tests/fe/cell_similarity_dgp_monomial_04.cc
tests/fe/cell_similarity_dgp_monomial_05.cc
tests/fe/cell_similarity_dgp_monomial_06.cc
tests/fe/cell_similarity_dgp_monomial_07.cc
tests/fe/cell_similarity_dgp_monomial_08.cc
tests/fe/cell_similarity_dgp_monomial_09.cc
tests/fe/cell_similarity_dgp_monomial_10.cc
tests/fe/cell_similarity_dgp_nonparametric_01.cc
tests/fe/cell_similarity_dgp_nonparametric_02.cc
tests/fe/cell_similarity_dgp_nonparametric_03.cc
tests/fe/cell_similarity_dgp_nonparametric_04.cc
tests/fe/cell_similarity_dgp_nonparametric_05.cc
tests/fe/cell_similarity_dgp_nonparametric_06.cc
tests/fe/cell_similarity_dgp_nonparametric_07.cc
tests/fe/cell_similarity_dgp_nonparametric_08.cc
tests/fe/cell_similarity_dgp_nonparametric_09.cc
tests/fe/cell_similarity_dgp_nonparametric_10.cc
tests/fe/jacobians.cc
tests/fe/mapping.cc
tests/fe/mapping_real_to_unit_02.cc
tests/fe/mapping_real_to_unit_q4_curved.cc
tests/fe/mapping_real_to_unit_q4_curved_codim.cc
tests/fe/mapping_real_to_unit_q4_sphere.cc
tests/fe/mapping_real_to_unit_q4_sphere_x.cc
tests/fe/mapping_real_to_unit_q4_sphere_y.cc
tests/fe/mapping_real_to_unit_q4_sphere_z.cc
tests/grid/grid_generator_01a.cc
tests/hp/project_01_curved_boundary.cc
tests/hp/step-10.cc
tests/hp/step-11.cc
tests/hp/step-11_compressed_set_sparsity.cc
tests/hp/step-5.cc
tests/hp/step-6.cc
tests/matrix_free/get_functions_circle.cc
tests/matrix_free/get_functions_gl.cc
tests/matrix_free/get_functions_mappingq.cc
tests/matrix_free/get_functions_q_hierarchical.cc
tests/matrix_free/get_functions_rect.cc
tests/matrix_free/get_values_plain.cc
tests/matrix_free/integrate_functions.cc
tests/matrix_free/integrate_functions_multife.cc
tests/matrix_free/matrix_vector_04.cc
tests/matrix_free/matrix_vector_07.cc
tests/matrix_free/matrix_vector_curl.cc
tests/matrix_free/matrix_vector_hp.cc
tests/matrix_free/matrix_vector_mg.cc
tests/matrix_free/matrix_vector_stokes_noflux.cc
tests/matrix_free/no_index_initialize.cc
tests/matrix_free/quadrature_points.cc
tests/mpi/multigrid_adaptive.cc
tests/mpi/multigrid_uniform.cc
tests/mpi/p4est_2d_constraintmatrix_03.cc
tests/mpi/p4est_2d_constraintmatrix_04.cc
tests/mpi/p4est_3d_constraintmatrix_03.cc
tests/multigrid/step-16-02.cc
tests/multigrid/step-16.cc

index b3f8ecf09500bddaffdcbe5159aa51ae78332eb5..ce008405b7565ac72a5b597ee2e734f30dafc2ee 100644 (file)
@@ -22,7 +22,7 @@
 #include <deal.II/base/logstream.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_generator.h>
@@ -45,7 +45,7 @@ void check ()
   p1[0] = -1;
   p2[0] = 1;
   static const ConeBoundary<dim> boundary (1, 0.5, p1, p2);
-  triangulation.set_boundary (0, boundary);
+  triangulation.set_manifold (0, boundary);
 
   triangulation.refine_global (2);
 
index 0a758c26b5e12273ccab9356850453da4ea78cc9..553a8d3bc1edb5ea24dcedc007eba02f8d9a24ae 100644 (file)
 
 
 
-// check the cells generated by the CylinderBoundary the default
+// check the cells generated by the CylinderManifold the default
 // axiparallel cylinder
 
 #include "../tests.h"
 #include <deal.II/base/logstream.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_generator.h>
@@ -41,8 +41,8 @@ void check ()
 {
   Triangulation<dim> triangulation;
   GridGenerator::cylinder (triangulation);
-  static const CylinderBoundary<dim> boundary (1,0);
-  triangulation.set_boundary (0, boundary);
+  static const CylinderManifold<dim> boundary (0);
+  triangulation.set_manifold (0, boundary);
   triangulation.refine_global (2);
 
   for (typename Triangulation<dim>::active_cell_iterator
index 00ad1f9c19d26a30db428520000f7fbc04b4e452..6430842c34b568a16ebfabd69f1ba667594ed917 100644 (file)
@@ -23,7 +23,7 @@
 #include <deal.II/base/logstream.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_generator.h>
@@ -54,8 +54,8 @@ void check ()
 
   GridTools::transform ((Point<dim> ( *)(const Point<dim> &))&rotate_to_y<dim>, triangulation);
 
-  static const CylinderBoundary<dim> boundary (1,1);
-  triangulation.set_boundary (0, boundary);
+  static const CylinderManifold<dim> boundary (1);
+  triangulation.set_manifold (0, boundary);
   triangulation.refine_global (2);
 
   for (typename Triangulation<dim>::active_cell_iterator
index b61ad237c3acad2bc8306b883e280cd717b62f5f..b7c93e5323a41b3f34a15e94fe16131737d19833 100644 (file)
@@ -23,7 +23,7 @@
 #include <deal.II/base/logstream.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_generator.h>
@@ -51,8 +51,8 @@ void check ()
 
   GridTools::transform ((Point<dim> ( *)(const Point<dim> &))&rotate_to_y<dim>, triangulation);
 
-  static const CylinderBoundary<dim> boundary (1,2);
-  triangulation.set_boundary (0, boundary);
+  static const CylinderManifold<dim> boundary (2);
+  triangulation.set_manifold (0, boundary);
   triangulation.refine_global (2);
 
   for (typename Triangulation<dim>::active_cell_iterator
index ea6ffa2f346556e19f7175509c27c4ba50b0d5d4..34572de73f8414fdddc9467e9f7328edde55b8b8 100644 (file)
@@ -23,7 +23,7 @@
 #include <deal.II/base/logstream.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_generator.h>
@@ -55,10 +55,9 @@ void check ()
 
   GridTools::transform ((Point<dim> ( *)(const Point<dim> &))&rotate_to_y<dim>, triangulation);
 
-  static const CylinderBoundary<dim> boundary (1,
-                                               Point<dim>(std::cos(xy_angle), std::sin(xy_angle), 0),
+  static const CylinderManifold<dim> boundary (Point<dim>(std::cos(xy_angle), std::sin(xy_angle), 0),
                                                Point<dim>());
-  triangulation.set_boundary (0, boundary);
+  triangulation.set_manifold (0, boundary);
   triangulation.refine_global (2);
 
   for (typename Triangulation<dim>::active_cell_iterator
index 0e16552ed8d458d40c53038e21e0e4c79937ee93..5f1eca12c36a1770262dc8b66111ff809dcdf0d7 100644 (file)
@@ -30,7 +30,7 @@
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_reordering.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
@@ -86,8 +86,8 @@ void check ()
   // set bottom face to use MyBoundary
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (coarse_grid.begin_active()->face(f)->center()[dim-1] == -1)
-      coarse_grid.begin_active()->face(f)->set_boundary_indicator (1);
-  coarse_grid.set_boundary (1, my_boundary);
+      coarse_grid.begin_active()->face(f)->set_manifold_indicator (1);
+  coarse_grid.set_manifold (1, my_boundary);
 
   // now try to refine this one
   // cell. we should get an exception
index 715cbd57596f475a18222c91f0ef996b9123d799..b6c693af6a85c1d82526bf39d6a3ea50bc3d4282 100644 (file)
@@ -27,7 +27,7 @@
 #include <deal.II/grid/grid_reordering.h>
 #include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
@@ -119,8 +119,8 @@ void check ()
   // set bottom face to use MyBoundary
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (coarse_grid.begin_active()->face(f)->center()[dim-1] == -1)
-      coarse_grid.begin_active()->face(f)->set_boundary_indicator (1);
-  coarse_grid.set_boundary (1, my_boundary);
+      coarse_grid.begin_active()->face(f)->set_manifold_indicator (1);
+  coarse_grid.set_manifold (1, my_boundary);
 
   // now try to refine this one
   // cell. we should get an exception
index 7fe258c2454d6b6a2e45e6dcee2a1ec2f040ec55..a719d8713b2e14ec65b4ff6bc90c43fa29061bdb 100644 (file)
@@ -121,8 +121,8 @@ void check ()
   // set bottom face to use MyBoundary
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (coarse_grid.begin_active()->face(f)->center()[dim-1] == -1)
-      coarse_grid.begin_active()->face(f)->set_boundary_indicator (1);
-  coarse_grid.set_boundary (1, my_boundary);
+      coarse_grid.begin_active()->face(f)->set_manifold_indicator (1);
+  coarse_grid.set_manifold (1, my_boundary);
 
   // now try to refine this one
   // cell. we should get an exception
index 7e8bae09bf46c59e0fca222cab2145aa3fc19205..8c433a0e499983f7f17fc31afc86490c38ef47c7 100644 (file)
@@ -28,7 +28,7 @@
 #include <deal.II/grid/grid_reordering.h>
 #include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
@@ -79,8 +79,8 @@ void check ()
       if (coarse_grid.begin_active()->face(f)->line(e)->center()[0] == 0)
         if (coarse_grid.begin_active()->face(f)->line(e)->center()[1] == 0.5)
           if (coarse_grid.begin_active()->face(f)->line(e)->center()[2] == 0)
-            coarse_grid.begin_active()->face(f)->line(e)->set_boundary_indicator (99);
-  coarse_grid.set_boundary (99, my_boundary);
+            coarse_grid.begin_active()->face(f)->line(e)->set_manifold_indicator (99);
+  coarse_grid.set_manifold (99, my_boundary);
 
   // now try to refine this one
   // cell. we should get an exception
index 14c012824191bb5e6395434624ddf63ab6a83ed0..dafdce0272aaebe40c2945bc8f7d4e930f8d12e9 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/fe/fe_values.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 
 void test()
@@ -35,8 +35,8 @@ void test()
 
   // shell is so narrow that MappingQ(2) distorts the cell
   GridGenerator::quarter_hyper_shell (tria, Point<dim>(), 0.95, 1, 1);
-  static HyperShellBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
 
   FE_Nothing<dim> dummy;
   MappingQ<dim> mapping(2);
index cfbe9d6906ecf194c9d1062b73059b7194e4e317..4db00c3473a29187de67e8bbef2e4497d2657cdb 100644 (file)
@@ -27,7 +27,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/numerics/fe_field_function.h>
 #include <deal.II/numerics/vector_tools.h>
@@ -48,11 +48,11 @@ public:
 template<int dim>
 void test()
 {
-  const HyperBallBoundary<dim> boundary_description;
+  const PolarManifold<dim> boundary_description;
 
   Triangulation<dim>   triangulation;
   GridGenerator::hyper_ball (triangulation);
-  triangulation.set_boundary (0, boundary_description);
+  triangulation.set_manifold (0, boundary_description);
   triangulation.refine_global (1);
 
   FE_Q<dim> fe(2);
index 1d5876637be0bb1357697dbe74ab4374c72a078c..321151aa241a3a0372b7fc31badd530bc3c4c7d9 100644 (file)
@@ -27,7 +27,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/numerics/fe_field_function.h>
 #include <deal.II/numerics/vector_tools.h>
@@ -51,11 +51,11 @@ public:
 template<int dim>
 void test()
 {
-  const HyperBallBoundary<dim> boundary_description;
+  const PolarManifold<dim> boundary_description;
 
   Triangulation<dim>   triangulation;
   GridGenerator::hyper_ball (triangulation);
-  triangulation.set_boundary (0, boundary_description);
+  triangulation.set_manifold (0, boundary_description);
   triangulation.refine_global (1);
 
   FESystem<dim> fe(FE_Q<dim> (2),2);
index 080b7034360381a06bd8d8d8c1164473b284bd91..230d70d612dea3293b1875af965ab7fa81d0cc74 100644 (file)
@@ -31,7 +31,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/numerics/fe_field_function.h>
 #include <deal.II/numerics/vector_tools.h>
@@ -52,11 +52,11 @@ public:
 template<int dim>
 void test()
 {
-  const HyperBallBoundary<dim> boundary_description;
+  const PolarManifold<dim> boundary_description;
 
   Triangulation<dim>   triangulation;
   GridGenerator::hyper_ball (triangulation);
-  triangulation.set_boundary (0, boundary_description);
+  triangulation.set_manifold (0, boundary_description);
   triangulation.refine_global (1);
 
   FE_Q<dim> fe(2);
index 3bbc04d690140917019cc9f6a445d492c8081409..e81bcc21a5e66ae84560615f82be3725e2ab51f9 100644 (file)
@@ -31,7 +31,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/numerics/fe_field_function.h>
 #include <deal.II/numerics/vector_tools.h>
@@ -56,11 +56,11 @@ public:
 template<int dim>
 void test()
 {
-  const HyperBallBoundary<dim> boundary_description;
+  const PolarManifold<dim> boundary_description;
 
   Triangulation<dim>   triangulation;
   GridGenerator::hyper_ball (triangulation);
-  triangulation.set_boundary (0, boundary_description);
+  triangulation.set_manifold (0, boundary_description);
   triangulation.refine_global (1);
 
   FESystem<dim> fe(FE_Q<dim> (2),2);
index dba992f3058850c7d442e1a7efff3e25f2e127e6..5adcc3bcdf42e70d82c78a662c575af8c222d035 100644 (file)
@@ -31,7 +31,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/numerics/fe_field_function.h>
 #include <deal.II/numerics/vector_tools.h>
@@ -56,11 +56,11 @@ public:
 template<int dim>
 void test()
 {
-  const HyperBallBoundary<dim> boundary_description;
+  const PolarManifold<dim> boundary_description;
 
   Triangulation<dim>   triangulation;
   GridGenerator::hyper_ball (triangulation);
-  triangulation.set_boundary (0, boundary_description);
+  triangulation.set_manifold (0, boundary_description);
   triangulation.refine_global (1);
 
   FESystem<dim> fe(FE_Q<dim> (2),2);
index 1b77cd426ab04b68376f356882a0119b756bccb5..fce7cb69cd36cb842deaba28fc4614e82150e307 100644 (file)
@@ -31,7 +31,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/numerics/fe_field_function.h>
 #include <deal.II/numerics/vector_tools.h>
@@ -65,11 +65,11 @@ public:
 template<int dim>
 void test()
 {
-  const HyperBallBoundary<dim> boundary_description;
+  const PolarManifold<dim> boundary_description;
 
   Triangulation<dim>   triangulation;
   GridGenerator::hyper_ball (triangulation);
-  triangulation.set_boundary (0, boundary_description);
+  triangulation.set_manifold (0, boundary_description);
   triangulation.refine_global (1);
 
   FESystem<dim> fe(FE_Q<dim> (2),2);
index 39d88e7be3dcae292d5ed200d17d94ecaf0922e0..c29a84c8c42cab1911536d71f75e60a683ae20f7 100644 (file)
@@ -31,7 +31,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/numerics/fe_field_function.h>
 #include <deal.II/numerics/vector_tools.h>
@@ -65,11 +65,11 @@ public:
 template<int dim>
 void test()
 {
-  const HyperBallBoundary<dim> boundary_description;
+  const PolarManifold<dim> boundary_description;
 
   Triangulation<dim>   triangulation;
   GridGenerator::hyper_ball (triangulation);
-  triangulation.set_boundary (0, boundary_description);
+  triangulation.set_manifold (0, boundary_description);
   triangulation.refine_global (1);
 
   FESystem<dim> fe(FE_Q<dim> (4),2);
index 2348eff411980447b60f55c8864dcf25cdc1be38..60dfd89fdd3da8eada47cd70589dd7665e8a1861 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/fe/fe_values.h>
@@ -1044,8 +1044,8 @@ void LaplaceProblem<dim>::run ()
         {
           GridGenerator::hyper_ball (triangulation);
 
-          static const HyperBallBoundary<dim> boundary;
-          triangulation.set_boundary (0, boundary);
+          static const PolarManifold<dim> boundary;
+          triangulation.set_manifold (0, boundary);
 
           triangulation.refine_global (1);
         }
index c88304ad58ce411f5763a46433a34ae6b107459d..ff422ac458160b1f68e5e50aa1d435d75644f840 100644 (file)
@@ -32,7 +32,7 @@ std::ofstream logfile("output");
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/hp/dof_handler.h>
@@ -62,8 +62,8 @@ void compute_pi_by_area ()
   Triangulation<dim> triangulation;
   GridGenerator::hyper_ball (triangulation);
 
-  static const HyperBallBoundary<dim> boundary;
-  triangulation.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  triangulation.set_manifold (0, boundary);
 
   MappingQ<dim>     mapping(degree);
   const FE_Q<dim>   dummy_fe (1);
@@ -134,8 +134,8 @@ void compute_pi_by_perimeter ()
   Triangulation<dim> triangulation;
   GridGenerator::hyper_ball (triangulation);
 
-  static const HyperBallBoundary<dim> boundary;
-  triangulation.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  triangulation.set_manifold (0, boundary);
 
   const MappingQ<dim> mapping (degree);
   const FE_Q<dim>     fe (1);
index d46fbf3034c7f206f9ce7313f290d0834e518f6b..09a0b5dbce2e5ffe4fe46cc6a7f95d796b2f2c74 100644 (file)
@@ -30,7 +30,7 @@ std::ofstream logfile("output");
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/hp/dof_handler.h>
@@ -56,8 +56,8 @@ void gnuplot_output()
 
   Triangulation<dim> triangulation;
   GridGenerator::hyper_ball (triangulation);
-  static const HyperBallBoundary<dim> boundary;
-  triangulation.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  triangulation.set_manifold (0, boundary);
 
   for (unsigned int refinement=0; refinement<2;
        ++refinement, triangulation.refine_global(1))
@@ -99,8 +99,8 @@ void compute_pi_by_area ()
       Triangulation<dim> triangulation;
       GridGenerator::hyper_ball (triangulation);
 
-      static const HyperBallBoundary<dim> boundary;
-      triangulation.set_boundary (0, boundary);
+      static const PolarManifold<dim> boundary;
+      triangulation.set_manifold (0, boundary);
 
       MappingQ<dim> mapping(degree);
 
@@ -165,8 +165,8 @@ void compute_pi_by_perimeter ()
       Triangulation<dim> triangulation;
       GridGenerator::hyper_ball (triangulation);
 
-      static const HyperBallBoundary<dim> boundary;
-      triangulation.set_boundary (0, boundary);
+      static const PolarManifold<dim> boundary;
+      triangulation.set_manifold (0, boundary);
 
       const MappingQ<dim> mapping (degree);
       const FE_Q<dim>     fe (1);
index 28d97a1eab0e7feac711b6d2857b286e0bd612b3..a88b70875c91422a8ca1b3df0afcab61575b64f8 100644 (file)
@@ -35,7 +35,7 @@ std::ofstream logfile("output");
 #include <deal.II/lac/precondition.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/dofs/dof_handler.h>
@@ -197,8 +197,8 @@ template <int dim>
 void LaplaceProblem<dim>::run ()
 {
   GridGenerator::hyper_ball (triangulation);
-  static const HyperBallBoundary<dim> boundary;
-  triangulation.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  triangulation.set_manifold (0, boundary);
 
   for (unsigned int cycle=0; cycle<6; ++cycle, triangulation.refine_global(1))
     {
index 2a825122fde4ca70cd8a9f4a06ca734d89ed676b..e183d669c00296752d57d658f48ac1a0bd4db88c 100644 (file)
@@ -48,7 +48,7 @@ std::ofstream logfile("output");
 
 #include <deal.II/grid/grid_in.h>
 
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 #include <fstream>
 #include <sstream>
@@ -306,8 +306,8 @@ void LaplaceProblem<dim>::run ()
 
           grid_in.read_ucd (input_file);
 
-          static const HyperBallBoundary<dim> boundary;
-          triangulation.set_boundary (0, boundary);
+          static const PolarManifold<dim> boundary;
+          triangulation.set_manifold (0, boundary);
         }
 
       deallog << "   Number of active cells: "
index 8fa10830c465e8b9f42a43539cdedab0d82bbaa4..d61122cc32b5c8e1e5a22dd6983564444203e9b0 100644 (file)
@@ -39,7 +39,7 @@ std::ofstream logfile("output");
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/hp/fe_values.h>
@@ -338,8 +338,8 @@ void LaplaceProblem<dim>::run ()
         {
           GridGenerator::hyper_ball (triangulation);
 
-          static const HyperBallBoundary<dim> boundary;
-          triangulation.set_boundary (0, boundary);
+          static const PolarManifold<dim> boundary;
+          triangulation.set_manifold (0, boundary);
 
           triangulation.refine_global (1);
         }
index 5f9d260a38564c6550133d4026577282c3f4edb9..762f75f524f073478f1a8ccec292053855eeaa79 100644 (file)
@@ -32,7 +32,7 @@
 
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_tools.h>
 
@@ -59,10 +59,10 @@ int main ()
     const int dim = 3;
     deallog << "Testing hyper_cube in dim: " << dim << "..."<< endl;
 
-    const HyperBallBoundary<dim> boundary_description;
+    const PolarManifold<dim> boundary_description;
     Triangulation<dim> volume_mesh;
     GridGenerator::hyper_ball(volume_mesh);
-    volume_mesh.set_boundary (0, boundary_description);
+    volume_mesh.set_manifold (0, boundary_description);
 
     // exclude one of the 6 faces
     // from the surface mesh
@@ -75,7 +75,7 @@ int main ()
         for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
           if (cell->at_boundary(f))
             {
-              cell->face(f)->set_boundary_indicator(1);
+              cell->face(f)->set_manifold_indicator(1);
               done = true;
               break;
             }
@@ -85,7 +85,7 @@ int main ()
 
     const HyperBallBoundary<dim-1,dim> surface_description;
     Triangulation<dim-1,dim> boundary_mesh;
-    boundary_mesh.set_boundary (0, surface_description);
+    boundary_mesh.set_manifold (0, surface_description);
 
     // now extract a mesh of the 5
     // surface faces
@@ -108,7 +108,7 @@ int main ()
          cell != boundary_mesh.end(); ++cell)
       for (unsigned int f=0; f<GeometryInfo<dim-1>::faces_per_cell; ++f)
         if (cell->at_boundary(f))
-          cell->face(f)->set_boundary_indicator(1);
+          cell->face(f)->set_manifold_indicator(1);
 
     boundary_mesh.refine_global (2);
 
index 573a45066d963b384260874568f3c94e7ee8718d..66d32023d7181b57f38f3596e766c46e2ea5e9fe 100644 (file)
@@ -31,7 +31,7 @@
 
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_tools.h>
 
@@ -96,15 +96,15 @@ int main ()
     map< Triangulation<dim-1,dim>::cell_iterator,
          Triangulation<dim,dim>::face_iterator>
          surface_to_volume_mapping;
-    const HyperBallBoundary<dim> boundary_description;
+    const PolarManifold<dim> boundary_description;
     Triangulation<dim> volume_mesh;
     GridGenerator::hyper_ball(volume_mesh);
-    volume_mesh.set_boundary (0, boundary_description);
+    volume_mesh.set_manifold (0, boundary_description);
     volume_mesh.refine_global (1);
 
     const HyperBallBoundary<dim-1,dim> surface_description;
     Triangulation<dim-1,dim> boundary_mesh;
-    boundary_mesh.set_boundary (0, surface_description);
+    boundary_mesh.set_manifold (0, surface_description);
 
     surface_to_volume_mapping
       = GridTools::extract_boundary_mesh (volume_mesh, boundary_mesh);
index b820cb439162a8734e23750515f35919012860a8..db3ac042456349736fc75a50d8aeecea4cd16348 100644 (file)
@@ -29,7 +29,7 @@
 
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_tools.h>
 
@@ -104,7 +104,7 @@ int main ()
     map< Triangulation<dim-1,dim>::cell_iterator,
          Triangulation<dim,dim>::face_iterator>
          surface_to_volume_mapping;
-    const HyperBallBoundary<dim> boundary_description;
+    const PolarManifold<dim> boundary_description;
     Triangulation<dim> volume_mesh;
     GridGenerator::hyper_ball(volume_mesh);
     for (Triangulation<dim>::active_cell_iterator
@@ -113,12 +113,12 @@ int main ()
       for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
         if (cell->at_boundary(f))
           cell->face(f)->set_all_boundary_indicators (1);
-    volume_mesh.set_boundary (1, boundary_description);
+    volume_mesh.set_manifold (1, boundary_description);
     volume_mesh.refine_global (1);
 
     const HyperBallBoundary<dim-1,dim> surface_description;
     Triangulation<dim-1,dim> boundary_mesh;
-    boundary_mesh.set_boundary (1, surface_description);
+    boundary_mesh.set_manifold (1, surface_description);
 
     surface_to_volume_mapping
       = GridTools::extract_boundary_mesh (volume_mesh, boundary_mesh);
index 942651a8081818a90211c803546567dc69ee2b40..39e02bccda6642793a56a8c7c13a265d168be3f0 100644 (file)
@@ -29,7 +29,7 @@
 
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_tools.h>
 
@@ -104,7 +104,7 @@ int main ()
     map< Triangulation<dim-1,dim>::cell_iterator,
          Triangulation<dim,dim>::face_iterator>
          surface_to_volume_mapping;
-    const HyperBallBoundary<dim> boundary_description;
+    const PolarManifold<dim> boundary_description;
     Triangulation<dim> volume_mesh;
     GridGenerator::hyper_ball(volume_mesh);
     for (Triangulation<dim>::active_cell_iterator
@@ -114,14 +114,14 @@ int main ()
         if (cell->at_boundary(f) && (cell->center()[dim-1]>0))
           cell->face(f)->set_all_boundary_indicators (1);
 
-    volume_mesh.set_boundary (1, boundary_description);
-    volume_mesh.set_boundary (0, boundary_description);
+    volume_mesh.set_manifold (1, boundary_description);
+    volume_mesh.set_manifold (0, boundary_description);
     volume_mesh.refine_global (1);
 
     const HyperBallBoundary<dim-1,dim> surface_description;
     Triangulation<dim-1,dim> boundary_mesh;
-    boundary_mesh.set_boundary (1, surface_description);
-    boundary_mesh.set_boundary (0, surface_description);
+    boundary_mesh.set_manifold (1, surface_description);
+    boundary_mesh.set_manifold (0, surface_description);
 
     set<types::boundary_id> boundary_ids;
     boundary_ids.insert(1);
index 9815e309650ca0981390bb6126f5afc1cfc64213..f74f631e62a7d13f80506d694c75e19bb6d1a52c 100644 (file)
@@ -27,7 +27,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_in.h>
 #include <deal.II/grid/grid_out.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 #include <fstream>
 #include <string>
@@ -39,7 +39,7 @@ void test(std::string filename)
 {
   HyperBallBoundary<dim, spacedim> boundary;
   Triangulation<dim, spacedim> tria;
-  tria.set_boundary(1, boundary);
+  tria.set_manifold(1, boundary);
   GridIn<dim, spacedim> gi;
   gi.attach_triangulation (tria);
   std::ifstream in (filename.c_str());
index fd527e3d1b6dc6ad2521371a6fada9315af101b4..b26fa257b6e973364df3f08203b4f597e48db1e4 100644 (file)
@@ -28,7 +28,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_tools.h>
 #include <deal.II/fe/fe_values.h>
@@ -44,9 +44,9 @@ void test ()
   Triangulation<dim> volume_mesh;
   GridGenerator::hyper_ball(volume_mesh);
 
-  const HyperBallBoundary<dim-1,dim> surface_description;
+  const PolarManifold<dim> surface_description;
   Triangulation<dim-1,dim> boundary_mesh;
-  boundary_mesh.set_boundary (0, surface_description);
+  boundary_mesh.set_manifold (0, surface_description);
 
   GridTools::extract_boundary_mesh (volume_mesh, boundary_mesh);
 
index 0849d7b8af78e56bf48b0630507991889d0503b0..cbe2afdf201861f092bbe70c2d32fefeccc91b23 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_tools.h>
 #include <deal.II/fe/fe_values.h>
@@ -43,7 +43,7 @@ void test ()
 
   const HyperBallBoundary<dim-1,dim> surface_description;
   Triangulation<dim-1,dim> boundary_mesh;
-  boundary_mesh.set_boundary (0, surface_description);
+  boundary_mesh.set_manifold (0, surface_description);
 
   GridTools::extract_boundary_mesh (volume_mesh, boundary_mesh);
 
index af63da020b28edd98000293075c80abbd3d9209f..1c4b86e7623b539b94c9a553874f2487d8f1389a 100644 (file)
@@ -27,7 +27,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 #include <fstream>
 #include <string>
@@ -45,7 +45,7 @@ int main ()
 
   TorusBoundary<dim, spacedim> boundary (1.5, .5);
   Triangulation<dim, spacedim> tria;
-  tria.set_boundary (0, boundary);
+  tria.set_manifold (0, boundary);
 
   GridGenerator::torus (tria, 1.5, .5);
   tria.refine_global(2);
index 60372d4f29eccef7f8f29a27c99c68df1a593e6a..a9188fa4fcfdd38d0f43952e2b67f97f362e37a8 100644 (file)
@@ -28,7 +28,7 @@
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/dofs/dof_handler.h>
@@ -99,10 +99,10 @@ check ()
     }
   else
     GridGenerator::hyper_cube(tr, -1./std::sqrt(static_cast<double>(dim)),1./std::sqrt(static_cast<double>(dim)));
-  static const HyperBallBoundary<dim> boundary;
+  static const PolarManifold<dim> boundary;
   if (dim != 1)
     {
-      tr.set_boundary (0, boundary);
+      tr.set_manifold (0, boundary);
     }
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
index bf41611661b081d593557798bb7b4afc17d349be..6b3c8d2a69332b87ca8aa604bfc8e35ee049c1c6 100644 (file)
@@ -21,7 +21,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/fe/fe_q.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/lac/sparse_matrix.h>
@@ -46,13 +46,13 @@ std::ofstream logfile("output");
 template <int dim>
 class Ball
   :
-  public StraightBoundary<dim>
+  public FlatManifold<dim>
 {
 public:
   virtual Point<dim>
   get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
   {
-    Point<dim> middle = StraightBoundary<dim>::get_new_point_on_line(line);
+    Point<dim> middle = FlatManifold<dim>::get_new_point_on_line(line);
 
     for (int i=0; i<dim; ++i)
       middle(i) -= .5;
@@ -67,7 +67,7 @@ public:
   virtual Point<dim>
   get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
   {
-    Point<dim> middle = StraightBoundary<dim>::get_new_point_on_quad(quad);
+    Point<dim> middle = FlatManifold<dim>::get_new_point_on_quad(quad);
 
     for (int i=0; i<dim; ++i)
       middle(i) -= .5;
@@ -83,7 +83,7 @@ public:
 template <int dim>
 class CurvedLine
   :
-  public StraightBoundary<dim>
+  public FlatManifold<dim>
 {
 public:
   virtual Point<dim>
@@ -98,7 +98,7 @@ template <int dim>
 Point<dim>
 CurvedLine<dim>::get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
 {
-  Point<dim> middle = StraightBoundary<dim>::get_new_point_on_line (line);
+  Point<dim> middle = FlatManifold<dim>::get_new_point_on_line (line);
 
   // if the line is at the top of bottom
   // face: do a special treatment on
@@ -147,7 +147,7 @@ template <int dim>
 Point<dim>
 CurvedLine<dim>::get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
 {
-  Point<dim> middle = StraightBoundary<dim>::get_new_point_on_quad (quad);
+  Point<dim> middle = FlatManifold<dim>::get_new_point_on_quad (quad);
 
   // if the face is at the top of bottom
   // face: do not move the midpoint in
@@ -266,12 +266,12 @@ void TestCases<dim>::run (const unsigned int test_case)
     {
       if (dim==3)
         {
-          tria->begin_active()->face(4)->set_boundary_indicator(1);
-          tria->begin_active()->face(5)->set_boundary_indicator(1);
+          tria->begin_active()->face(4)->set_manifold_indicator(1);
+          tria->begin_active()->face(5)->set_manifold_indicator(1);
         };
 
       // set the boundary function
-      tria->set_boundary(1, (test_case==2)
+      tria->set_manifold(1, (test_case==2)
                          ? ((Boundary<dim> &)ball) : ((Boundary<dim> &)curved_line));
 
       // refine once
index e0f3e1d58a4a7492d829f49b2dc43a80ae700b42..b669f007bd8aa474f5397c41435cf946a6200a9b 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -95,8 +95,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index 32c32e9a233dadbde0a40e0ae164ae24e71e3390..d1eaf982482f796cbf9900a6f8496eb2e5a58368 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -120,8 +120,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index f0e57a117cc4cf0d2b908f6a8c336ffdd9422b31..8d1d98949b624490073debde792a9eb73ebddb19 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_raviart_thomas.h>
@@ -101,8 +101,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index 74eee718072d09ef78d9ace665bd95eab5a5468f..66634f246f9d8e92992b4f8c5da09c5db6d00bf4 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -121,8 +121,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index 97fc058ebd9116e37a2d39671bc89c455b6b6da4..3280b8e23b93f405135d4080d4b989f11dd9cb43 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -87,8 +87,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index a4fb1c78ce7d65857ad72afd054dfcf7866c9d59..6c2e7ebacf29e62ba3e37d2fcd39a441e08f70f0 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -88,8 +88,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index d0ce6a04efb6bc37edfc9b8a68e880b141bdb0d4..c2f2ef5019e1f7a691dce91e0a88e7c589ba602a 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -88,8 +88,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index 9215f4c530afe7bfd1f2adb911e97c3fac401982..75c9e11a4c0e5644184c75d770a1ad4663e364c9 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -89,8 +89,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index 2f2cff48c0f32db6c7b5a2033242ff75baa3bac9..32598a876194148ec5ddaef2144bddfd1b333f6a 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -92,8 +92,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index b1436f6df953645d3dc59946e5e47f968b35a0a1..357f764086340129a005a22b2f7a4e78f39bc8ca 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -93,8 +93,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index 4ac711fffa0ffb2dc911981359a42982b8c313e9..ffd8f400711f76e1363ebe703740703fb91c410a 100644 (file)
@@ -28,7 +28,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -97,8 +97,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FE_RaviartThomas<dim> fe(1);
   test(tr, fe);
index 20ecbbe0738c48675b6446aa5a2de0768d307d65..253cda5821f8b5a4156d924f0961e3c62b543458 100644 (file)
@@ -28,7 +28,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -97,8 +97,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FE_Nedelec<dim> fe(1);
   test(tr, fe);
index 0d0cb82d17e7a77706c19e24e68db9c9ebcafe7e..4806ae1e991c6eff651728d338ccb90a292a24c1 100644 (file)
@@ -28,7 +28,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -97,8 +97,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe(FE_RaviartThomas<dim>(1),1);
   test(tr, fe);
index 923ab6b65ac5a842913fe68a4cff2a80cf09ddfe..c26ee8e24cddd3a1ffcfa92ee4f9260fef541458 100644 (file)
@@ -30,7 +30,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -95,8 +95,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index 42328f17e872f431aeb9505aa15c45eb6c36d892..eaa03e4eac9e919473fb601eab62b7508d8ff79c 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -93,8 +93,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index 3a26c309e540c1301285027b7ce6856b1b1e0c8b..42fc0cc381b1657c06afef76493e7aefc4409386 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -94,8 +94,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index e357f019042fa52c019770c80205971499a6b53f..f0ce11914c8b6cdc699f6d95804a8f317cd9c565 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -94,8 +94,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index 834662625389a2bec37d3827e4b864f1f1854689..d543530d1f733368844dccd35147ecc0ff03674a 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -97,8 +97,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index 7104b9ef3519ff92f8da4b44b3b1782757fc7cc4..4aa6b44bee851cf2fea866f1aca21ede3713c1a9 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -98,8 +98,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index e906d800df9e1b49cd21c72215c405ce14936f2e..099e49abfd8fc2fb6b7e638cd0643e7f7c72f011 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -99,8 +99,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index 498fa3dae79d14aacbb9fc78ff98996172f96b32..beb2352df7f835a2bb2c626882d8e8cf588a80fa 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -98,8 +98,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index 59dec4c5d354a6303f07b6237fcc91fa96fbbc28..0fcc9bf268fc5bf23e9acdc37371888671d95253 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -103,8 +103,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index 25e40dac85e0acfcda4b6bd14eaf519e4d65ae4e..a22f30773c7f387a1ede25160f0f49996bc8cb74 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -98,8 +98,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_Q<dim>(2), 2,
index 06efab414689c8d3f4fb7c8420115de31817915a..8a954b58a820cc9ab322e17c0c99bc58e513790e 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -102,8 +102,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1), 1,
                     FE_RaviartThomas<dim>(1), 1,
index 3f99ca5cebcbfd88000ca7e5ac624014f1dd866a..071d9cd11fb3ee6ab0b4421604f5aed0b39df883 100644 (file)
@@ -26,7 +26,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_system.h>
@@ -91,8 +91,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1),
                     SymmetricTensor<2,dim>::n_independent_components);
index 2e17287429f8c48604d5e08db14f11001ffd2526..775dec95b655299c9d067f3eea121a7b9f3f045c 100644 (file)
@@ -26,7 +26,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_system.h>
@@ -91,8 +91,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1),
                     Tensor<2,dim>::n_independent_components);
index dbf9c01df3f6dc7f9e178370bfdd6150b8e6b6fa..1a8e9c45b48d100c13585dde80e6dab1af278fb6 100644 (file)
@@ -26,7 +26,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_system.h>
@@ -85,8 +85,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1),
                     SymmetricTensor<2,dim>::n_independent_components);
index fc0f020cb33289ed2b6dc197c666204c47c4d17c..c37d411f5fe871d47272179d7b7d3533b96ab5fe 100644 (file)
@@ -26,7 +26,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_system.h>
@@ -85,8 +85,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   FESystem<dim> fe (FE_Q<dim>(1),
                     Tensor<2,dim>::n_independent_components);
index 0b196be527e6f25146062f386bec730f6b0c9eab..5d8054282ce296c8b3137b82e483082a4680d3d0 100644 (file)
@@ -20,7 +20,7 @@
 // refinement
 
 #include "../tests.h"
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria.h>
@@ -46,8 +46,8 @@ void check (double r1, double r2, unsigned int n)
   Point<dim> center;
   Triangulation<dim> tria (Triangulation<dim>::none, true);
   GridGenerator::hyper_shell (tria, center, r1, r2, n, true);
-  static const HyperShellBoundary<dim> boundary(center);
-  tria.set_boundary(0, boundary);
+  static const PolarManifold<dim> boundary(center);
+  tria.set_manifold(0, boundary);
 
   for (unsigned int i=0; i<2; ++i)
     {
index 9ae80cc5438a7ab32a856950e1eebc4ea123de88..322b3452c1d529d41a8893f9d203cd65754024c4 100644 (file)
@@ -19,7 +19,7 @@
 // this series test the other possibilities for the number of cells
 
 #include "../tests.h"
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria.h>
@@ -43,8 +43,8 @@ void check (double r1, double r2, unsigned int n)
   Point<dim> center;
   Triangulation<dim> tria (Triangulation<dim>::none);
   GridGenerator::hyper_shell (tria, center, r1, r2, n);
-  static const HyperShellBoundary<dim> boundary(center);
-  tria.set_boundary(0, boundary);
+  static const PolarManifold<dim> boundary(center);
+  tria.set_manifold(0, boundary);
   tria.refine_global(1);
 
   GridOut grid_out;
index ae1ac3d012479ea411d219ed276769ee8b10b764..27561c9449a565f9d449227d6392b327ecc6244b 100644 (file)
@@ -19,7 +19,7 @@
 // this series test the other possibilities for the number of cells
 
 #include "../tests.h"
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria.h>
@@ -43,8 +43,8 @@ void check (double r1, double r2, unsigned int n)
   Point<dim> center;
   Triangulation<dim> tria (Triangulation<dim>::none);
   GridGenerator::hyper_shell (tria, center, r1, r2, n);
-  static const HyperShellBoundary<dim> boundary(center);
-  tria.set_boundary(0, boundary);
+  static const PolarManifold<dim> boundary(center);
+  tria.set_manifold(0, boundary);
   tria.refine_global(1);
 
   GridOut grid_out;
index e51778c7283e4888e633eeae1a8c115ead6f420d..7c802f2af64a4da01b74d64d9aeb07afa881ccd4 100644 (file)
@@ -19,7 +19,7 @@
 // this series test the other possibilities for the number of cells
 
 #include "../tests.h"
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria.h>
@@ -43,8 +43,8 @@ void check (double r1, double r2, unsigned int n)
   Point<dim> center;
   Triangulation<dim> tria (Triangulation<dim>::none);
   GridGenerator::hyper_shell (tria, center, r1, r2, n);
-  static const HyperShellBoundary<dim> boundary(center);
-  tria.set_boundary(0, boundary);
+  static const PolarManifold<dim> boundary(center);
+  tria.set_manifold(0, boundary);
 
   GridOut grid_out;
   grid_out.write_gnuplot (tria, logfile);
index be84422aae26bb59c4019467a6167541aed05582..d2a65a538c9281e1e94115905c414d1f9dd24fb9 100644 (file)
@@ -19,8 +19,8 @@
 #include "../tests.h"
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/base/logstream.h>
@@ -36,10 +36,10 @@ template <int dim>
 void test ()
 {
   Triangulation<dim> tria;
-  static const HyperBallBoundary<dim> x;
+  static const PolarManifold<dim> x;
   if (dim == 2)
     {
-      tria.set_boundary (0, x);
+      tria.set_manifold (0, x);
       GridGenerator::hyper_ball (tria);
     }
   else
index 7ec3897f182764c96aab47c9e5245809bdda8ce4..d8e2ee95c66844da65669646c8063aab211e1487 100644 (file)
@@ -23,7 +23,7 @@
 #include <deal.II/base/logstream.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_in.h>
 
@@ -40,8 +40,8 @@ void test(std::ostream &out)
   Triangulation<dim> tr;
 
   GridGenerator::hyper_cube_with_cylindrical_hole(tr, .3, .4, 1, 1, false);
-  CylinderBoundary<dim> boundary(.3, 2);
-  tr.set_boundary(1, boundary);
+  CylinderManifold<dim> boundary(.3, 2);
+  tr.set_manifold(1, boundary);
   {
     std::ofstream grid_file("coarse_grid.inp");
     go.write_ucd(tr, grid_file);
@@ -51,7 +51,7 @@ void test(std::ostream &out)
   deallog << "Writing refined from constructor" << std::endl;
   go.write_ucd(tr, out);
 
-  tr.set_boundary(1);
+  tr.set_manifold(1);
   tr.clear();
 
   GridIn<dim> gi;
@@ -63,7 +63,7 @@ void test(std::ostream &out)
     grid_file.close();
   }
 
-  tr.set_boundary(1, boundary);
+  tr.set_manifold(1, boundary);
   tr.refine_global(1);
   deallog << "Writing refined from file" << std::endl;
   go.write_ucd(tr, out);
@@ -74,7 +74,7 @@ void test(std::ostream &out)
     grid_file.close();
   }
 
-  tr.set_boundary(1);
+  tr.set_manifold(1);
   tr.clear();
 }
 
index 1b333ea2500c32971e2c584355b9f09871ce4a2a..f6b447073dbbf506a47aeeee4013124c41a16a15 100644 (file)
@@ -17,7 +17,7 @@
 
 
 #include "../tests.h"
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria.h>
@@ -40,13 +40,13 @@ std::ofstream logfile("output");
 
 template <int dim>
 class Ball :
-  public StraightBoundary<dim>
+  public FlatManifold<dim>
 {
 public:
   virtual Point<dim>
   get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
   {
-    Point<dim> middle = StraightBoundary<dim>::get_new_point_on_line(line);
+    Point<dim> middle = FlatManifold<dim>::get_new_point_on_line(line);
 
     for (int i=0; i<dim; ++i)
       middle(i) -= .5;
@@ -61,7 +61,7 @@ public:
   virtual Point<dim>
   get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
   {
-    Point<dim> middle = StraightBoundary<dim>::get_new_point_on_quad(quad);
+    Point<dim> middle = FlatManifold<dim>::get_new_point_on_quad(quad);
 
     for (int i=0; i<dim; ++i)
       middle(i) -= .5;
@@ -76,7 +76,7 @@ public:
 
 template <int dim>
 class CurvedLine :
-  public StraightBoundary<dim>
+  public FlatManifold<dim>
 {
 public:
   virtual Point<dim>
@@ -91,7 +91,7 @@ template <int dim>
 Point<dim>
 CurvedLine<dim>::get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
 {
-  Point<dim> middle = StraightBoundary<dim>::get_new_point_on_line (line);
+  Point<dim> middle = FlatManifold<dim>::get_new_point_on_line (line);
 
   // if the line is at the top of bottom
   // face: do a special treatment on
@@ -139,7 +139,7 @@ template <int dim>
 Point<dim>
 CurvedLine<dim>::get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
 {
-  Point<dim> middle = StraightBoundary<dim>::get_new_point_on_quad (quad);
+  Point<dim> middle = FlatManifold<dim>::get_new_point_on_quad (quad);
 
   // if the face is at the top of bottom
   // face: do not move the midpoint in
@@ -222,8 +222,8 @@ void test (const int test_case)
     {
       if (dim==3)
         {
-          tria.begin_active()->face(4)->set_boundary_indicator(1);
-          tria.begin_active()->face(5)->set_boundary_indicator(1);
+          tria.begin_active()->face(4)->set_manifold_indicator(1);
+          tria.begin_active()->face(5)->set_manifold_indicator(1);
         };
 
 
@@ -231,9 +231,9 @@ void test (const int test_case)
       Ball<dim>       ball;
       CurvedLine<dim> curved_line;
       if (test_case==2)
-        tria.set_boundary (1, ball);
+        tria.set_manifold (1, ball);
       else
-        tria.set_boundary (1, curved_line);
+        tria.set_manifold (1, curved_line);
 
       // refine once
       tria.begin_active()->set_refine_flag();
@@ -255,7 +255,7 @@ void test (const int test_case)
           tria.execute_coarsening_and_refinement();
         };
 
-      tria.set_boundary (1);
+      tria.set_manifold (1);
       break;
     }
     }
index b342d00ef71bd7e36a8c6c4035d4f30933502789..795bf7d07b69aca2cc06cc6e1396349d7a0c21fc 100644 (file)
@@ -21,7 +21,7 @@
 #include "../tests.h"
 #include <deal.II/base/logstream.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/grid_out.h>
@@ -60,8 +60,8 @@ void test1 ()
       Triangulation<dim> tria;
       GridGenerator::hyper_ball(tria, Point<dim>(), 1);
 
-      static const HyperBallBoundary<dim> boundary;
-      tria.set_boundary (0, boundary);
+      static const PolarManifold<dim> boundary;
+      tria.set_manifold (0, boundary);
 
       for (unsigned int i=0; i<4; ++i)
         {
index b417626514e733cd8d22e39095cce5f71b437cd8..cf76d1edc7d19438ea81cad6de4b843dfe3c7465 100644 (file)
@@ -21,7 +21,7 @@
 #include <deal.II/base/logstream.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/tria_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/fe/mapping_q.h>
@@ -34,10 +34,10 @@ int main ()
 {
   const unsigned int dim=2;
   Point<dim> origin;
-  HyperShellBoundary<dim> boundary(origin);
+  PolarManifold<dim> boundary(origin);
   MappingQ<dim> mapping(2);
   Triangulation<dim> tria;
-  tria.set_boundary(0, boundary);
+  tria.set_manifold(0, boundary);
   const double inner_radius=1.;
   const double outer_radius=5.;
   GridGenerator::hyper_shell(tria, origin, inner_radius, outer_radius, 8);
@@ -88,7 +88,7 @@ int main ()
                       // circle.
                       new_points.insert(std::pair<types::global_dof_index, Point<dim> > (
                                           face->vertex_index(vertex_no), n_radius/inner_radius*v+n_center));
-                      face->set_boundary_indicator(1);
+                      face->set_manifold_indicator(1);
                     }
                   else
                     Assert(false, ExcInternalError());
@@ -97,8 +97,8 @@ int main ()
     }
 
   GridGenerator::laplace_transformation (tria, new_points);
-  HyperBallBoundary<dim> inner_ball(n_center, n_radius);
-  tria.set_boundary(1, inner_ball);
+  PolarManifold<dim> inner_ball(n_center, n_radius);
+  tria.set_manifold(1, inner_ball);
 
   GridOut grid_out;
   std::ofstream eps_stream2("output");
index 72f86fd92887f1d6c45a5ed3546555893248834d..3a312a240f97e339aaf0f63050a8445867b34107 100644 (file)
@@ -23,7 +23,7 @@
 #include <deal.II/base/logstream.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/tria_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
 
@@ -35,9 +35,9 @@ int main ()
 {
   const unsigned int dim=3;
   Point<dim> origin;
-  CylinderBoundary<dim> boundary;
+  CylinderManifold<dim> boundary;
   Triangulation<dim> tria;
-  tria.set_boundary(0, boundary);
+  tria.set_manifold(0, boundary);
   GridGenerator::cylinder(tria, 1, .7);
   tria.refine_global(2);
 
index 804e2d24e4a0577691d584447ae9ba535ed8e221..49e8dde817724ad6f8999254165eac785857fa21 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/fe/fe_values.h>
@@ -319,8 +319,8 @@ void Step6<dim>::run ()
             {
               GridGenerator::hyper_ball (triangulation);
 
-              static const HyperBallBoundary<dim> boundary;
-              triangulation.set_boundary (0, boundary);
+              static const PolarManifold<dim> boundary;
+              triangulation.set_manifold (0, boundary);
 
               triangulation.refine_global (1);
             }
index b1aa01fbd4b5a08e5616583a0c1c69a067c344e8..b0f590dcbc276a899ccfb75f6a38bf2d67e8785d 100644 (file)
@@ -28,7 +28,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/lac/constraint_matrix.h>
 #include <deal.II/fe/fe_q.h>
@@ -66,8 +66,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   tr.refine_global(2);
 
index 19a06cd992df293317313fceedba45f1a5a31d4c..8f8c9a634a7dd372a9090244d2da2cc6466d8628 100644 (file)
@@ -29,7 +29,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/lac/constraint_matrix.h>
 #include <deal.II/fe/fe_q.h>
@@ -125,8 +125,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   tr.refine_global(2);
 
index 6050fc89af07509db24662271027ae1bda0314b3..5325aee1bd3bfff06da2efb7ed66a317c39193d1 100644 (file)
@@ -32,7 +32,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/dofs/dof_renumbering.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -79,8 +79,8 @@ void test_hyper_cube()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   tr.refine_global(1);
 
index 347e2f79db483460538cf918439af87e504c6914..5e93d8eb94521ddeee4496fcc073936c71cc5fd7 100644 (file)
@@ -34,7 +34,7 @@
 
 #include <grid/tria.h>
 #include <grid/grid_generator.h>
-#include <grid/tria_boundary_lib.h>
+#include <grid/manifold_lib.h>
 #include <grid/grid_out.h>
 #include <dofs/dof_handler.h>
 #include <fe/fe_system.h>
@@ -68,51 +68,51 @@ colorize_sixty_deg_hyper_shell(Triangulation<3> &tria,
         double radius = cell->face(f)->center().norm() - center.norm();
         if (std::fabs(cell->face(f)->center()(2) - sqrt(3.)*cell->face(f)->center()(0)) < eps ) // z = sqrt(3)x set boundary 2
           {
-            cell->face(f)->set_boundary_indicator(2);
+            cell->face(f)->set_manifold_indicator(2);
             for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
               if (cell->face(f)->line(j)->at_boundary())
                 if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) > eps)
-                  cell->face(f)->line(j)->set_boundary_indicator(2);
+                  cell->face(f)->line(j)->set_manifold_indicator(2);
           }
         else if (std::fabs(cell->face(f)->center()(2) + sqrt(3.)*cell->face(f)->center()(0)) < eps) // z = -sqrt(3)x set boundary 3
           {
-            cell->face(f)->set_boundary_indicator(3);
+            cell->face(f)->set_manifold_indicator(3);
             for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
               if (cell->face(f)->line(j)->at_boundary())
                 if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) > eps)
-                  cell->face(f)->line(j)->set_boundary_indicator(3);
+                  cell->face(f)->line(j)->set_manifold_indicator(3);
           }
         else if (std::fabs(cell->face(f)->center()(2) - sqrt(3.)*cell->face(f)->center()(1)) < eps ) // z = sqrt(3)y set boundary 4
           {
-            cell->face(f)->set_boundary_indicator(4);
+            cell->face(f)->set_manifold_indicator(4);
             for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
               if (cell->face(f)->line(j)->at_boundary())
                 if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) > eps)
-                  cell->face(f)->line(j)->set_boundary_indicator(4);
+                  cell->face(f)->line(j)->set_manifold_indicator(4);
           }
         else if (std::fabs(cell->face(f)->center()(2) + sqrt(3.)*cell->face(f)->center()(1)) < eps ) // z = -sqrt(3)y set boundary 5
           {
-            cell->face(f)->set_boundary_indicator(5);
+            cell->face(f)->set_manifold_indicator(5);
             for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
               if (cell->face(f)->line(j)->at_boundary())
                 if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) > eps)
-                  cell->face(f)->line(j)->set_boundary_indicator(5);
+                  cell->face(f)->line(j)->set_manifold_indicator(5);
           }
         else if (radius < middle) // inner radius set boundary 0
           {
-            cell->face(f)->set_boundary_indicator(0);
+            cell->face(f)->set_manifold_indicator(0);
             for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
               if (cell->face(f)->line(j)->at_boundary())
                 if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) < eps)
-                  cell->face(f)->line(j)->set_boundary_indicator(0);
+                  cell->face(f)->line(j)->set_manifold_indicator(0);
           }
         else if (radius > middle) // outer radius set boundary 1
           {
-            cell->face(f)->set_boundary_indicator(1);
+            cell->face(f)->set_manifold_indicator(1);
             for (unsigned int j=0; j<GeometryInfo<3>::lines_per_face; ++j)
               if (cell->face(f)->line(j)->at_boundary())
                 if (std::fabs(cell->face(f)->line(j)->vertex(0).norm() - cell->face(f)->line(j)->vertex(1).norm()) < eps)
-                  cell->face(f)->line(j)->set_boundary_indicator(1);
+                  cell->face(f)->line(j)->set_manifold_indicator(1);
           }
         else
           AssertThrow (false, ExcInternalError());
@@ -171,8 +171,8 @@ void run()
                          0.5,
                          1.0);
 
-  static HyperShellBoundary<dim> boundary((Point<dim>()));
-  triangulation.set_boundary (0, boundary);
+  static PolarManifold<dim> boundary((Point<dim>()));
+  triangulation.set_manifold (0, boundary);
 
   // write out the mesh. may not be
   // strictly necessary here, but is
index f304cf93876bb5da99f01b1852433e875cd61b50..a264c80bc19a60a1a5444ca81ff2a5c522603f83 100644 (file)
@@ -21,7 +21,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/hp/dof_handler.h>
 #include <deal.II/hp/fe_collection.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -60,8 +60,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   tr.refine_global(2);
 
index 9fdda0c54975cce92631bb6bad5aa48e1e476e05..42fb80ada10176efafe60fb2b0694ea8c0d9ffcd 100644 (file)
@@ -21,7 +21,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/hp/dof_handler.h>
 #include <deal.II/hp/fe_collection.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -118,8 +118,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   tr.refine_global(2);
 
index a4e7e9e4a89960e4df9b0eba6c8eab14a9106be8..9f06c311aded2000f93f104d6328c480f3dd9412 100644 (file)
@@ -22,7 +22,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/hp/dof_handler.h>
 #include <deal.II/dofs/dof_renumbering.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -70,8 +70,8 @@ void test_hyper_cube()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   tr.refine_global(1);
 
index 152aef1b23fbec7cefeaa53bc497ca900394fcd1..114b91766f892e70379f44cc04f29de2d7cca09d 100644 (file)
@@ -23,7 +23,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/lac/constraint_matrix.h>
 #include <deal.II/fe/fe_q.h>
@@ -61,8 +61,8 @@ void test_hyper_sphere()
   Triangulation<dim> tr;
   GridGenerator::hyper_ball(tr);
 
-  static const HyperBallBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
 
   tr.refine_global(2);
 
index 13a3d0a59283d87b729dcb3443984586b66dd4d1..5bb0cd2ec2f490e1060b824d8a8b6fd143a2522e 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_tools.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/hp/dof_handler.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -63,14 +63,14 @@ public:
 template<int dim>
 void test()
 {
-  HyperBallBoundary<dim> boundary;
+  PolarManifold<dim> boundary;
 
   // create 2 triangulations with the
   // same coarse grid, and refine
   // them differently
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  tria.set_boundary (0, boundary);
+  tria.set_manifold (0, boundary);
 
 
   FE_Q<dim> fe(1);
index 7fda38c7ae8818fa8974f088a36a1270838d7189..3c57a88d24397a95635544337dd7913f8b954160 100644 (file)
@@ -24,7 +24,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/logstream.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/grid_generator.h>
@@ -86,7 +86,7 @@ void test ()
   deallog << "dim=" << dim << std::endl;
 
   Triangulation<dim> tria;
-  StraightBoundary<dim> boundary;
+  FlatManifold<dim> boundary;
 
   for (unsigned int case_no=0; case_no<2; ++case_no)
     {
index 9e42fc3a6c7900c945330aa07744a2b08b63ca22..bdfabddcf95427d77fba02330119c511d4250a8b 100644 (file)
@@ -24,7 +24,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/logstream.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/grid_generator.h>
@@ -86,7 +86,7 @@ void test ()
   deallog << "dim=" << dim << std::endl;
 
   Triangulation<dim> tria;
-  StraightBoundary<dim> boundary;
+  FlatManifold<dim> boundary;
 
   for (unsigned int case_no=0; case_no<2; ++case_no)
     {
index f99a6de525510ada54786841a6869e3ab67521fc..7ee11c5231d686b5da07db96b2ea2a06be3bb75a 100644 (file)
@@ -24,7 +24,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/logstream.h>
 #include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/manifold.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/grid_generator.h>
@@ -40,7 +40,7 @@ void test ()
   deallog << "dim=" << dim << std::endl;
 
   Triangulation<dim> tria;
-  StraightBoundary<dim> boundary;
+  FlatManifold<dim> boundary;
 
   GridGenerator::hyper_cube(tria, 0, 1);
 
index a775f42e01579700cf52cc89707facbb076981f7..f3c349b0f73e1c4fcfa84cb49defd02efd2537ba 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index f3abedbaa3d82db2679c0a8142a2ed63ffc066c1..853870ea9fa3be0b827e9538516669bfe8ee086a 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 560e370884bb8a8ea2c675ce231c6294882aa417..e9df1e8eb0782cbefe54cc8f0f16cb07601107e6 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 06a44e2fa38c0bedb2ff6d285f7cbeeb42bfc361..0563c9b404726372bee5695cf631e78a4a231489 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 92f3e798a935603ec764baf69cb912c5c3bb917a..84007365d9339b07e7f4b373250045e7dd3076e2 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index ea25342768957b848fe4c4f1518ed08d2c2022c0..a901b9e4a8b6780513b2b644f52730c0d7d4ad70 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index d68ce7406bccad5bf9961a009dcb0a159413169d..3eb1c9760e8d4098bb4b07a62188bd34bdc53b63 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 5ec35b03d1c82ab0eec17a264515dc49876dc9fb..abb4f77398a6e42c54b4fb81d860694c629befeb 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index f168c575efed2a04e157da1224d6ad4816a02eac..7eaafa9eac27eebecec2754be39a6e5e057836bf 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index d595473a6dbe99844a04e6ee53bd949019ea04a2..a274ab7503848102c173310753b4e09a60f7c2e4 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 2d17109e8b65128a14513445d1cba22dc687cf9e..b9acb42c19e6d8f2d78371a73defe090d07feb72 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index c92552b924f77b9f8b683b81143f2a1021f3aba3..e15c9eccd34418c0a4d9828022bcf0bd8fa3ed5e 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 667829ef93366579e73e5b674bb7fadb2b05fb04..f704a30cff99cc3d94db321bbe8edfdada1a6d66 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index cc293f730e930e6b2cc9e98f91a292db059b7805..2b409d59bc2d39959c55bfd5e596808231102cc9 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index d0d8b7308ded1c11382f3e3b966e57a09ca22eca..6f6b6383a785232776f5efcd386afaec688ac8b1 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 320a751f4f19a1563b54c8dd6a2f9886d227f5a7..e16816757871cffd8af05a41ff55694d601201d8 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 2e7ea87b06d1a11459c54a2fc7882653b1977bcd..4ecbfe877beb840935a3228898ebef4f95a0d2ab 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 5a7a1c0241650702028e05970d2375b02c23e113..14a208626857a9ca5e8ae8ca00d7718414ddd2e1 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 70fc0c2cd2028eef63f50bc995817104e3153241..742e259968d3c8b9c4f5aa24c4af5e26c5c9d4c6 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 2674c78f0772bf9ef350bb75f3cb1172ff6632d0..8ba9536a05abc8b1b2c0bb2e195ca1baa079c3f0 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_monomial.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 9f8408ecf8debf43e1385e78091b5f9eca0c2328..195619b92ebab19567e09b33c523fbea3bc1d924 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index a0cedea1dce99b5ef538f1dc7a5f76cb3b4f7e21..cf446f5ffcc75f307646482c3c1860cadc321b38 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index c9b59764137c895c53e31206e72d99c6ee5fbc51..faeb144e61e015bb2bdb330ca210035f631dc4da 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 3a61ff3d760d790e903bcd349cd9ff6a07e319d0..2dd9fda0aef0da21c02dd2c55b20b767ff9d9a78 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 0a97684c97b199cceffe1d9f66e8829bc438b952..c824830b5c43ddcc928ba10fac37ddfb1549d3e4 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index e8524ccbf470f6958d46ca0471078f4688d416e3..f670050322412033592da01f2cb3898c33c4ccbc 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 956719459cbfcc613dbf37c638db2fb509b3aafc..87728dd148dbd3cae5f4ff4bce09288e0fd0fc37 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 75e708153b5ccff2df9d2a56285462222d8356e7..f0fe0a5150aea239063fe1971ba0d77cac14334c 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 79262451cb4107e8483a7f6fa7430488a698d44d..feac5e995bdb1dce56092e6b8db4697f99048c3f 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 73262ee96512e29264df1a429b635b0b3137708d..3d82e22eed860c68bb76573106dd4a27dd39af7d 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgp_nonparametric.h>
@@ -152,13 +152,13 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
-  tr.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary(tr.begin_active()->center());
+  tr.set_manifold (1, boundary);
 
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_indicator (1);
+      tr.begin_active()->face(f)->set_manifold_indicator (1);
 
   test(tr);
 }
index 4351b190548954c63817b23ff7501356750a0a56..fe77dd4677b7bfd5ec67bac8744632e6eb089ed4 100644 (file)
 #include <deal.II/fe/mapping_q.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 template<int dim>
 void test()
 {
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
 
   MappingQ<dim> mapping(3);
   FE_Nothing<dim> dummy;
index d646edf3aa9bce6c5d29e0f14c9a3059b029c593..80cfdc3e324f0aaee59946dd5ffcafcc875fe9f8 100644 (file)
@@ -26,7 +26,7 @@
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/mapping_cartesian.h>
 #include <deal.II/fe/mapping_q1.h>
 #include <deal.II/fe/mapping_q.h>
@@ -306,10 +306,10 @@ void create_triangulations(std::vector<Triangulation<2> *> &tria_ptr,
       v2(1) = 3.;
       v3(0) = 3.;
       v3(1) = 3.;
-      tria->set_boundary(1,*boundary1);
-      tria->set_boundary(2,*boundary2);
-      tria->begin_active()->face(1)->set_boundary_indicator(1);
-      tria->begin_active()->face(3)->set_boundary_indicator(2);
+      tria->set_manifold(1,*boundary1);
+      tria->set_manifold(2,*boundary2);
+      tria->begin_active()->face(1)->set_manifold_indicator(1);
+      tria->begin_active()->face(3)->set_manifold_indicator(2);
       double pi=std::acos(-1.);
       double alpha=2*std::atan(0.5);
       exact_areas.push_back(4+pi-2.5*(alpha-std::sin(alpha)));
@@ -346,8 +346,8 @@ void create_triangulations(std::vector<Triangulation<2> *> &tria_ptr,
       v2(1) = 1.;
       v3(0) = 0.5;
       v3(1) = 1.5;
-      tria->set_boundary(1,*boundary1);
-      tria->begin_active()->face(1)->set_boundary_indicator(1);
+      tria->set_manifold(1,*boundary1);
+      tria->begin_active()->face(1)->set_manifold_indicator(1);
       exact_areas.push_back(0.);
       for (unsigned int i=0; i<=4; ++i)
         show[4][i]=1;
@@ -401,8 +401,8 @@ void create_triangulations(std::vector<Triangulation<3> *> &tria_ptr,
       tria=new Triangulation<3>();
       tria_ptr.push_back(tria);
       GridGenerator::hyper_cube(*tria, 1., 3.);
-      tria->set_boundary(1,*boundary1);
-      tria->begin_active()->face(1)->set_boundary_indicator(1);
+      tria->set_manifold(1,*boundary1);
+      tria->begin_active()->face(1)->set_manifold_indicator(1);
       exact_areas.push_back(8.+pi/3*h*h*(3*r-h));
     }
 
index 9c9af6b097d2806572be94f5ae1f0e5ba73d7ae6..0df81eebdc667042d8ab39ec92c3a8d25f97f4f6 100644 (file)
@@ -26,7 +26,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/fe/mapping_q1.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/numerics/fe_field_function.h>
 
 template<int dim>
 void test2()
 {
-  const HyperBallBoundary<dim> boundary_description;
+  const PolarManifold<dim> boundary_description;
 
   Triangulation<dim>   triangulation;
   GridGenerator::hyper_ball (triangulation);
-  triangulation.set_boundary (0, boundary_description);
+  triangulation.set_manifold (0, boundary_description);
   triangulation.refine_global (1);
   MappingQ1< dim > mapping;
 
@@ -70,7 +70,7 @@ void test2()
   points[95]=p;
   fe_function.value_list (points, m); // <<<< this fails at point[95] but only if the other points are filled in?!
 
-  triangulation.set_boundary (0);
+  triangulation.set_manifold (0);
 
   deallog << "OK" << std::endl;
 }
index 3eb7e558cf48a58133e26d1ee8438a7ac874b8d0..666b9c4682068e501d8657bd375beb56423c6fcc 100644 (file)
@@ -30,7 +30,7 @@
 #include <deal.II/base/utilities.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/mapping_q.h>
 
 
@@ -50,8 +50,8 @@ void test_real_to_unit_cell()
 
   // set the boundary indicator for
   // one face of the single cell
-  triangulation.set_boundary (1, boundary);
-  triangulation.begin_active()->face(0)->set_boundary_indicator (1);
+  triangulation.set_manifold (1, boundary);
+  triangulation.begin_active()->face(0)->set_manifold_indicator (1);
 
   const unsigned int n_points = 5;
   std::vector< Point<dim> > unit_points(Utilities::fixed_power<dim>(n_points));
index fd2c553cf41d3ba9db997f6fb9bb9b4dc78a39b3..9c1b8d09a9d105c9c7fcd29eada063aa11c79df1 100644 (file)
@@ -30,7 +30,7 @@
 #include <deal.II/base/utilities.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/mapping_q.h>
 
 
@@ -50,8 +50,8 @@ void test_real_to_unit_cell()
 
   // set the boundary indicator for
   // one face of the single cell
-  triangulation.set_boundary (1, boundary);
-  triangulation.begin_active()->face(0)->set_boundary_indicator (1);
+  triangulation.set_manifold (1, boundary);
+  triangulation.begin_active()->face(0)->set_manifold_indicator (1);
 
   const unsigned int n_points = 5;
   std::vector< Point<dim> > unit_points(Utilities::fixed_power<dim>(n_points));
index fb667f8e921335eb7503d01496c0666dcd329c92..dacae7c7204f4e1598ee47d95aad8b2ccacdb49b 100644 (file)
@@ -35,7 +35,7 @@
 #include <deal.II/base/utilities.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/mapping_q.h>
 
 
@@ -46,7 +46,7 @@ void test_real_to_unit_cell()
   // define a spherical cap boundary
   // to be used for one of the faces
   const double radius = Point<dim>(1.43757e-10, 4.48023e+06, 4.48023e+06).norm();
-  HyperBallBoundary<dim> boundary (Point<dim>(), radius);
+  PolarManifold<dim> boundary (Point<dim>(), radius);
 
   // create the mesh: a single cell
   // with the following coordinates:
@@ -77,7 +77,7 @@ void test_real_to_unit_cell()
   // set the boundary indicator for
   // one face and adjacent edges of
   // the single cell
-  triangulation.set_boundary (1, boundary);
+  triangulation.set_manifold (1, boundary);
   triangulation.begin_active()->face(5)->set_all_boundary_indicators (1);
 
   // now try to find the coordinates
index 8af5d89bc982befd26633bc6d664c39c0c971487..e01f579e2e351b9a8df4766220a4fb5402242c29 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/utilities.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/mapping_q.h>
 
 
@@ -36,7 +36,7 @@ void test_real_to_unit_cell()
   // define a spherical cap boundary
   // to be used for one of the faces
   const double radius = Point<dim>(1.43757e-10, 4.48023e+06, 4.48023e+06).norm();
-  HyperBallBoundary<dim> boundary (Point<dim>(), radius);
+  PolarManifold<dim> boundary (Point<dim>(), radius);
 
   // create the mesh: a single cell
   // with the following coordinates:
@@ -67,7 +67,7 @@ void test_real_to_unit_cell()
   // set the boundary indicator for
   // one face and adjacent edges of
   // the single cell
-  triangulation.set_boundary (1, boundary);
+  triangulation.set_manifold (1, boundary);
   triangulation.begin_active()->face(5)->set_all_boundary_indicators (1);
 
   // now try to find the coordinates
index a96c6bc26c37daa6f2167f50bd04fd76c1709e5e..45b4f613d66e041556a439ed97c8b69743b23950 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/utilities.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/mapping_q.h>
 
 
@@ -36,7 +36,7 @@ void test_real_to_unit_cell()
   // define a spherical cap boundary
   // to be used for one of the faces
   const double radius = Point<dim>(1.43757e-10, 4.48023e+06, 4.48023e+06).norm();
-  HyperBallBoundary<dim> boundary (Point<dim>(), radius);
+  PolarManifold<dim> boundary (Point<dim>(), radius);
 
   // create the mesh: a single cell
   // with the following coordinates:
@@ -67,7 +67,7 @@ void test_real_to_unit_cell()
   // set the boundary indicator for
   // one face and adjacent edges of
   // the single cell
-  triangulation.set_boundary (1, boundary);
+  triangulation.set_manifold (1, boundary);
   triangulation.begin_active()->face(5)->set_all_boundary_indicators (1);
 
   // now try to find the coordinates
index 3a839f676494f4ec9cee612b53890c3b8c2a170d..e106281c382fd75d49dedd706229d34f9cdbccfd 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/base/utilities.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/fe/mapping_q.h>
 
 
@@ -36,7 +36,7 @@ void test_real_to_unit_cell()
   // define a spherical cap boundary
   // to be used for one of the faces
   const double radius = Point<dim>(1.43757e-10, 4.48023e+06, 4.48023e+06).norm();
-  HyperBallBoundary<dim> boundary (Point<dim>(), radius);
+  PolarManifold<dim> boundary (Point<dim>(), radius);
 
   // create the mesh: a single cell
   // with the following coordinates:
@@ -67,7 +67,7 @@ void test_real_to_unit_cell()
   // set the boundary indicator for
   // one face and adjacent edges of
   // the single cell
-  triangulation.set_boundary (1, boundary);
+  triangulation.set_manifold (1, boundary);
   triangulation.begin_active()->face(5)->set_all_boundary_indicators (1);
 
   // now see if the point is inside
index d2f29ff513bf5258be78e5acfcc190559e16d736..99749476e2031819e74d6422b7943c9d322b55da 100644 (file)
@@ -24,7 +24,7 @@
 #include <deal.II/base/tensor.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 
 #include <fstream>
@@ -47,7 +47,7 @@ void test(std::ostream &out)
   if (dim>1) p3[1] = 1.;
   if (dim>2) p3[2] = 4.;
 
-  HalfHyperBallBoundary<dim> boundary_description (p1, 3);
+  HalfPolarManifold<dim> boundary_description (p1, 3);
   GridOut go;
   GridOut::OutputFormat format = GridOut::gnuplot;
 
@@ -55,7 +55,7 @@ void test(std::ostream &out)
     deallog << "half_hyper_ball" << std::endl;
     Triangulation<dim> tr;
     GridGenerator::half_hyper_ball(tr, p1, 3.);
-    tr.set_boundary (0, boundary_description);
+    tr.set_manifold (0, boundary_description);
 
     tr.refine_global (2);
     go.write(tr, out, format);
index 2bfed59373b498c05d682ca3ffd2dd24b7c845fd..cb3d056b2da949f7b96a47fb5320ca96c5a6a59b 100644 (file)
@@ -25,7 +25,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_tools.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/hp/dof_handler.h>
 #include <deal.II/hp/q_collection.h>
@@ -64,14 +64,14 @@ public:
 template<int dim>
 void test()
 {
-  HyperBallBoundary<dim> boundary;
+  PolarManifold<dim> boundary;
 
   // create 2 triangulations with the
   // same coarse grid, and refine
   // them differently
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  tria.set_boundary (0, boundary);
+  tria.set_manifold (0, boundary);
 
 
   hp::FECollection<dim> fe;
index ed6f02f7cf387ad608e7ba6915179eb9389517bb..330075c526d1a3dad43135ac738f922ad0741d9f 100644 (file)
@@ -30,7 +30,7 @@ std::ofstream logfile("output");
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/hp/dof_handler.h>
@@ -56,8 +56,8 @@ void gnuplot_output()
 
   Triangulation<dim> triangulation;
   GridGenerator::hyper_ball (triangulation);
-  static const HyperBallBoundary<dim> boundary;
-  triangulation.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  triangulation.set_manifold (0, boundary);
 
   for (unsigned int refinement=0; refinement<2;
        ++refinement, triangulation.refine_global(1))
@@ -99,8 +99,8 @@ void compute_pi_by_area ()
       Triangulation<dim> triangulation;
       GridGenerator::hyper_ball (triangulation);
 
-      static const HyperBallBoundary<dim> boundary;
-      triangulation.set_boundary (0, boundary);
+      static const PolarManifold<dim> boundary;
+      triangulation.set_manifold (0, boundary);
 
       MappingQ<dim> m(degree);
       const hp::MappingCollection<dim> mapping (m);
@@ -166,8 +166,8 @@ void compute_pi_by_perimeter ()
       Triangulation<dim> triangulation;
       GridGenerator::hyper_ball (triangulation);
 
-      static const HyperBallBoundary<dim> boundary;
-      triangulation.set_boundary (0, boundary);
+      static const PolarManifold<dim> boundary;
+      triangulation.set_manifold (0, boundary);
 
       const MappingQ<dim> m (degree);
       const hp::MappingCollection<dim> mapping(m);
index 994d146876c4f69a6f6095e55212f7c121d1f16a..ab02f94c779cf4e92194041b9dd39d1832899896 100644 (file)
@@ -35,7 +35,7 @@ std::ofstream logfile("output");
 #include <deal.II/lac/precondition.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/dofs/dof_handler.h>
@@ -197,8 +197,8 @@ template <int dim>
 void LaplaceProblem<dim>::run ()
 {
   GridGenerator::hyper_ball (triangulation);
-  static const HyperBallBoundary<dim> boundary;
-  triangulation.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  triangulation.set_manifold (0, boundary);
 
   for (unsigned int cycle=0; cycle<6; ++cycle, triangulation.refine_global(1))
     {
index d0c6e36b2f759c5f05be2db49e2d0595ab9b1b48..5169a69e76783a40434c8995cb643d0e64240380 100644 (file)
@@ -38,7 +38,7 @@ std::ofstream logfile("output");
 #include <deal.II/lac/precondition.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/dofs/dof_handler.h>
@@ -200,8 +200,8 @@ template <int dim>
 void LaplaceProblem<dim>::run ()
 {
   GridGenerator::hyper_ball (triangulation);
-  static const HyperBallBoundary<dim> boundary;
-  triangulation.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  triangulation.set_manifold (0, boundary);
 
   for (unsigned int cycle=0; cycle<6; ++cycle, triangulation.refine_global(1))
     {
index 3bff1ee6187a04d0424d758e30f99f022eb795e5..2964877e6082a2cafa77fa45826666248991f8a9 100644 (file)
@@ -47,7 +47,7 @@ std::ofstream logfile("output");
 
 #include <deal.II/grid/grid_in.h>
 
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 #include <fstream>
 #include <sstream>
@@ -304,8 +304,8 @@ void LaplaceProblem<dim>::run ()
 
           grid_in.read_ucd (input_file);
 
-          static const HyperBallBoundary<dim> boundary;
-          triangulation.set_boundary (0, boundary);
+          static const PolarManifold<dim> boundary;
+          triangulation.set_manifold (0, boundary);
         }
 
       deallog << "   Number of active cells: "
index ed805aa9cd262112daebbf4def1b6a93ce554790..2a34966443210c3b7dd9b629a4eaa45cd184f20e 100644 (file)
@@ -38,7 +38,7 @@ std::ofstream logfile("output");
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/hp/fe_values.h>
@@ -337,8 +337,8 @@ void LaplaceProblem<dim>::run ()
         {
           GridGenerator::hyper_ball (triangulation);
 
-          static const HyperBallBoundary<dim> boundary;
-          triangulation.set_boundary (0, boundary);
+          static const PolarManifold<dim> boundary;
+          triangulation.set_manifold (0, boundary);
 
           triangulation.refine_global (1);
         }
index eab595cb6fa8c3817610e5afc8695d6b1eb18890..79a9f4791f64dfbf32260a436ae6f6d7ba3e7516 100644 (file)
@@ -34,8 +34,8 @@ void test ()
 {
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   // refine first and last cell
   tria.begin(tria.n_levels()-1)->set_refine_flag();
   tria.last()->set_refine_flag();
index a5b2dd1b90c2f0cc3f1fe810a8664f6cfbe45098..350ab0f43839dd172e1714dc173e8194a3cb7ec4 100644 (file)
@@ -91,8 +91,8 @@ void test ()
   typedef double number;
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   // refine first and last cell
   tria.begin(tria.n_levels()-1)->set_refine_flag();
   tria.last()->set_refine_flag();
index 866eea77c355bf39aaecddb3c7a5beb8d87b2703..dfa41da404d77aa3a1c1e4237caf209042047713 100644 (file)
@@ -35,8 +35,8 @@ void test ()
   typedef double number;
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   // refine first and last cell
   tria.begin(tria.n_levels()-1)->set_refine_flag();
   tria.last()->set_refine_flag();
index 01601f2944368d11b0db341e80ea43a3a36064ac..760f5956b99de5174a9e1bd87ad6d311889ead55 100644 (file)
@@ -94,8 +94,8 @@ void test ()
   typedef double number;
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   // refine first and last cell
   tria.begin(tria.n_levels()-1)->set_refine_flag();
   tria.last()->set_refine_flag();
index 881125ad796eba688cc86ad8e6ee92bda5e9b4e6..13af19b4d793e523b2022bc868fe1cc55fcb232f 100644 (file)
@@ -56,8 +56,8 @@ void test ()
   typedef double number;
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   // refine first and last cell
   tria.begin(tria.n_levels()-1)->set_refine_flag();
   tria.last()->set_refine_flag();
index 378006b9e15ed76023561e247ea51dd45ed97400..57e9c9b2345b8895a23fd9a33105a99e926fd849 100644 (file)
@@ -31,7 +31,7 @@
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -143,9 +143,9 @@ void test ()
 {
   Triangulation<dim> tria;
   GridGenerator::hyper_shell (tria, Point<dim>(), 1., 2., 96, true);
-  static const HyperShellBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
-  tria.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
+  tria.set_manifold (1, boundary);
 
   // refine a few cells
   for (unsigned int i=0; i<11-3*dim; ++i)
index 313a249d87c9b33796bbda3b9b4be619169b3499..51ce6406b82fc77139f1a4d78cbc8f9b5505fec4 100644 (file)
@@ -30,7 +30,7 @@
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -148,8 +148,8 @@ void test ()
   typedef double number;
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   typename Triangulation<dim>::active_cell_iterator
   cell = tria.begin_active (),
   endc = tria.end();
index 2f3f7fe0bb9b3d11d5ac2e481253f1a3d9c4ab78..8797ba311e94b05e87a37226db5fbe636e825775 100644 (file)
@@ -31,7 +31,7 @@
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -229,8 +229,8 @@ void test ()
   // cells
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   typename Triangulation<dim>::active_cell_iterator
   cell = tria.begin_active (),
   endc = tria.end();
index 36b71b3aaa5546c8a5e670ee43eef50a3975356c..0fdc5a1e62ec08fbfde279d3aecb56c206e1b35b 100644 (file)
@@ -34,9 +34,9 @@ void test ()
   Triangulation<dim> tria;
   GridGenerator::hyper_shell (tria, Point<dim>(),
                               0.5, 1., 96, true);
-  static const HyperShellBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
-  tria.set_boundary (1, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
+  tria.set_manifold (1, boundary);
   if (dim == 2)
     tria.refine_global (2);
 
index 663b35ff9824216378ca2ae021daf47568dea188..20bd36a4bd7d48b7d1cdcfd774e25c94c2b1d984 100644 (file)
@@ -38,8 +38,8 @@ void test ()
 {
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   typename Triangulation<dim>::active_cell_iterator
   cell = tria.begin_active (),
   endc = tria.end();
index 2993a1f2baa233e9c36a223b36154141ef86f07b..f67eea6757d841974a9f3754fa6f8c75a6ebcb9e 100644 (file)
@@ -34,7 +34,7 @@ std::ofstream logfile("output");
 #include <deal.II/lac/block_vector.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/dofs/dof_renumbering.h>
@@ -107,8 +107,8 @@ void test ()
 {
   Triangulation<dim>   tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   tria.refine_global(4-dim);
 
   // refine a few cells
index a56cfd84b5d5f8643f50f67a01368129861d0ece..f3c105a094fd9c68f2e6fa70cb86f17e52865978 100644 (file)
@@ -99,8 +99,8 @@ void test ()
   typedef double number;
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   tria.refine_global(1);
 
   // refine a few cells
index c86f675a4089dccdf940f220f09457430c7ab301..c0ad6c688e440e74d3c4cde7f4bf9bb8f59d6adb 100644 (file)
@@ -37,8 +37,8 @@ void test ()
 {
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   tria.refine_global(5-dim);
 
   FE_Q<dim> fe (fe_degree);
index c4cfd300ae7e5b4482cb6f725d76b9c7afa3e48b..759e6265a7ffc48cc75c6e37eb77eaee242e8afa 100644 (file)
@@ -34,7 +34,7 @@ std::ofstream logfile("output");
 #include <deal.II/lac/block_vector.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/dofs/dof_renumbering.h>
@@ -125,9 +125,9 @@ void test ()
   Triangulation<dim>   triangulation;
   GridGenerator::hyper_shell (triangulation, Point<dim>(),
                               0.5, 1., 96, true);
-  static HyperShellBoundary<dim> boundary;
-  triangulation.set_boundary (0, boundary);
-  triangulation.set_boundary (1, boundary);
+  static PolarManifold<dim> boundary;
+  triangulation.set_manifold (0, boundary);
+  triangulation.set_manifold (1, boundary);
   triangulation.begin_active()->set_refine_flag();
   triangulation.last()->set_refine_flag();
   triangulation.execute_coarsening_and_refinement();
index 3c2af197fccc697df80bad9f92b51885bf1b3b91..82018416e179cd716d8307878a16210eab0e4a0c 100644 (file)
@@ -32,7 +32,7 @@ std::ofstream logfile("output");
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -174,8 +174,8 @@ void test ()
 {
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   // refine first and last cell
   tria.begin(tria.n_levels()-1)->set_refine_flag();
   tria.last()->set_refine_flag();
index 206b90ef730cd41e6dd0e26aff6e2177f32824a6..245b651716ad30cc502f019497adaf76424fd10c 100644 (file)
@@ -28,7 +28,7 @@
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/lac/constraint_matrix.h>
@@ -48,8 +48,8 @@ void test ()
   typedef double number;
   Triangulation<dim> tria;
   GridGenerator::hyper_ball (tria);
-  static const HyperBallBoundary<dim> boundary;
-  tria.set_boundary (0, boundary);
+  static const PolarManifold<dim> boundary;
+  tria.set_manifold (0, boundary);
   tria.refine_global(5-dim);
 
   MappingQ<dim> mapping (4);
index 95d256af9ed16f01c283118e356fe30c5cdbbeed..a90a20f392973ff72d41621127dd1f758539873f 100644 (file)
@@ -36,7 +36,7 @@
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_refinement.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/dofs/dof_tools.h>
@@ -463,8 +463,8 @@ namespace Step50
           {
             GridGenerator::hyper_cube (triangulation);
 
-            // static const HyperBallBoundary<dim> boundary;
-            // triangulation.set_boundary (0, boundary);
+            // static const PolarManifold<dim> boundary;
+            // triangulation.set_manifold (0, boundary);
 
             triangulation.refine_global (2);
           }
index 50af594234927564f4c1beba567b6b8bd66233d0..4181eddc8e34d88f347f489cdbd2d6460ff96218 100644 (file)
@@ -36,7 +36,7 @@
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_refinement.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/dofs/dof_tools.h>
@@ -456,8 +456,8 @@ namespace Step50
           {
             GridGenerator::hyper_cube (triangulation);
 
-            // static const HyperBallBoundary<dim> boundary;
-            // triangulation.set_boundary (0, boundary);
+            // static const PolarManifold<dim> boundary;
+            // triangulation.set_manifold (0, boundary);
 
             triangulation.refine_global (3);
           }
index c7bab7f9a3ce572d62d1aa2f6e362281adb5c480..badf5d2492b732e5737788f7fb5de536b603dadc 100644 (file)
@@ -38,7 +38,7 @@
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/dofs/dof_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/numerics/data_out.h>
 #include <deal.II/lac/trilinos_vector.h>
 
@@ -111,9 +111,9 @@ void test()
                               R1,
                               12,
                               true);
-  static HyperShellBoundary<dim> boundary;
-  tr.set_boundary (0, boundary);
-  tr.set_boundary (1, boundary);
+  static PolarManifold<dim> boundary;
+  tr.set_manifold (0, boundary);
+  tr.set_manifold (1, boundary);
 
   tr.refine_global (1);
   for (unsigned int step=0; step<20; ++step)
@@ -218,8 +218,8 @@ void test()
       cat_file((std::string("dat.") + Utilities::int_to_string(0)).c_str());
     }
 
-  tr.set_boundary (0);
-  tr.set_boundary (1);
+  tr.set_manifold (0);
+  tr.set_manifold (1);
 }
 
 
index d5940f956f3066fb8475ea8da3f4c3543cce1650..d3a30412f11720a7750e65d7a077a995a43c4a63 100644 (file)
@@ -39,7 +39,7 @@
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/dofs/dof_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/numerics/data_out.h>
 #include <deal.II/lac/trilinos_vector.h>
 
@@ -147,9 +147,9 @@ void test()
                               R1,
                               12,
                               true);
-  static HyperShellBoundary<dim> boundary;
-//  tr.set_boundary (0, boundary);
-// tr.set_boundary (1, boundary);
+  static PolarManifold<dim> boundary;
+//  tr.set_manifold (0, boundary);
+// tr.set_manifold (1, boundary);
 
   tr.refine_global (3);
   if (1)
@@ -311,8 +311,8 @@ void test()
   std::ofstream output (filename.c_str());
   data_out.write_deal_II_intermediate (output);
 
-  tr.set_boundary (0);
-  tr.set_boundary (1);
+  tr.set_manifold (0);
+  tr.set_manifold (1);
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
     deallog << "OK" << std::endl;
 }
index b30e4008e3efbd4624b8f03506617a176e7f0d81..1c5465f5f338af1fb035a4ca94bc55ed77e49948 100644 (file)
@@ -38,7 +38,7 @@
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/dofs/dof_accessor.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/numerics/data_out.h>
 #include <deal.II/lac/trilinos_block_vector.h>
 
@@ -146,9 +146,9 @@ void test()
                               R1,
                               12,
                               true);
-  static HyperShellBoundary<dim> boundary;
-//  tr.set_boundary (0, boundary);
-//tr.set_boundary (1, boundary);
+  static PolarManifold<dim> boundary;
+//  tr.set_manifold (0, boundary);
+//tr.set_manifold (1, boundary);
 
   tr.refine_global (1);
   if (1)
@@ -306,8 +306,8 @@ void test()
   std::ofstream output (filename.c_str());
   data_out.write_vtu (output);
 
-  tr.set_boundary (0);
-  tr.set_boundary (1);
+  tr.set_manifold (0);
+  tr.set_manifold (1);
 
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
     deallog << "OK" << std::endl;
index 66a9e22e6c423b63c65bf4f38d80627f9b1617cf..cece429dadd5ae732ee8255557042a9bcb7e05a0 100644 (file)
@@ -40,7 +40,7 @@
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_refinement.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 #include <deal.II/integrators/laplace.h>
 
@@ -566,8 +566,8 @@ void LaplaceProblem<dim>::run ()
         {
           GridGenerator::hyper_cube (triangulation, -1, 1);
 
-//          static const HyperBallBoundary<dim> boundary;
-//          triangulation.set_boundary (0, boundary);
+//          static const PolarManifold<dim> boundary;
+//          triangulation.set_manifold (0, boundary);
 
           triangulation.refine_global (1);
         }
index ab5c4b41dbccf9e7779e0ac7f71e837f4f199969..1e5dcad9e82b9b0fc667a5c3d90ace46c3c08bff 100644 (file)
@@ -40,7 +40,7 @@
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_refinement.h>
-#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/manifold_lib.h>
 
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/dofs/dof_tools.h>
@@ -561,8 +561,8 @@ void LaplaceProblem<dim>::run ()
         {
           GridGenerator::hyper_ball (triangulation);
 
-          static const HyperBallBoundary<dim> boundary;
-          triangulation.set_boundary (0, boundary);
+          static const PolarManifold<dim> boundary;
+          triangulation.set_manifold (0, boundary);
 
           triangulation.refine_global (1);
         }

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.