Fix creating of boost point

diff --git a/source/grid/grid_tools.cc b/source/grid/grid_tools.cc
index fcbd2772d1bb05240f3275d7ee60a439efe1ba73..990ada7cacaf00b9757d88bba874cb786f60c5b4 100644 (file)
--- a/source/grid/grid_tools.cc
+++ b/source/grid/grid_tools.cc
@@ -2746,6 +2746,29 @@ namespace GridTools
       // return if the scalar product of a is larger.
       return (scalar_product_a > scalar_product_b);
     }
+
+
+
+    boost::geometry::model::point<double, 1, boost::geometry::cs::cartesian>
+    create_boost_point(const Point<1> &p)
+    {
+      return boost::geometry::model::
+        point<double, 1, boost::geometry::cs::cartesian>(p[0]);
+    }
+
+    boost::geometry::model::point<double, 2, boost::geometry::cs::cartesian>
+    create_boost_point(const Point<2> &p)
+    {
+      return {p[0], p[1]};
+    }
+
+    boost::geometry::model::point<double, 3, boost::geometry::cs::cartesian>
+    create_boost_point(const Point<3> &p)
+    {
+      return {p[0], p[1], p[2]};
+    }
+
+
   } // namespace internal
 
   template <int dim, template <int, int> class MeshType, int spacedim>
@@ -2785,15 +2808,9 @@ namespace GridTools
     if (relevant_cell_bounding_boxes_rtree != nullptr &&
         !relevant_cell_bounding_boxes_rtree->empty())
       {
-        using point_t = boost::geometry::model::
-          point<double, spacedim, boost::geometry::cs::cartesian>;
-        const auto pt =
-          spacedim == 1 ?
-            point_t(p[0]) :
-            (spacedim == 2 ? point_t(p[0], p[1]) : point_t(p[0], p[1], p[2]));
         if (relevant_cell_bounding_boxes_rtree->qbegin(
-              boost::geometry::index::intersects(pt)) ==
-            relevant_cell_bounding_boxes_rtree->qend())
+              boost::geometry::index::intersects(internal::create_boost_point(
+                p))) == relevant_cell_bounding_boxes_rtree->qend())
           return cell_and_position;
       }