From 699ff83712b5796ac50ed332d9dad55d38450e81 Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Wed, 4 Dec 2013 14:39:37 +0000 Subject: Fix Bessel function error handling (bug 6807, bug 15901). --- ChangeLog | 28 ++++++++++++++++++++++++++++ 1 file changed, 28 insertions(+) (limited to 'ChangeLog') diff --git a/ChangeLog b/ChangeLog index dc2fd4c9e9..cf60c7741d 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,3 +1,31 @@ +2013-12-04 Joseph Myers + + [BZ #6807] + [BZ #15901] + * math/w_j0.c (y0): Raise FE_DIVBYZERO on zero argument. + * math/w_j0f.c (y0f): Likewise. + * math/w_j0l.c (__y0l): Likewise. + * math/w_j1.c (y1): Likewise. + * math/w_j1f.c (y1f): Likewise. + * math/w_j1l.c (__y1l): Likewise + * math/w_jn.c (yn): Likewise. + * math/w_jnf.c (ynf): Likewise. + * sysdeps/ieee754/k_standard.c (__kernel_standard): Use ERANGE for + Bessel function pole errors in _POSIX_ mode. Use NAN as return + value for Bessel function domain errors outside _SVID_ mode. + Adjust sign of return value for yn (negative integer, 0). + * sysdeps/ieee754/ldbl-128/e_jnl.c (__ieee754_ynl): Use division + by zero in return for negative x and set sign appropriately for + negative n. + * sysdeps/ieee754/ldbl-128ibm/e_jnl.c (__ieee754_ynl): Likewise. + * sysdeps/ieee754/ldbl-96/e_jnl.c (__ieee754_ynl): Likewise. + * math/libm-test.inc (y0_test_data): Add more tests and adjust + expectations in error cases. + (y1_test_data): Likewise. + (yn_test_data): Likewise. + * sysdeps/i386/fpu/libm-test-ulps: Update. + * sysdeps/x86_64/fpu/libm-test-ulps: Likewise. + 2013-11-12 Ulrich Weigand * sysdeps/unix/sysv/linux/powerpc/Makefile (abi-variants): Rename -- cgit 1.4.1