correctly rounded mathematical functions

Paul Zimmermann Paul.Zimmermann@inria.fr
Thu Jan 6 15:33:11 GMT 2022


       Dear Joel,

> > as far as I know, four libraries provide long double functions:
> > GNU libc, Intel Math Library, OpenLibm and Musl. See Table 11 from
> > https://members.loria.fr/PZimmermann/papers/accuracy.pdf.
> 
> Is the FreeBSD long double support related to any of those libraries?
> 
> https://github.com/freebsd/freebsd-src/tree/main/lib/msun

after a quick look at the cosl code
(https://github.com/freebsd/freebsd-src/blob/main/lib/msun/ld80/k_cosl.c)
it seems to be the same as the OpenLibm implementation
(openlibm-0.7.5/ld80/k_cosl.c).

> Any comments on their implementation?

the OpenLibm implementation has huge errors (up to 533 ulps for powl),
but they have fixed some of the issues I have reported to them.

See also https://github.com/JuliaMath/openlibm/issues/236.

Best regards,
Paul


More information about the Newlib mailing list