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

CGAL · Dec 9, 2024 · ddbfcd0 · ddbfcd0
1 parent 55a3cb5
commit ddbfcd0
Showing 1 changed file with 5 additions and 0 deletions.
diff --git 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;