From a3fe6f2fdb2d5c67b2c08c852452aa02f180e46c Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Sat, 5 Aug 2023 04:25:46 -0600
Subject: [PATCH] Make sure we don't match points on the perimeter of the
 domain as 'inside'.

---
 source/grid/grid_generator.cc | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

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);
-- 
2.39.5