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_jnf.c | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) (limited to 'math/w_jnf.c') diff --git a/math/w_jnf.c b/math/w_jnf.c index c9114e7b09..cb1aab8e73 100644 --- a/math/w_jnf.c +++ b/math/w_jnf.c @@ -49,8 +49,11 @@ ynf (int n, float x) return __kernel_standard_f (n, x, 113); } else if (x == 0.0) - /* d = -one/(x-x) */ - return __kernel_standard_f (n, x, 112); + { + /* d = -one/(x-x) */ + feraiseexcept (FE_DIVBYZERO); + return __kernel_standard_f (n, x, 112); + } else if (_LIB_VERSION != _POSIX_) /* yn(n,x>X_TLOSS) */ return __kernel_standard_f (n, x, 139); -- cgit 1.4.1