[PATCH] Fix the inaccuracy of j0f (BZ 14469) and y0f (BZ 14471) [v2]
Paul Zimmermann
Paul.Zimmermann@inria.fr
Mon Jan 25 09:59:14 GMT 2021
Dear Adhemerval,
> So if I understand correctly, you patch increase the ulps for non-default
> rounding for this specific inputs? If it were the case, could be fix it
> as well?
I guess the previous code already had an error > 9 ulps on this specific
input:
Failure: Test: j0_downward (0x3.17bcfcp+4)
Result:
is: 1.1670220923447945e-04 0x1.e97c0b0296dc8p-14
should be: 1.1670220923447976e-04 0x1.e97c0b0296ddfp-14
difference: 3.1170812458958252e-19 0x1.7000000000000p-62
ulp : 23.0000
max.ulp : 2.0000
However on my machine I have an error of only 9 ulps for that input.
Maybe your machine doesn't have fma? Maybe I should develop my code
on a machine without fma, assuming the errors will be smaller with fma.
Paul
More information about the Libc-alpha
mailing list