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

Better solution 

There was still 'double' variables that were not forced to 64 bits. I
opt for another solution: the result of myrand() itself is forced to
64 bits.
This commit is contained in:
Laurent Rineau 2014-07-17 14:32:10 +02:00
parent 4689ac4db5
commit d2bd75d8a6
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
Triangulation_3/test/Triangulation_3/test_static_filters.cpp
View File

@ -70,16 +70,16 @@ double rand_base()
double my_rand()
{
// Ensure 53 random bits, not 48.
return rand_base() + rand_base()/1024;
return CGAL::IA_force_to_double(rand_base() + rand_base()/1024);
}
// Random point in unit cube.
Weighted_point_3 my_rand_wp3()
{
double x = CGAL::IA_force_to_double(my_rand());
double y = CGAL::IA_force_to_double(my_rand());
double z = CGAL::IA_force_to_double(my_rand());
double r = CGAL::IA_force_to_double(my_rand());
double x = my_rand();
double y = my_rand();
double z = my_rand();
double r = my_rand();
return Weighted_point_3( FTr_with_SF::Weighted_point_3(FTr_with_SF::Bare_point(x, y, z),r) , FTr_without_SF::Weighted_point_3(FTr_without_SF::Bare_point(x, y, z),r) );
}
@ -89,7 +89,7 @@ void perturb(Weighted_point_3 &p, double rel_eps)
double x = p.first.x()*(1+rand_base()*rel_eps);
double y = p.first.y()*(1+rand_base()*rel_eps);
double z = p.first.z()*(1+rand_base()*rel_eps);
double r= p.first.weight()*(1+rand_base()*rel_eps);
double r = p.first.weight()*(1+rand_base()*rel_eps);
p=Weighted_point_3( FTr_with_SF::Weighted_point_3(FTr_with_SF::Bare_point(x, y, z),r) , FTr_without_SF::Weighted_point_3(FTr_without_SF::Bare_point(x, y, z),r) );
}