songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Fix bug : check the length of the new segment inserted in the streamline

This commit is contained in:
Abdelkrim Mebarki 2007-11-15 16:13:32 +00:00
parent 7b314c3825
commit 6d6970d681
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Stream_lines_2/include/CGAL/Stream_lines_2.h
View File

@ -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);