Re: issue with tgammaf

Paul Zimmermann <[email protected]> writes:

>> You might have more luck with lgammaf_r as that avoids use of
>> _impure_ptr in newlib.
>
> thank you very much Keith, that made it. The results for lgammaf_r are
> the same as those with OpenLibm and Musl:

That isn't surprising; so many libraries math support can be traced back
to the original SunPro code.

Thanks much for doing this analysis; it's nice to see careful
measurements of how accurate each of these C libraries is. Have you
published the test framework so we can replicate these results and
monitor their change over time?

Do you know that newlib has two different sets of 32-bit math functions?
Arm supplied a new set of many common math functions that use 64-bit
doubles for the intermediate computations? I'd be very interested in
measuring the effect of that code on accuracy. The new paths are
selected by compiling newlib with -D__OBSOLETE_MATH_DEFAULT=0 while the
old can be selected with -D__OBSOLETE_MATH_DEFAULT=1

-- 
-keith

-----BEGIN PGP SIGNATURE-----

iQIzBAEBCAAdFiEEw4O3eCVWE9/bQJ2R2yIaaQAAABEFAl9jiscACgkQ2yIaaQAA
ABGGDw//RneKzvPM2J0PcKOwDxFa/4LxrPxb805qWeEcwP1rTeRObpQh+wR5bOzC
kpPaBiA24ujHwD2eZ1n7LB3CnrXFNK5ibHWORJry9Os/d4PlwlqLfJKh6daZyISb
5KlbEajledA9G39MaSQJVHCka+hqKvzaXu0is/di+pjtUlKJGmNo0ARafrRjUvTh
6tZhiBVuhoyqvzrGOaQWtep09Z72h10SvPP4xF2ua3NrgppcwRyA5mbN0aUvfWzb
MTZR4fxKaBxobMtbxrAC1nJIyp0pIoPayc5Q69ShhbsO3wIdwcltZpYhjvtcPsZA
2DrunGR8NIt+np4LZfpFAhL9SUsBqZeT4+oyp6K29VcDgOt9PSxv9R47irM/Srf2
2ewuMnLaWGRRW+QWb97Zts27rZ6iZTGm7IgI5t0HHKIOEHa6A1rPTqLCUX8255l4
kEq+gsvG+rweWFuwE2+2XLwR8MPtJs5fz6MfnF13hHcWZXA6kqcr46kBW8KCLKUp
KCLFn41CRvUJGj9x+8hctW1uYrgGdqmYBVSUlX3Id6U4n252hwNvYDCxyOmbGXfX
q7BgAISPnwuvPF7HwCIcmLyZVcgdP1hPxUDZr/avcbI1duG8yhSxTdPJNCvmkxxl
c9i0iSFpISTL4p++YccfJsNRNLORWA67M97+V4nK9QaEqKKEk9A=
=QIpv
-----END PGP SIGNATURE-----