]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Let find_active_cell_around_point() return mesh.end() if point is not found in BBs 12599/head
authorPeter Munch <peterrmuench@gmail.com>
Sun, 25 Jul 2021 13:00:47 +0000 (15:00 +0200)
committerPeter Munch <peterrmuench@gmail.com>
Sun, 25 Jul 2021 13:00:47 +0000 (15:00 +0200)
source/grid/grid_tools.cc

index 5c9b76652778317dbbe95a3c47029087d0400878..b2d78e2a70627f878554e0d22cd94306ed9c6d47 100644 (file)
@@ -2774,6 +2774,8 @@ namespace GridTools
     std::pair<typename MeshType<dim, spacedim>::active_cell_iterator,
               Point<dim>>
       cell_and_position;
+    cell_and_position.first = mesh.end();
+
     // To handle points at the border we keep track of points which are close to
     // the unit cell:
     std::pair<typename MeshType<dim, spacedim>::active_cell_iterator,

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.