system/glibc
3
0
Fork 0
Code Issues Projects Releases Wiki Activity

Clean up ldbl-128 / ldbl-128ibm expm1l dead code (bug 16415).

Browse Source 
The ldbl-128 and ldbl-128ibm expm1l implementations have code to
handle +Inf and finite arguments above an overflow threshold.  Since
they now use __expl for large positive arguments to fix other
problems, this code is unreachable; this patch removes it.

Tested for mips64 and powerpc.

	[BZ #16415]
	* sysdeps/ieee754/ldbl-128/s_expm1l.c (maxlog): Remove variable.
	(__expm1l): Remove code to handle positive infinity and overflow.
	* sysdeps/ieee754/ldbl-128ibm/s_expm1l.c (maxlog): Remove
	variable.
	(__expm1l): Remove code to handle positive infinity and overflow.
This commit is contained in:
Joseph Myers 2015-09-16 16:42:46 +00:00
parent 560b04462f
commit 8124ac3e73
4 changed files with 19 additions and 38 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
ChangeLog
View File

@ -1,3 +1,12 @@
2015-09-16 Joseph Myers <joseph@codesourcery.com>
[BZ #16415]
* sysdeps/ieee754/ldbl-128/s_expm1l.c (maxlog): Remove variable.
(__expm1l): Remove code to handle positive infinity and overflow.
* sysdeps/ieee754/ldbl-128ibm/s_expm1l.c (maxlog): Remove
variable.
(__expm1l): Remove code to handle positive infinity and overflow.
2015-09-16 Paul Eggert <eggert@cs.ucla.edu>
Paul Pluzhnikov <ppluzhnikov@google.com>

12
NEWS
View File

@ -9,12 +9,12 @@ Version 2.23
* The following bugs are resolved with this release:
2542, 2543, 2558, 2898, 14341, 14912, 15786, 15918, 16141, 16296, 16517,
16519, 16520, 16521, 16734, 16973, 16985, 17243, 17244, 17787, 17905,
18084, 18086, 18240, 18265, 18370, 18421, 18480, 18525, 18595, 18610,
18618, 18647, 18661, 18674, 18675, 18681, 18757, 18778, 18781, 18787,
18789, 18790, 18795, 18796, 18820, 18823, 18824, 18857, 18863, 18870,
18872, 18873, 18875, 18887, 18921, 18952, 18961, 18966, 18967.
2542, 2543, 2558, 2898, 14341, 14912, 15786, 15918, 16141, 16296, 16415,
16517, 16519, 16520, 16521, 16734, 16973, 16985, 17243, 17244, 17787,
17905, 18084, 18086, 18240, 18265, 18370, 18421, 18480, 18525, 18595,
18610, 18618, 18647, 18661, 18674, 18675, 18681, 18757, 18778, 18781,
18787, 18789, 18790, 18795, 18796, 18820, 18823, 18824, 18857, 18863,
18870, 18872, 18873, 18875, 18887, 18921, 18952, 18961, 18966, 18967.
* The obsolete header <regexp.h> has been removed. Programs that require
this header must be updated to use <regex.h> instead.

18
sysdeps/ieee754/ldbl-128/s_expm1l.c
View File

@ -84,8 +84,6 @@ static const long double
C1 = 6.93145751953125E-1L,
C2 = 1.428606820309417232121458176568075500134E-6L,
/* ln (2^16384 * (1 - 2^-113)) */
maxlog = 1.1356523406294143949491931077970764891253E4L,
/* ln 2^-114 */
minarg = -7.9018778583833765273564461846232128760607E1L, big = 1e4932L;
@ -110,14 +108,9 @@ __expm1l (long double x)
}
if (ix >= 0x7fff0000)
{
/* Infinity. */
/* Infinity (which must be negative infinity). */
if (((ix & 0xffff) | u.parts32.w1 | u.parts32.w2 | u.parts32.w3) == 0)
{
if (sign)
return -1.0L;
else
return x;
}
return -1.0L;
/* NaN. No invalid exception. */
return x;
}
@ -126,13 +119,6 @@ __expm1l (long double x)
if ((ix == 0) && (u.parts32.w1 | u.parts32.w2 | u.parts32.w3) == 0)
return x;
/* Overflow. */
if (x > maxlog)
{
__set_errno (ERANGE);
return (big * big);
}
/* Minimum value. */
if (x < minarg)
return (4.0/big - 1.0L);

18
sysdeps/ieee754/ldbl-128ibm/s_expm1l.c
View File

@ -82,8 +82,6 @@ static const long double
C1 = 6.93145751953125E-1L,
C2 = 1.428606820309417232121458176568075500134E-6L,
/* ln (2^16384 * (1 - 2^-113)) */
maxlog = 1.1356523406294143949491931077970764891253E4L,
/* ln 2^-114 */
minarg = -7.9018778583833765273564461846232128760607E1L, big = 1e290L;
@ -105,14 +103,9 @@ __expm1l (long double x)
return __expl (x);
if (ix >= 0x7ff00000)
{
/* Infinity. */
/* Infinity (which must be negative infinity). */
if (((ix - 0x7ff00000) | lx) == 0)
{
if (sign)
return -1.0L;
else
return x;
}
return -1.0L;
/* NaN. No invalid exception. */
return x;
}
@ -121,13 +114,6 @@ __expm1l (long double x)
if ((ix | lx) == 0)
return x;
/* Overflow. */
if (x > maxlog)
{
__set_errno (ERANGE);
return (big * big);
}
/* Minimum value. */
if (x < minarg)
return (4.0/big - 1.0L);