Fixed inaccuracy of j0f (BZ #28185)
commit6bbf7298323bf31bc43494b2201465a449778e10
authorPaul Zimmermann <Paul.Zimmermann@inria.fr>
Tue, 5 Oct 2021 08:32:36 +0000 (5 10:32 +0200)
committerPaul Zimmermann <Paul.Zimmermann@inria.fr>
Tue, 5 Oct 2021 11:45:37 +0000 (5 13:45 +0200)
tree4dabd416c1e18864b18b09549d494f101df36420
parenta312e8fe6d89f5eae6a4583d5db577121e61c0b5
Fixed inaccuracy of j0f (BZ #28185)

The largest errors over the full binary32 range are after this
patch (on x86_64):

RNDN: libm wrong by up to 9.00e+00 ulp(s) [9] for x=0x1.04c39cp+6
RNDZ: libm wrong by up to 9.00e+00 ulp(s) [9] for x=0x1.04c39cp+6
RNDU: libm wrong by up to 9.00e+00 ulp(s) [9] for x=0x1.04c39cp+6
RNDD: libm wrong by up to 8.98e+00 ulp(s) [9] for x=0x1.4b7066p+7

Inputs that were yielding huge errors have been added to "make check".
Reviewed-by: Adhemeral Zanella <adhemerval.zanella@linaro.org>
math/auto-libm-test-in
math/auto-libm-test-out-j0
sysdeps/ieee754/flt-32/e_j0f.c