From ff079642b66eddb9f06bd306585657150e1b23b0 Mon Sep 17 00:00:00 2001
From: Andreas Fabri <andreas.fabri@geometryfactory.com>
Date: Fri, 6 Dec 2024 14:15:05 +0000
Subject: [PATCH] Add two curves to dD examples (we get an assertion)

---
 .../examples/Frechet_distance/Frechet_distance_d.cpp     | 9 ++++++---
 Kernel_23/include/CGAL/Bbox.h                            | 7 +++++++
 2 files changed, 13 insertions(+), 3 deletions(-)

diff --git a/Frechet_distance/examples/Frechet_distance/Frechet_distance_d.cpp b/Frechet_distance/examples/Frechet_distance/Frechet_distance_d.cpp
index 8acde5a1db2..1726b42c21d 100644
--- a/Frechet_distance/examples/Frechet_distance/Frechet_distance_d.cpp
+++ b/Frechet_distance/examples/Frechet_distance/Frechet_distance_d.cpp
@@ -10,9 +10,12 @@ using Kernel = CGAL::Epick_d<CGAL::Dimension_tag<4>>;
 using Traits = CGAL::Frechet_distance_traits_d<Kernel>;
 using Point = Traits::Point_d;
 
-int main(int argc, char* argv[])
+int main(int , char*)
 {
-    std::vector<Point> A, B;
-    bool res = CGAL::is_Frechet_distance_larger<Traits>(A, B, 0.001);
+    std::array<Point,4> A = { Point(0,0,0,0), Point(1,0,0,0), Point(1,1,0,1),Point(1,1,1,0)};
+    std::array<Point,4> B = { Point(0,0,0,0), Point(1,0,0,1), Point(1,1,0,0),Point(1,1,1,1)};
+
+    std::pair<double, double> res = CGAL::approximate_Frechet_distance<Traits>(A, B, 0.000001);
+    std::cout << "The Frechet distance between the polylines is between " <<  res.first << " and " << res.second << std::endl;
     return 0;
 }
diff --git a/Kernel_23/include/CGAL/Bbox.h b/Kernel_23/include/CGAL/Bbox.h
index b0fd9df0647..9210af5441f 100644
--- a/Kernel_23/include/CGAL/Bbox.h
+++ b/Kernel_23/include/CGAL/Bbox.h
@@ -175,6 +175,13 @@ std::istream& operator>>(std::istream& in, Impl::Bbox<Container, Derived>& bbox)
     return in;
 }
 
+template<int N, typename T>
+Bbox<Dimension_tag<N>, T> operator+(Bbox<Dimension_tag<N>, T> bbox, const Bbox<Dimension_tag<N>, T>& other)
+{
+    bbox += other;
+    return bbox;
+}
+
 
 } // namespace CGAL