mirror of
git://sourceware.org/git/glibc.git
synced 2025-02-05 12:40:55 +08:00
e4e52ff059
Bug 21094 reports 3ulp errors of cosf and tanf for certain arguments near pi/2 arising from the use of an insufficiently accurate range reduction. (To be clear, this is a quality-of-implementation issue relating to the apparent intent of those particular cosf and tanf implementations; 3ulp is within the general glibc accuracy goals, so not inherently a bug.) This patch fixes that error by making a wider range of cases use the existing more accurate range reduction for arguments close to pi/2. The wider range of values is still narrow enough for the "z -= pio2_2;" in the more accurate case to be exact, as the code expects. Tested for x86_64, x86 and mips64; no ulps updates needed (but at least on mips64, the larger ulps were seen if the tests were added without the substantive fix). [BZ #21094] * sysdeps/ieee754/flt-32/e_rem_pio2f.c (__ieee754_rem_pio2f): Use 24+24+24-bit pi for wider range of values around pi/2. * math/auto-libm-test-in: Add more tests of cos and tan. * math/auto-libm-test-out-cos: Regenerated. * math/auto-libm-test-out-tan: Likewise. |
||
|---|---|---|
| .. | ||
| e_acosf.c | ||
| e_acoshf.c | ||
| e_asinf.c | ||
| e_atan2f.c | ||
| e_atanhf.c | ||
| e_coshf.c | ||
| e_exp2f.c | ||
| e_expf.c | ||
| e_fmodf.c | ||
| e_gammaf_r.c | ||
| e_hypotf.c | ||
| e_ilogbf.c | ||
| e_j0f.c | ||
| e_j1f.c | ||
| e_jnf.c | ||
| e_lgammaf_r.c | ||
| e_log2f.c | ||
| e_log10f.c | ||
| e_logf.c | ||
| e_powf.c | ||
| e_rem_pio2f.c | ||
| e_remainderf.c | ||
| e_sinhf.c | ||
| e_sqrtf.c | ||
| k_cosf.c | ||
| k_rem_pio2f.c | ||
| k_sinf.c | ||
| k_tanf.c | ||
| lgamma_negf.c | ||
| lgamma_productf.c | ||
| mpn2flt.c | ||
| s_asinhf.c | ||
| s_atanf.c | ||
| s_cbrtf.c | ||
| s_ceilf.c | ||
| s_copysignf.c | ||
| s_cosf.c | ||
| s_erff.c | ||
| s_expm1f.c | ||
| s_fabsf.c | ||
| s_finitef.c | ||
| s_floorf.c | ||
| s_fpclassifyf.c | ||
| s_frexpf.c | ||
| s_fromfpf_main.c | ||
| s_fromfpf.c | ||
| s_fromfpxf.c | ||
| s_getpayloadf.c | ||
| s_isinff.c | ||
| s_isnanf.c | ||
| s_issignalingf.c | ||
| s_llrintf.c | ||
| s_llroundf.c | ||
| s_log1pf.c | ||
| s_logbf.c | ||
| s_lrintf.c | ||
| s_lroundf.c | ||
| s_modff.c | ||
| s_nearbyintf.c | ||
| s_nextafterf.c | ||
| s_nextupf.c | ||
| s_remquof.c | ||
| s_rintf.c | ||
| s_roundevenf.c | ||
| s_roundf.c | ||
| s_scalblnf.c | ||
| s_scalbnf.c | ||
| s_setpayloadf_main.c | ||
| s_setpayloadf.c | ||
| s_setpayloadsigf.c | ||
| s_signbitf.c | ||
| s_sincosf.c | ||
| s_sinf.c | ||
| s_tanf.c | ||
| s_tanhf.c | ||
| s_totalorderf.c | ||
| s_totalordermagf.c | ||
| s_truncf.c | ||
| s_ufromfpf.c | ||
| s_ufromfpxf.c | ||
| t_exp2f.h | ||
| w_expf_compat.c | ||