]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Remove HyperShellBoundary from step-32.
authorDavid Wells <wellsd2@rpi.edu>
Thu, 5 Apr 2018 21:40:43 +0000 (17:40 -0400)
committerDavid Wells <wellsd2@rpi.edu>
Fri, 20 Apr 2018 21:09:06 +0000 (17:09 -0400)
As of ff7cc4c8075 this functionality is also in Manifold, so we do not
need to use a Boundary object here.

examples/step-32/step-32.cc

index c39623cefb3f43c7a0a80c53c3b9d23cbd40f851..765075000fba84fb35bbf81e5b3352fac09e511b 100644 (file)
@@ -51,7 +51,6 @@
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/filtered_iterator.h>
 #include <deal.II/grid/manifold_lib.h>
-#include <deal.II/grid/tria_boundary_lib.h>
 #include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/grid_refinement.h>
 
@@ -3580,15 +3579,6 @@ namespace Step32
   // instead of a simple cube geometry), and use the
   // <code>project_temperature_field()</code> function instead of the library
   // function <code>VectorTools::project</code>.
-  // In this example, however, we define both a SphericalManifold() and a
-  // HyperShellBoundary() object to describe the geometry of the domain.
-  // The reason we do so here, is because we want to impose no normal flux
-  // boundary conditions, and they require knowledge of the normals to a boundary,
-  // which a SphericalManifold() alone cannot compute. Consequently, we set
-  // all manifold indicators of cells and adjacent edges to zero, then overwrite
-  // the manifold indicators of all boundary objects by one. We then associate
-  // a SphericalManifold object with the former, and the HyperShellBoundary
-  // object that can also provide normal vectors with the latter.
   template <int dim>
   void BoussinesqFlowProblem<dim>::run ()
   {
@@ -3598,12 +3588,6 @@ namespace Step32
                                 EquationData::R1,
                                 (dim==3) ? 96 : 12,
                                 true);
-    triangulation.set_all_manifold_ids(0);
-    triangulation.set_all_manifold_ids_on_boundary(1);
-    static SphericalManifold<dim> manifold;
-    static HyperShellBoundary<dim> boundary;
-    triangulation.set_manifold (0, manifold);
-    triangulation.set_manifold (1, boundary);
 
     global_Omega_diameter = GridTools::diameter (triangulation);
 

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.