]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Simplify manifold use in matrix-free tests 18165/head
authorMartin Kronbichler <martin.kronbichler@rub.de>
Tue, 25 Feb 2025 08:10:41 +0000 (09:10 +0100)
committerMartin Kronbichler <martin.kronbichler@rub.de>
Tue, 25 Feb 2025 08:12:20 +0000 (09:12 +0100)
38 files changed:
tests/matrix_free/assemble_matrix_01.cc
tests/matrix_free/assemble_matrix_02.cc
tests/matrix_free/assemble_matrix_03.cc
tests/matrix_free/copy_feevaluation.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_notempl.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_hessians.cc
tests/matrix_free/integrate_functions_multife.cc
tests/matrix_free/inverse_mass_01.cc
tests/matrix_free/inverse_mass_03.cc
tests/matrix_free/inverse_mass_05.cc
tests/matrix_free/inverse_mass_06.cc
tests/matrix_free/inverse_mass_08.cc
tests/matrix_free/matrix_vector_04.cc
tests/matrix_free/matrix_vector_04b.cc
tests/matrix_free/matrix_vector_07.cc
tests/matrix_free/matrix_vector_14.cc
tests/matrix_free/matrix_vector_14_notempl.cc
tests/matrix_free/matrix_vector_14b.cc
tests/matrix_free/matrix_vector_15.cc
tests/matrix_free/matrix_vector_16.cc
tests/matrix_free/matrix_vector_17.cc
tests/matrix_free/matrix_vector_curl.cc
tests/matrix_free/matrix_vector_hp.cc
tests/matrix_free/matrix_vector_hp_no_template.cc
tests/matrix_free/matrix_vector_hp_parallel_01.cc
tests/matrix_free/matrix_vector_mg.cc
tests/matrix_free/matrix_vector_stokes_noflux.cc
tests/matrix_free/matrix_vector_stokes_notempl.cc
tests/matrix_free/matrix_vector_stokes_onedof.cc
tests/matrix_free/no_index_initialize.cc
tests/matrix_free/quadrature_points.cc

index 5bdf3b40577381e4c1e53576800748ab3ca469ab..d786bdd202d76bbdd856b324e1410d8272ebdd6b 100644 (file)
@@ -24,7 +24,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -128,16 +127,8 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
index f6727aac21b337860c21a34dfea40a538299f1d4..8b17dc02d3e9c3dce9e17536425a732552af6460 100644 (file)
@@ -24,7 +24,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -184,16 +183,8 @@ template <int dim>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   if (dim == 2)
     tria.refine_global(1);
index dc4c5043e26164eedd5d1dd379ec61ec3752ebb3..eb6b812f72dfbb75b18b5b97a7a46fa0a2f50c77 100644 (file)
@@ -27,7 +27,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -168,16 +167,8 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   tria.refine_global(5 - dim);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
index 85acdf22fa48874922565cdefff89d5cf13dffed..1622dc5f180720dce641259b2158cad27ec7f192 100644 (file)
@@ -31,7 +31,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -131,11 +130,8 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  SphericalManifold<dim> manifold;
-  Triangulation<dim>     triangulation;
+  Triangulation<dim> triangulation;
   GridGenerator::hyper_shell(triangulation, Point<dim>(), 0.5, 1., 96, true);
-  triangulation.set_all_manifold_ids(0);
-  triangulation.set_manifold(0, manifold);
 
   triangulation.begin_active()->set_refine_flag();
   triangulation.last()->set_refine_flag();
