From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Wed, 1 May 2024 15:38:17 +0000 (+0530)
Subject: Ensure we only set boundary ids on boundary faces.
X-Git-Tag: v9.6.0-rc1~326^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=2a52a6a0d33fc2d2aa6c4493ac50a4cfb8f7c5ad;p=dealii.git

Ensure we only set boundary ids on boundary faces.
---

diff --git a/examples/step-3/doc/results.dox b/examples/step-3/doc/results.dox
index d5baf27f87..7aade07768 100644
--- a/examples/step-3/doc/results.dox
+++ b/examples/step-3/doc/results.dox
@@ -122,13 +122,14 @@ suggestions:
   For example, we can label all of the cells along the top and
   bottom boundaries with a boundary indicator 1 by checking to
   see if the cell centers' y-coordinates are within a tolerance
-  (here 1e-12) of -1 and 1. Try this immediately after calling
+  (here `1e-12`) of -1 and 1. Try this immediately after calling
   GridGenerator::hyper_cube(), as before:
   @code
   for (auto &face : triangulation.active_face_iterators())
-    if (std::fabs(face->center()(1) - (-1.0)) < 1e-12 ||
-        std::fabs(face->center()(1) - (1.0)) < 1e-12)
-      face->set_boundary_id(1);
+    if (face->at_boundary())
+      if (std::fabs(face->center()(1) - (-1.0)) < 1e-12 ||
+          std::fabs(face->center()(1) - (1.0)) < 1e-12)
+        face->set_boundary_id(1);
   @endcode
   Although this code is a bit longer than before, it is useful for
   complex geometries, as it does not require knowledge of face labels.