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

fix corner case in round, where the float is representing 

a power of 2 in the mantissa
This commit is contained in:
Michael Hemmer 2010-05-11 11:41:12 +00:00
parent 7d57644b59
commit 1ab0e40c06
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Number_types/include/CGAL/CORE_BigFloat.h
View File

@ -279,11 +279,12 @@ round(const CORE::BigFloat& x, long rel_prec = CORE::defRelPrec.toLong() ){
// else
// std::cout <<"current prec: " << " SINGLETON " << std::endl;
// std::cout <<"desired prec: " << rel_prec << std::endl;
// std::cout <<"bitLength: " << CORE::bitLength(m) << std::endl;
// long shift = ::CORE::bitLength(m) - rel_prec - 1;
long shift ;
if (err == 0)
shift = ::CORE::bitLength(m) - rel_prec - 2;
shift = ::CORE::bitLength(m) - rel_prec - 3;
else
shift = CGAL::relative_precision(x) - rel_prec -1;
@ -330,8 +331,9 @@ public:
CGAL_precondition(!Singleton()(x));
CGAL_precondition(!CGAL::zero_in(x));
x = x.abs();
NT w = Width()(x);
w /= ::CORE::BigFloat(x.m()-x.err(),0,x.exp());
w /= ::CORE::BigFloat(x.m()-x.err(),0,x.exp());
w = w.abs();
return -(CORE::ceilLg(w.m()+w.err())+w.exp()*14);
}