From: David Wells <wellsd2@rpi.edu>
Date: Fri, 27 Apr 2018 17:39:37 +0000 (-0400)
Subject: Remove a static manifold.
X-Git-Tag: v9.0.0-rc1~88^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F6364%2Fhead;p=dealii.git

Remove a static manifold.
---

diff --git a/examples/step-5/step-5.cc b/examples/step-5/step-5.cc
index 9c57640eda..ef62241bc9 100644
--- a/examples/step-5/step-5.cc
+++ b/examples/step-5/step-5.cc
@@ -441,12 +441,14 @@ void Step5<dim>::run ()
   // <code>GridIn</code> class to find out what input formats are presently
   // supported.)
 
-  // The grid in the file describes a circle. Therefore we have to use
-  // a manifold object which tells the triangulation where to put new
-  // points on the boundary when the grid is refined. This works in
-  // the same way as in the first example, but in this case we only
-  // set the manifold ids of the boundary.
-  static const SphericalManifold<dim> boundary;
+  // The grid in the file describes a circle. Therefore we have to use a
+  // manifold object which tells the triangulation where to put new points on
+  // the boundary when the grid is refined. Unlike step-1, since GridIn does
+  // not know that the domain has a circular boundary (unlike
+  // GridGenerator::hyper_shell) we have to explicitly attach a manifold to
+  // the boundary after creating the triangulation to get the correct result
+  // when we refine the mesh.
+  const SphericalManifold<dim> boundary;
   triangulation.set_all_manifold_ids_on_boundary(0);
   triangulation.set_manifold (0, boundary);