d/freebsd-skq
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

ld80/s_expl.c:

Browse Source 
* In the special case x = -Inf or -NaN, use a micro-optimization
  to eliminate the need to access u.xbits.man.

* Fix an off-by-one for small arguments |x| < 0x1p-65.

ld128/s_expl.c:

* In the special case x = -Inf or -NaN, use a micro-optimization
  to eliminate the need to access u.xbits.manh and u.xbits.manl.

* Fix an off-by-one for small arguments |x| < 0x1p-114.

Obtained from:	bde
This commit is contained in:
Steve Kargl 2013-06-03 18:51:34 +00:00
parent 31407861b8
commit 1783063f18
2 changed files with 9 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
lib/msun/ld128/s_expl.c
View File

@ -234,18 +234,16 @@ expl(long double x)
ix = hx & 0x7fff;
if (ix >= BIAS + 13) { /* |x| >= 8192 or x is NaN */
if (ix == BIAS + LDBL_MAX_EXP) {
if (hx & 0x8000 && u.xbits.manh == 0 &&
u.xbits.manl == 0)
return (0.0L); /* x is -Inf */
if (hx & 0x8000) /* x is -Inf or -NaN */
return (-1 / x);
return (x + x); /* x is +Inf or NaN */
}
if (x > o_threshold)
return (huge * huge);
if (x < u_threshold)
return (tiny * tiny);
} else if (ix < BIAS - 115) { /* |x| < 0x1p-115 */
if (huge + x > 1.0L) /* trigger inexact iff x != 0 */
return (1.0L + x);
} else if (ix < BIAS - 114) { /* |x| < 0x1p-114 */
return (1 + x); /* 1 with inexact iff x != 0 */
}
/* Reduce x to (k*ln2 + endpoint[n2] + r1 + r2). */

12
lib/msun/ld80/s_expl.c
View File

@ -246,18 +246,16 @@ expl(long double x)
ix = hx & 0x7fff;
if (ix >= BIAS + 13) { /* |x| >= 8192 or x is NaN */
if (ix == BIAS + LDBL_MAX_EXP) {
if (hx & 0x8000 && u.xbits.man == 1ULL << 63)
return (0.0L); /* x is -Inf */
return (x + x); /* x is +Inf, NaN or unsupported */
if (hx & 0x8000) /* x is -Inf, -NaN or unsupported */
return (-1 / x);
return (x + x); /* x is +Inf, +NaN or unsupported */
}
if (x > o_threshold)
return (huge * huge);
if (x < u_threshold)
return (tiny * tiny);
} else if (ix < BIAS - 66) { /* |x| < 0x1p-66 */
/* includes pseudo-denormals */
if (huge + x > 1.0L) /* trigger inexact iff x != 0 */
return (1.0L + x);
} else if (ix < BIAS - 65) { /* |x| < 0x1p-65 (includes pseudos) */
return (1 + x); /* 1 with inexact iff x != 0 */
}
ENTERI();