]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Make sure we don't match points on the perimeter of the domain as 'inside'.
authorWolfgang Bangerth <bangerth@colostate.edu>
Sat, 5 Aug 2023 10:25:46 +0000 (04:25 -0600)
committerWolfgang Bangerth <bangerth@colostate.edu>
Sat, 5 Aug 2023 10:25:46 +0000 (04:25 -0600)
source/grid/grid_generator.cc

index de4ce69097473666723f3ca4e8aad97905ab1d6f..6c25951c6a46fa40054e1d794f093dd1408adcd5 100644 (file)
@@ -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);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.