From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Sat, 5 Aug 2023 10:25:46 +0000 (-0600)
Subject: Make sure we don't match points on the perimeter of the domain as 'inside'.
X-Git-Tag: relicensing~617^2~2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a3fe6f2fdb2d5c67b2c08c852452aa02f180e46c;p=dealii.git

Make sure we don't match points on the perimeter of the domain as 'inside'.
---

diff --git a/source/grid/grid_generator.cc b/source/grid/grid_generator.cc
index de4ce69097..6c25951c6a 100644
--- a/source/grid/grid_generator.cc
+++ b/source/grid/grid_generator.cc
@@ -3450,7 +3450,7 @@ namespace GridGenerator
                 if (face->at_boundary() &&
                     internal::point_in_2d_box(face->center(),
                                               center,
-                                              outer_radius))
+                                              outer_radius * (1. - 1e-12)))
                   face->set_manifold_id(polar_manifold_id);
                 else
                   face->set_manifold_id(tfi_manifold_id);