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_j1.c | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) (limited to 'math/w_j1.c') diff --git a/math/w_j1.c b/math/w_j1.c index 743c940..a9fb7ae 100644 --- a/math/w_j1.c +++ b/math/w_j1.c @@ -51,8 +51,11 @@ y1 (double x) return __kernel_standard (x, x, 11); } else if (x == 0.0) - /* d = -one/(x-x) */ - return __kernel_standard (x, x, 10); + { + /* d = -one/(x-x) */ + feraiseexcept (FE_DIVBYZERO); + return __kernel_standard (x, x, 10); + } else if (_LIB_VERSION != _POSIX_) /* y1(x>X_TLOSS) */ return __kernel_standard (x, x, 37); -- cgit v1.1