From 59e53a78985be92134f126c5597d15a9998001d1 Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Thu, 9 Jun 2016 17:23:51 +0000 Subject: Fix ldbl-128 expm1l (sNaN) (bug 20232). The ldbl-128 version of expm1l returns sNaN for sNaN input. This patch fixes it to add such inputs to themselves so that qNaN is returned in this case. Tested for mips64. [BZ #20232] * sysdeps/ieee754/ldbl-128/s_expm1l.c (__expm1l): Add NaN input to itself. --- ChangeLog | 6 ++++++ sysdeps/ieee754/ldbl-128/s_expm1l.c | 4 ++-- 2 files changed, 8 insertions(+), 2 deletions(-) diff --git a/ChangeLog b/ChangeLog index f72f22e2df..2d88f74d6b 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,3 +1,9 @@ +2016-06-09 Joseph Myers + + [BZ #20232] + * sysdeps/ieee754/ldbl-128/s_expm1l.c (__expm1l): Add NaN input to + itself. + 2016-06-09 H.J. Lu * sysdeps/x86_64/start.S (_start): Always indirect branch to diff --git a/sysdeps/ieee754/ldbl-128/s_expm1l.c b/sysdeps/ieee754/ldbl-128/s_expm1l.c index b1100a9d0e..da902efafd 100644 --- a/sysdeps/ieee754/ldbl-128/s_expm1l.c +++ b/sysdeps/ieee754/ldbl-128/s_expm1l.c @@ -111,8 +111,8 @@ __expm1l (long double x) /* Infinity (which must be negative infinity). */ if (((ix & 0xffff) | u.parts32.w1 | u.parts32.w2 | u.parts32.w3) == 0) return -1.0L; - /* NaN. No invalid exception. */ - return x; + /* NaN. Invalid exception if signaling. */ + return x + x; } /* expm1(+- 0) = +- 0. */ -- cgit v1.2.3