glibc/sysdeps/ieee754/ldbl-128ibm
Joseph Myers b72592e75f Fix log10 (1) in round-downward mode (bug 16977).
As with various other issues of this kind, bug 16977 is log10 (1)
wrongly returning -0 rather than +0 in round-downward mode because of
an implementation effectively in terms of log1p (x - 1).  This patch
fixes the issue in the same way used for log.

Tested x86_64 and x86 and ulps updated accordingly.  Also tested for
mips64 to confirm a fix was needed for ldbl-128 and to validate that
fix (also applied to ldbl-128ibm since that version of logl is
essentially the same as the ldbl-128 one).

	[BZ #16977]
	* sysdeps/i386/fpu/e_log10.S (__ieee754_log10): Take absolute
	value when x - 1 is zero.
	* sysdeps/i386/fpu/e_log10f.S (__ieee754_log10f): Likewise.
	* sysdeps/i386/fpu/e_log10l.S (__ieee754_log10l): Likewise.
	* sysdeps/ieee754/ldbl-128/e_log10l.c (__ieee754_log10l): Return
	0.0L for an argument of 1.0L.
	* sysdeps/ieee754/ldbl-128ibm/e_log10l.c (__ieee754_log10l):
	Likewise.
	* sysdeps/x86_64/fpu/e_log10l.S (__ieee754_log10l): Take absolute
	value when x - 1 is zero.
	* math/libm-test.inc (log10_test): Use ALL_RM_TEST.
	* sysdeps/i386/fpu/libm-test-ulps: Update.
	* sysdeps/x86_64/fpu/libm-test-ulps: Likewise.
2014-05-23 12:07:50 +00:00
..
e_acoshl.c
e_acosl.c
e_asinl.c
e_atan2l.c
e_atanhl.c
e_coshl.c
e_exp10l.c
e_expl.c
e_fmodl.c
e_gammal_r.c
e_hypotl.c
e_ilogbl.c
e_j0l.c
e_j1l.c
e_jnl.c
e_lgammal_r.c
e_log2l.c
e_log10l.c Fix log10 (1) in round-downward mode (bug 16977). 2014-05-23 12:07:50 +00:00
e_logl.c
e_powl.c
e_rem_pio2l.c
e_remainderl.c
e_sinhl.c
e_sqrtl.c
gamma_productl.c
ieee754.h
k_cosl.c
k_sincosl.c
k_sinl.c
k_tanl.c
ldbl2mpn.c
Makefile
math_ldbl.h
mpn2ldbl.c
printf_fphex.c
s_asinhl.c
s_atanl.c
s_cbrtl.c
s_ceill.c
s_copysignl.c
s_cosl.c
s_cprojl.c
s_ctanhl.c
s_ctanl.c
s_erfl.c
s_expm1l.c
s_fabsl.c
s_finitel.c
s_floorl.c
s_fmal.c
s_fpclassifyl.c
s_frexpl.c
s_isinf_nsl.c
s_isinfl.c
s_isnanl.c
s_issignalingl.c
s_llrintl.c
s_llroundl.c
s_log1pl.c
s_logbl.c
s_lrintl.c
s_lroundl.c
s_modfl.c
s_nearbyintl.c
s_nextafterl.c
s_nexttoward.c
s_nexttowardf.c
s_remquol.c
s_rintl.c
s_roundl.c
s_scalblnl.c
s_scalbnl.c
s_signbitl.c
s_sincosl.c
s_sinl.c
s_tanhl.c
s_tanl.c
s_truncl.c
strtold_l.c
t_sincosl.c
w_expl.c
x2y2m1l.c