From: Luca Heltai Date: Thu, 14 May 2020 20:57:21 +0000 (+0200) Subject: Dirty fix. X-Git-Tag: v9.2.0-rc2~4^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F10218%2Fhead;p=dealii.git Dirty fix. --- diff --git a/source/grid/grid_tools.cc b/source/grid/grid_tools.cc index 5b41068c71..ce8fb651a4 100644 --- a/source/grid/grid_tools.cc +++ b/source/grid/grid_tools.cc @@ -1724,9 +1724,13 @@ namespace GridTools // the user did not provide a hint_cell, and at the beginning of the // while loop we performed an actual global search on the mesh // vertices. Not finding the point then means the point is outside the - // domain. - AssertThrow(current_cell.state() == IteratorState::valid, - ExcPointNotFound(p)); + // domain, or that we've had problems with the algorithm above. Try as a + // last resort the other (simpler) algorithm. + if (current_cell.state() != IteratorState::valid) + return find_active_cell_around_point(mapping, + mesh, + p, + marked_vertices); current_cell = typename MeshType::active_cell_iterator(); }