From a3fe6f2fdb2d5c67b2c08c852452aa02f180e46c Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth 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