[PATCH 09/17] math: Use exp10m1f from CORE-MATH

Florian Weimer fweimer@redhat.com
Wed Nov 6 10:29:04 GMT 2024


* Adhemerval Zanella:

> The CORE-MATH implementation is correctly rounded (for any rounding mode)
> and shows better performance compared to the generic exp10m1f.

When compiled with GCC 11, it's not correctly rounded on x86-64 (both
current AMD and Intel CPUs, with and without FMA):

Failure: Test: exp10m1_towardzero (-0x1.1p+4)                                   
Result:                                                                         
 is:         -1.00000000e+00  -0x1.000000p+0                                    
 should be:  -9.99999940e-01  -0x1.fffffep-1                                    
 difference:  5.96046447e-08   0x1.000000p-24                                   
 ulp       :  1.0000                                                            
 max.ulp   :  0.0000                                                            

I don't see this with GCC 14.

Any suggestions how to debug this further?

Thanks,
Florian



More information about the Libc-alpha mailing list