[Boost-users] math tools roots

Subject: [Boost-users] math tools roots
From: Matwey V. Kornilov (matwey.kornilov_at_[hidden])
Date: 2011-10-09 09:47:54

Hi all,

I believe that there is a bug in eps_tolerance implementation. At line 34 of 
toms748_solve.hpp one can see:

return (fabs(a - b) / (std::min)(fabs(a), fabs(b))) <= eps;

When both a and b equal zero this condition is evaluated to a false instead 
of a true. Left hand of condition is NaN and the right one is a floating-
point number hence the result is false.

Well obviously, a user wants to see true here because the relative length of 
interval in this case tend to zero.