From 729a5a2d6d2da82df2d770be5f5a7f8db342146f Mon Sep 17 00:00:00 2001
From: Abdelkrim Mebarki <Abdelkrim.Mebarki@sophia.inria.fr>
Date: Mon, 27 Feb 2006 23:04:10 +0000
Subject: [PATCH] update

---
 Stream_lines_2/include/CGAL/Stream_lines_2.h | 45 +++++++++++++++++++-
 1 file changed, 43 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 af4ee5a7fee..d2434e50223 100644
--- a/Stream_lines_2/include/CGAL/Stream_lines_2.h
+++ b/Stream_lines_2/include/CGAL/Stream_lines_2.h
@@ -96,7 +96,43 @@ protected:
   FT                                                             separating_distance;
   FT                                                             saturation_ratio;
   Point_2                                                        seed_point;
-  unsigned int                                                   _number_of_lines;
+  int                                                            samp_step;
+  unsigned int _number_of_lines;
+  Vector_field_2 * vf_2;
+  Integrator_2 * int_2;
+public:
+  void set_separating_distance(FT new_value){separating_distance = new_value;}
+  void set_saturation_ratio(FT new_value){ saturation_ratio = new_value;}
+  void update()
+    {
+      ir = il = 0; // initialization
+      fSepStl_seed = separating_distance*saturation_ratio;
+      m_DT.clear();
+      Point_2 pPoint;
+      pPoint = Point_2(min_x-separating_distance,min_y-separating_distance);
+      m_DT.insert(pPoint);
+      pPoint = Point_2(min_x-separating_distance,max_y+separating_distance);
+      m_DT.insert(pPoint);
+      pPoint = Point_2(max_x+separating_distance,min_y-separating_distance);
+      m_DT.insert(pPoint);
+      pPoint = Point_2(max_x+separating_distance,max_y+separating_distance);
+      m_DT.insert(pPoint);
+      for (int i=(int) (min_x-separating_distance);i<max_x+(int)
+	     separating_distance;i=i+(int) (fSepStl_seed))
+	{
+	  pPoint = Point_2((FT)i,(FT)max_y+separating_distance);
+	  m_DT.insert(pPoint);
+	  pPoint = Point_2((FT)max_x+separating_distance,(FT)i);
+	  m_DT.insert(pPoint);
+	  pPoint = Point_2((FT)i,min_y-separating_distance);
+	  m_DT.insert(pPoint);
+	  pPoint = Point_2(min_x-separating_distance,(FT)i);
+	  m_DT.insert(pPoint);
+	}
+      _number_of_lines = 0;
+      place_stream_lines(*vf_2, *int_2,
+			 samp_step);     
+    }
 protected:
   void place_stream_lines(const Vector_field_2 & vector_field_2, const Integrator_2 & integrator,
 			  const int & sampling_step, const bool & step_by_step = false);
@@ -196,6 +232,10 @@ vector_field_2, const Integrator_2 & m_integrator, const FT & m_separating_dista
       m_DT.insert(pPoint);
     }
   _number_of_lines = 0;
+  vf_2 = new Vector_field_2(vector_field_2);
+  int_2 = new Integrator_2(m_integrator);
+  samp_step = sampling_step;
+  stl_container.clear();
   place_stream_lines(vector_field_2, m_integrator,
 		     sampling_step, step_by_step);
 }
@@ -761,7 +801,8 @@ Stream_lines_2<VectorField_2, Integrator_2>::end()
 template <class VectorField_2, class Integrator_2>
 inline
 void Stream_lines_2<VectorField_2, Integrator_2>::make_iterator()
-{ 
+{
+  iterator_container.clear();
   for(typename Stream_line_container_2::iterator
 	begin=stl_container.begin();
       begin!=stl_container.end();begin++)