From a1e46031c0f5d36b7de8d7a2d170d16355b4ea64 Mon Sep 17 00:00:00 2001 From: David Wells Date: Fri, 27 Apr 2018 13:39:37 -0400 Subject: [PATCH] Remove a static manifold. --- examples/step-5/step-5.cc | 14 ++++++++------ 1 file changed, 8 insertions(+), 6 deletions(-) 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); -- 2.39.5