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). --- math/w_j1l.c | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) (limited to 'math/w_j1l.c') diff --git a/math/w_j1l.c b/math/w_j1l.c index c4f8111..e4242ec 100644 --- a/math/w_j1l.c +++ b/math/w_j1l.c @@ -49,8 +49,11 @@ __y1l (long double x) return __kernel_standard_l (x, x, 211); } else if (x == 0.0L) - /* d = -one/(x-x) */ - return __kernel_standard_l (x, x, 210); + { + /* d = -one/(x-x) */ + feraiseexcept (FE_DIVBYZERO); + return __kernel_standard_l (x, x, 210); + } else if (_LIB_VERSION != _POSIX_) /* y1(x>X_TLOSS) */ return __kernel_standard_l (x, x, 237); -- cgit v1.1