Module Name: src
Committed By: christos
Date: Wed Nov 6 16:49:21 UTC 2013
Modified Files:
src/lib/libm/src: s_exp2.c
Log Message:
simplify.
To generate a diff of this commit:
cvs rdiff -u -r1.2 -r1.3 src/lib/libm/src/s_exp2.c
Please note that diffs are not public domain; they are subject to the
copyright notices on the relevant files.
Modified files:
Index: src/lib/libm/src/s_exp2.c
diff -u src/lib/libm/src/s_exp2.c:1.2 src/lib/libm/src/s_exp2.c:1.3
--- src/lib/libm/src/s_exp2.c:1.2 Mon Jan 11 18:38:24 2010
+++ src/lib/libm/src/s_exp2.c Wed Nov 6 11:49:21 2013
@@ -25,7 +25,7 @@
*/
#include <sys/cdefs.h>
-__RCSID("$NetBSD: s_exp2.c,v 1.2 2010/01/11 23:38:24 christos Exp $");
+__RCSID("$NetBSD: s_exp2.c,v 1.3 2013/11/06 16:49:21 christos Exp $");
#ifdef __FBSDID
__FBSDID("$FreeBSD: src/lib/msun/src/s_exp2.c,v 1.7 2008/02/22 02:27:34 das Exp $");
#endif
@@ -343,9 +343,9 @@ static const double tbl[TBLSIZE * 2] = {
double
exp2(double x)
{
- double r, t, twopk, twopkp1000, z;
+ double r, t, twopk, z;
uint32_t hx, ix, lx, i0;
- int k;
+ int k, big;
/* Filter out exceptional cases. */
GET_HIGH_WORD(hx,x);
@@ -378,19 +378,20 @@ exp2(double x)
/* Compute r = exp2(y) = exp2t[i0] * p(z - eps[i]). */
t = tbl[i0]; /* exp2t[i0] */
z -= tbl[i0 + 1]; /* eps[i0] */
- if (k >= -1021 << 20)
+ big = k >= -1021 << 20;
+ if (big)
INSERT_WORDS(twopk, 0x3ff00000 + k, 0);
else
- INSERT_WORDS(twopkp1000, 0x3ff00000 + k + (1000 << 20), 0);
+ INSERT_WORDS(twopk, 0x3ff00000 + k + (1000 << 20), 0);
r = t + t * z * (P1 + z * (P2 + z * (P3 + z * (P4 + z * P5))));
/* Scale by 2**(k>>20). */
- if(k >= -1021 << 20) {
+ if (big) {
if (k == 1024 << 20)
return (r * 2.0 * 0x1p1023);
return (r * twopk);
} else {
- return (r * twopkp1000 * twom1000);
+ return (r * twopk * twom1000);
}
}
