From 474e221c6ed91ba9715ce6ba826c5422f364d47b Mon Sep 17 00:00:00 2001 From: Peter Munch 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::active_cell_iterator, Point> 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::active_cell_iterator, -- 2.39.5