index b6f5aed05208336840b22c075dbcc0b4f4092eb8..73da0f3ee15b1fec31e49895ddb668279cb93919 100644 (file)
@@ -29,16 +29,8 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   // refine first and last cell
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
index 4330c683c2c2b4a12bc78d0507089fa99ae8581b..b10893f46621f80009eea4cb1ab80de3284ba6c2 100644 (file)
@@ -31,16 +31,8 @@ void
 test()
 {
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   // refine first and last cell
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
index b06e1fcac106ca31699876b3c0b276ca191ffe53..7926f5dd106e5a8ebda6ce54de428d6dc4f2749f 100644 (file)
@@ -30,16 +30,8 @@ void
 test()
 {
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   // refine first and last cell
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
index 56b02c760d80bc24521f2c471c4896d51944c79d..d14ecdc0d7e8aec727db88d4ceed342caff523be 100644 (file)
@@ -52,16 +52,8 @@ void
 test()
 {
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   // refine first and last cell
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
index 62d92148afe6fb2bf848b45240d3bd332aba6b23..b1184380d587bc7a129970a579aa33b180d67d63 100644 (file)
@@ -38,16 +38,8 @@ void
 test()
 {
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   // refine first and last cell
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
index 2bcb688262731639ad33ffae41d18b5338def3bb..a4fbeb9bbb08d2fd198038775e1eebfa4871eb96 100644 (file)
@@ -52,16 +52,8 @@ void
 test()
 {
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   // refine first and last cell
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
index 6eb3d4094d3f15e2ab799473053c58b8c31d8a9a..72a81a13c3fb72ae33a105e1a256a6d662f12f7e 100644 (file)
@@ -28,7 +28,6 @@
 #include <deal.II/fe/fe_values.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -151,22 +150,19 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_shell(tria, Point<dim>(), 1., 2., 96, true);
-  tria.set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   // refine a few cells
   for (unsigned int i = 0; i < 11 - 3 * dim; ++i)
     {
-      typename Triangulation<dim>::active_cell_iterator cell =
-                                                          tria.begin_active(),
-                                                        endc = tria.end();
-      unsigned int counter                                   = 0;
-      for (; cell != endc; ++cell, ++counter)
-        if (counter % (7 - i) == 0)
-          cell->set_refine_flag();
+      unsigned int counter = 0;
+      for (const auto &cell : tria.active_cell_iterators())
+        {
+          if (counter % (7 - i) == 0)
+            cell->set_refine_flag();
+          ++counter;
+        }
       tria.execute_coarsening_and_refinement();
     }
 
index 378af1e13091408caa7aee7a63f64f7387ed1bb2..4f2d7217c5381fe4ba317a0fda1a551909153430 100644 (file)
@@ -27,7 +27,6 @@
 #include <deal.II/fe/fe_values.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -152,24 +151,14 @@ void
 test()
 {
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
-
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 0.2)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
@@ -178,14 +167,15 @@ test()
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
   for (unsigned int i = 0; i < 7 - 2 * dim; ++i)
     {
-      cell                 = tria.begin_active();
       unsigned int counter = 0;
-      for (; cell != endc; ++cell, ++counter)
-        if (counter % (7 - i) == 0)
-          cell->set_refine_flag();
+      for (const auto &cell : tria.active_cell_iterators())
+        {
+          if (counter % (7 - i) == 0)
+            cell->set_refine_flag();
+          ++counter;
+        }
       tria.execute_coarsening_and_refinement();
     }
 
index 0c5c93b02fb9f901f6184dc7a8cadc312b973dcc..54fd25f0a398c9cac92cac691e5802d7416ecd40 100644 (file)
@@ -27,7 +27,6 @@
 #include <deal.II/fe/fe_values.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -168,24 +167,14 @@ void
 test()
 {
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
-
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 0.2)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
@@ -194,14 +183,15 @@ test()
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
   for (unsigned int i = 0; i < 7 - 2 * dim; ++i)
     {
-      cell                 = tria.begin_active();
       unsigned int counter = 0;
-      for (; cell != endc; ++cell, ++counter)
-        if (counter % (7 - i) == 0)
-          cell->set_refine_flag();
+      for (const auto &cell : tria.active_cell_iterators())
+        {
+          if (counter % (7 - i) == 0)
+            cell->set_refine_flag();
+          ++counter;
+        }
       tria.execute_coarsening_and_refinement();
     }
 
index 6cc01f8cb0b6cbfd2a6255240c2967bf5afa6844..47fb6cb94a9ae08ebbbb22967f1b9004e5e38d0e 100644 (file)
@@ -28,7 +28,6 @@
 #include <deal.II/fe/fe_values.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -240,24 +239,14 @@ test()
 {
   // create hyper ball geometry and refine some
   // cells
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
-
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 0.2)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
@@ -266,14 +255,15 @@ test()
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
   for (unsigned int i = 0; i < 7 - 2 * dim; ++i)
     {
-      cell                 = tria.begin_active();
       unsigned int counter = 0;
-      for (; cell != endc; ++cell, ++counter)
-        if (counter % (7 - i) == 0)
-          cell->set_refine_flag();
+      for (const auto &cell : tria.active_cell_iterators())
+        {
+          if (counter % (7 - i) == 0)
+            cell->set_refine_flag();
+          ++counter;
+        }
       tria.execute_coarsening_and_refinement();
     }
 
index 4aa6e562975e91c417fe7f969a4f54917dfc36dc..a2951e18f96cf8f1ad9f488775c5f4b347dd746b 100644 (file)
@@ -25,7 +25,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/precondition.h>
@@ -173,24 +172,15 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index 25869436a09f3681f0d522a02c4064e33e2704b8..cd0f4a0b3a53cd6a95c0b68dffb991ddf547146c 100644 (file)
@@ -26,7 +26,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/precondition.h>
@@ -188,24 +187,15 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index bc51704c03950c33722367aff83122af062848c8..ab2dc3e06bf486a3dda2b73ce137a72a12f14636 100644 (file)
@@ -26,7 +26,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/precondition.h>
@@ -174,24 +173,15 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index b1c4158dbd81d12735dbad72c34773d2b870c08a..07f8fb2330f8ec3d5ccd4c862e4ebf00063fbaba 100644 (file)
@@ -25,7 +25,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/precondition.h>
@@ -165,24 +164,15 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index 0ffa5d859748e8afd57f136dbf8901bf87fa319b..349d626290522e0ee60e369f7e286d9ea24b63a0 100644 (file)
@@ -26,7 +26,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/precondition.h>
@@ -163,14 +162,8 @@ template <int dim>
 void
 test(const unsigned int fe_degree)
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  for (const auto &cell : tria.active_cell_iterators())
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
index 74de2c40e2af182b6b64763f9774435de1c2924f..dd3214687869d60cefb931dd2f7c92db48215900 100644 (file)
@@ -28,11 +28,8 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_shell(tria, Point<dim>(), 0.5, 1., 96, true);
-  tria.set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   if (dim == 2)
     tria.refine_global(2);
 
index f2b06c95bdff1aa88b316aa337b70cead5eecb23..e8d95823cb1a22c526192a1735062728e3f5b196 100644 (file)
@@ -26,11 +26,8 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_shell(tria, Point<dim>(), 0.5, 1., 96, true);
-  tria.set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   if (dim == 2)
     tria.refine_global(2);
 
index cbb2762bea846cbd00111cba806e34fdf78b1349..87e7491fa151722e25f373c03994a754ee80f0ab 100644 (file)
@@ -35,23 +35,13 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 0.2)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
@@ -63,14 +53,15 @@ test()
   tria.execute_coarsening_and_refinement();
   tria.refine_global(4 - dim);
 #endif
-  cell = tria.begin_active();
   for (unsigned int i = 0; i < 10 - 3 * dim; ++i)
     {
-      cell                 = tria.begin_active();
       unsigned int counter = 0;
-      for (; cell != endc; ++cell, ++counter)
-        if (counter % (7 - i) == 0)
-          cell->set_refine_flag();
+      for (const auto &cell : tria.active_cell_iterators())
+        {
+          if (counter % (7 - i) == 0)
+            cell->set_refine_flag();
+          ++counter;
+        }
       tria.execute_coarsening_and_refinement();
     }
 
index 3974da0ab510f728e326e88684b49bf44c42fa85..48636162d09dabe52f961f2786f3a2e630d1a0f7 100644 (file)
@@ -32,23 +32,14 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index 4739dc90462566f75a79ebf839fb9338d1eca15d..55bdc4d2cb159a55e7127358c7569a39aa2c4908 100644 (file)
@@ -29,23 +29,14 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index fa705dc0af242e5a1dc354558156d87e17217e60..95d507b9c25c6f6ca23dad412d6dfb90112bfaf7 100644 (file)
@@ -29,23 +29,14 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index bb70ef75349e180768e6981df789980ebec63cda..cf439fd4499599a4caabc90587a6628a330071a8 100644 (file)
@@ -27,7 +27,6 @@
 #include <deal.II/fe/fe_values.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -183,24 +182,15 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index e6455e9d4387e195ced69aa8be52822178007bf8..b8588888c60ec613195158a3507bbc445ecb6923 100644 (file)
@@ -27,7 +27,6 @@
 #include <deal.II/fe/fe_values.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -189,24 +188,14 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
-
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index c0209fd03adbc1d9caf1e6f9e007d8f474335715..a65652d3042d68b96fa2934b6e570036b9ff8ff3 100644 (file)
@@ -32,23 +32,14 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   if (dim < 3 || fe_degree < 2)
     tria.refine_global(1);
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  cell = tria.begin_active();
-  for (; cell != endc; ++cell)
+  for (const auto &cell : tria.active_cell_iterators())
     if (cell->center().norm() < 1e-8)
       cell->set_refine_flag();
   tria.execute_coarsening_and_refinement();
index f681dfd0e952435bfabee74f444a0101119ecd71..98bee529f180eda9ab448ccd6f8453aa6a8938a1 100644 (file)
@@ -31,7 +31,6 @@
 #include <deal.II/fe/fe_values.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -106,16 +105,8 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   tria.refine_global(4 - dim);
 
   // refine a few cells
index 81d272b191af33bb4546ecd656b00a96cae611d1..9cedcb38262c8748b26d8d31644cd9615dcc5fc7 100644 (file)
@@ -25,7 +25,6 @@
 #include <deal.II/dofs/dof_tools.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 
 #include <deal.II/hp/fe_values.h>
 
@@ -92,14 +91,8 @@ void
 test()
 {
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  for (const auto &cell : tria.active_cell_iterators())
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   tria.refine_global(1);
 
   // refine a few cells
index 7f258fe33ed73f7044f2cae275446b2aaffa7a90..37d90c0713dc9c733806f4f255db4d11edcecd83 100644 (file)
@@ -77,28 +77,20 @@ test()
     return;
 
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   tria.refine_global(1);
 
   // refine a few cells
   for (unsigned int i = 0; i < 11 - 3 * dim; ++i)
     {
-      typename Triangulation<dim>::active_cell_iterator cell =
-                                                          tria.begin_active(),
-                                                        endc = tria.end();
-      unsigned int counter                                   = 0;
-      for (; cell != endc; ++cell, ++counter)
-        if (counter % (7 - i) == 0)
-          cell->set_refine_flag();
+      unsigned int counter = 0;
+      for (const auto &cell : tria.active_cell_iterators())
+        {
+          if (counter % (7 - i) == 0)
+            cell->set_refine_flag();
+          ++counter;
+        }
       tria.execute_coarsening_and_refinement();
     }
 
index 1621dec4ccd9b936da0f785e9c44fd1468cd30cb..dc270a51492dc61c69ca7e8e4fc37228c7610d94 100644 (file)
@@ -29,7 +29,6 @@
 #include <deal.II/fe/fe_values.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/hp/fe_values.h>
@@ -88,14 +87,8 @@ test()
     return;
 
   using number = double;
-  const SphericalManifold<dim>              manifold;
   parallel::distributed::Triangulation<dim> tria(MPI_COMM_WORLD);
   GridGenerator::hyper_ball(tria);
-  for (const auto &cell : tria.active_cell_iterators())
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   tria.refine_global(1);
 
   // refine a few cells
index 35d56c014e105b7e6256da49dff31ef6973a3c9a..f32ca508ae80b5d1fee3bdeccb95382306f631e1 100644 (file)
@@ -34,17 +34,9 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria(
+  Triangulation<dim> tria(
     Triangulation<dim>::limit_level_difference_at_vertices);
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   tria.refine_global(5 - dim);
 
   FE_Q<dim> fe(fe_degree);
index 00487f48171b8646979adca260a70c8d786f6999..5b88acbc9bcef8021535fa36ee15bd7b212d94d5 100644 (file)
@@ -32,7 +32,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -125,11 +124,8 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  SphericalManifold<dim> manifold;
-  Triangulation<dim>     triangulation;
+  Triangulation<dim> triangulation;
   GridGenerator::hyper_shell(triangulation, Point<dim>(), 0.5, 1., 96, true);
-  triangulation.set_all_manifold_ids(0);
-  triangulation.set_manifold(0, manifold);
 
   triangulation.begin_active()->set_refine_flag();
   triangulation.last()->set_refine_flag();
index ea8dac3c6ffdb4152c5528faee83dd1a723884c5..4e1230c477d177f7ac63f35a17fda756f7f1a5cf 100644 (file)
@@ -29,7 +29,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -121,11 +120,8 @@ template <int dim>
 void
 test(const unsigned int fe_degree)
 {
-  SphericalManifold<dim> manifold;
-  Triangulation<dim>     triangulation;
+  Triangulation<dim> triangulation;
   GridGenerator::hyper_shell(triangulation, Point<dim>(), 0.5, 1., 96, true);
-  triangulation.set_all_manifold_ids(0);
-  triangulation.set_manifold(0, manifold);
 
   triangulation.begin_active()->set_refine_flag();
   triangulation.last()->set_refine_flag();
index 0984a9c94b7694b6cd08ffcd1bc28f48d798edc9..295c6c053f35a467c78778fc36f9f11e9863d9e9 100644 (file)
@@ -33,7 +33,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -132,11 +131,8 @@ template <int dim, int fe_degree>
 void
 test(const FESystem<dim> &fe)
 {
-  SphericalManifold<dim> manifold;
-  Triangulation<dim>     triangulation;
+  Triangulation<dim> triangulation;
   GridGenerator::hyper_shell(triangulation, Point<dim>(), 0.5, 1., 96, true);
-  triangulation.set_all_manifold_ids(0);
-  triangulation.set_manifold(0, manifold);
   triangulation.begin_active()->set_refine_flag();
   triangulation.last()->set_refine_flag();
   triangulation.execute_coarsening_and_refinement();
index 7b069025419e21e42d4b0e23fa1cdf418701cefc..3c3aec06321e823da046557ce3619ed1291ce673 100644 (file)
@@ -26,7 +26,6 @@
 #include <deal.II/fe/fe_values.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -179,16 +178,8 @@ template <int dim, int fe_degree>
 void
 test()
 {
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
 
   // refine first and last cell
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
index 0c29cd525c06e6f11a5846024dfcf9f376ab22ca..c5305d8d4c2196c3f9ecbbb97b5ac56176fa64cf 100644 (file)
@@ -26,7 +26,6 @@
 #include <deal.II/fe/mapping_q.h>
 
 #include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/tria.h>
 
 #include <deal.II/lac/affine_constraints.h>
@@ -48,16 +47,8 @@ void
 test()
 {
   using number = double;
-  const SphericalManifold<dim> manifold;
-  Triangulation<dim>           tria;
+  Triangulation<dim> tria;
   GridGenerator::hyper_ball(tria);
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    for (const unsigned int f : GeometryInfo<dim>::face_indices())
-      if (cell->at_boundary(f))
-        cell->face(f)->set_all_manifold_ids(0);
-  tria.set_manifold(0, manifold);
   tria.refine_global(5 - dim);
 
   MappingQ<dim>   mapping(4);

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.