Re: [libm] Implementations for rsqrt[fl].
Steve Kargl <[email protected]>
| Newsgroups | gmane.os.freebsd.devel.hackers |
|---|---|
| Message-ID | <[email protected]> |
On 5/8/26 09:38, Steve Kargl wrote: > On 5/7/26 22:55, Paul Zimmermann wrote: > >>> Max ulp: 0.500000 at 1.17549421e-38 >>> %./tlibm rsqrt -fa 1.17549421e-38 >>> x = 1.17549421e-38f, /* 0x007fffff */ >>> libm = 9.22337204e+18f, /* 0x5f000000 */ >>> mpfr = 9.22337314e+18f, /* 0x5f000001 */ >>> ULP = 0.50000 >> >> it seems libm is not correctly rounded in this example. >> Do you plan to fix that? > > If I knew how, yes. :-). > I found a/the fix. -- steve