Add assert to compute_point_locations_try_all

diff --git a/include/deal.II/grid/grid_tools.h b/include/deal.II/grid/grid_tools.h
index 3ce447ad3657933484ecfde0389dd60819b869bc..939f93c31668b14808273b870a615758991fa5ed 100644 (file)
--- a/include/deal.II/grid/grid_tools.h
+++ b/include/deal.II/grid/grid_tools.h
@@ -889,6 +889,8 @@ namespace GridTools
    * The algorithm builds an rtree of @p points to sort them spatially, before
    * attempting to call find_active_cell_around_point().
    *
+   * @note This function is not implemented for the codimension one case (<tt>spacedim != dim</tt>).
+   *
    * @note If a point is not found inside the mesh, or is lying inside an
    * artificial cell of a parallel::TriangulationBase, the point is silently
    * ignored. If you want to infer for which points the search failed, use the
@@ -951,6 +953,8 @@ namespace GridTools
    *   >
    * @endcode
    *
+   * @note This function is not implemented for the codimension one case (<tt>spacedim != dim</tt>).
+   *
    * @note This function optimizes the search by making use of
    * GridTools::Cache::get_cell_bounding_boxes_rtree(), which either returns
    * a cached rtree or builds and stores one. Building an rtree might hinder

diff --git a/source/grid/grid_tools.cc b/source/grid/grid_tools.cc
index 94abc37f78ab354ea8a7b1f4899bbb3b94553f5a..c557384fc05e8f18c64db3406033ba1bd586e01d 100644 (file)
--- a/source/grid/grid_tools.cc
+++ b/source/grid/grid_tools.cc
@@ -5457,6 +5457,9 @@ namespace GridTools
     const typename Triangulation<dim, spacedim>::active_cell_iterator
       &cell_hint)
   {
+    Assert((dim == spacedim),
+           ExcMessage("Only implemented for dim==spacedim."));
+
     // Alias
     namespace bgi = boost::geometry::index;