From: David Wells 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=a1e46031c0f5d36b7de8d7a2d170d16355b4ea64;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::run () // GridIn 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 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 boundary; triangulation.set_all_manifold_ids_on_boundary(0); triangulation.set_manifold (0, boundary);