https://gcc.gnu.org/bugzilla/show_bug.cgi?id=63488 --- Comment #7 from Paul Zimmermann <zimmerma+gcc at loria dot fr> --- I agree that near zeroes we can expect large errors. However for other functions I got only small errors in ulps, maybe I was unlucky. Also the ultimate goal is to get correct rounding, even near zeroes.