From 474e221c6ed91ba9715ce6ba826c5422f364d47b Mon Sep 17 00:00:00 2001
From: Peter Munch <peterrmuench@gmail.com>
Date: Sun, 25 Jul 2021 15:00:47 +0200
Subject: [PATCH] Let find_active_cell_around_point() return mesh.end() if
 point is not found in BBs

---
 source/grid/grid_tools.cc | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/source/grid/grid_tools.cc b/source/grid/grid_tools.cc
index 5c9b766527..b2d78e2a70 100644
--- a/source/grid/grid_tools.cc
+++ b/source/grid/grid_tools.cc
@@ -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,
-- 
2.39.5