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_j1f.c | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) (limited to 'math/w_j1f.c') diff --git a/math/w_j1f.c b/math/w_j1f.c index bf7deb09f7..f70913d5ae 100644 --- a/math/w_j1f.c +++ b/math/w_j1f.c @@ -49,8 +49,11 @@ y1f (float x) return __kernel_standard_f (x, x, 111); } else if (x == 0.0f) - /* d = -one/(x-x) */ - return __kernel_standard_f (x, x, 110); + { + /* d = -one/(x-x) */ + feraiseexcept (FE_DIVBYZERO); + return __kernel_standard_f (x, x, 110); + } else if (_LIB_VERSION != _POSIX_) /* y1(x>X_TLOSS) */ return __kernel_standard_f (x, x, 137); -- cgit v1.2.3