From 6d6970d68168ded8c55dab2ccf133e96e6421d8a Mon Sep 17 00:00:00 2001
From: Abdelkrim Mebarki <Abdelkrim.Mebarki@sophia.inria.fr>
Date: Thu, 15 Nov 2007 16:13:32 +0000
Subject: [PATCH] Fix bug : check the length of the new segment inserted in the
 streamline

---
 Stream_lines_2/include/CGAL/Stream_lines_2.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/Stream_lines_2/include/CGAL/Stream_lines_2.h b/Stream_lines_2/include/CGAL/Stream_lines_2.h
index 9d76e080874..30fe0fb360e 100644
--- a/Stream_lines_2/include/CGAL/Stream_lines_2.h
+++ b/Stream_lines_2/include/CGAL/Stream_lines_2.h
@@ -386,7 +386,7 @@ Stream_lines_2<VectorField_2, Integrator_2>::integrate_forward(const Vector_fiel
       if(number_of_points > 30)
   bEnd = bEnd || ((distance(stl.front(), stl.back()))<vector_field_2.get_integration_step());
       FT dist_ = distance(ex_old_point,new_point);
-      bEnd = bEnd || dist_ <= vector_field_2.get_integration_step();
+      bEnd = bEnd || dist_ < 0.9*vector_field_2.get_integration_step();
       if (!bEnd)
   {
     if(sampling != sampling_step)
@@ -512,7 +512,7 @@ void Stream_lines_2<VectorField_2, Integrator_2>::integrate_backward(const Vecto
       new_point = integrator(old_point,vector_field_2,false);
       bEnd = !vector_field_2.is_in_domain(new_point);
       FT dist_ = distance(ex_old_point,new_point);
-      bEnd = bEnd || dist_ <= vector_field_2.get_integration_step() || (new_point == old_point);/* to review */ 
+      bEnd = bEnd || dist_ < 0.9*vector_field_2.get_integration_step() || (new_point == old_point);/* to review */ 
       if(number_of_points > 30)
   bEnd = bEnd || ((distance(stl.front(), stl.back()))<vector_field_2.get_integration_step());
       //    bEnd = bEnd || (number_of_points > 3000);