Fix undue underflow in ldbl-128ibm version of expl

git-mirror · Mar 26, 2012 · a3f6131 · a3f6131
1 parent d6a403f
commit a3f6131
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 2 deletions.
diff --git a/ChangeLog b/ChangeLog
@@ -1,3 +1,8 @@
+2012-03-26  Andreas Schwab  <schwab@linux-m68k.org>
+
+	* sysdeps/ieee754/ldbl-128ibm/e_expl.c (lomark): Adjust to take
+	denormal result into account.
+
 2012-03-25  Roland McGrath  <roland@hack.frob.com>
 
 	* posix/confstr.c (confstr): Lift RESTENVS definition to function scope.

diff --git a/sysdeps/ieee754/ldbl-128ibm/e_expl.c b/sysdeps/ieee754/ldbl-128ibm/e_expl.c
@@ -1,5 +1,5 @@
 /* Quad-precision floating point e^x.
-   Copyright (C) 1999,2004,2006, 2008, 2011 Free Software Foundation, Inc.
+   Copyright (C) 1999-2012 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Jakub Jelinek <jj@ultra.linux.cz>
    Partly based on double-precision code
@@ -74,7 +74,7 @@ static const long double C[] = {
 
 /* Largest integer x for which e^x underflows.  */
 #define lomark C[1]
--709.08956571282405153382846025171462914L,
+-744.44007192138121808966388925909996033L,
 
 /* 3x2^96 */
 #define THREEp96 C[2]