9 Oct
2011
9 Oct
'11
5:21 p.m.
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.
Will fix: can you please file a bug report at svn.boost.org so I don't forget? Thanks, John.