From: Luca Heltai <luca.heltai@sissa.it>
Date: Thu, 14 May 2020 20:57:21 +0000 (+0200)
Subject: Dirty fix.
X-Git-Tag: v9.3.0-rc1~1630^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=82a2e4b0d74dc64d574c442216c8f95822823fda;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<spacedim>(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<dim, spacedim>::active_cell_iterator();
       }