Fix for two identical points. Must be checked with André

2024-12-09 14:51:42 +00:00 · 2024-12-09 14:51:42 +00:00 · 4671ed9ebd
parent 016dfbfe07
commit 4671ed9ebd
1 changed files with 5 additions and 0 deletions
--- a/Frechet_distance/include/CGAL/Frechet_distance/internal/high_level_predicates.h
+++ b/Frechet_distance/include/CGAL/Frechet_distance/internal/high_level_predicates.h
@ -76,6 +76,11 @@ fill_lambda(const Point& circle_center,
  }
  c -= CGAL::square(typename Traits::FT(radius));

+  if(is_zero(a)){
+    I = std::make_pair(Lambda<C>(0), Lambda<C>(1));
+    return true;
+  }
+
  FT minus_b_div_a = b / a;
  FT d = CGAL::square(minus_b_div_a) - c / a;