Source-Changes-HG archive
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index][Old Index]
[src/trunk]: src/lib/libm/src simplify.
details: https://anonhg.NetBSD.org/src/rev/dcd7cbbad5d7
branches: trunk
changeset: 791136:dcd7cbbad5d7
user: christos <christos%NetBSD.org@localhost>
date: Wed Nov 06 16:49:21 2013 +0000
description:
simplify.
diffstat:
lib/libm/src/s_exp2.c | 15 ++++++++-------
1 files changed, 8 insertions(+), 7 deletions(-)
diffs (49 lines):
diff -r 5911069f7c3b -r dcd7cbbad5d7 lib/libm/src/s_exp2.c
--- a/lib/libm/src/s_exp2.c Wed Nov 06 16:31:08 2013 +0000
+++ b/lib/libm/src/s_exp2.c Wed Nov 06 16:49:21 2013 +0000
@@ -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 @@
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 @@
/* 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);
}
}
Home |
Main Index |
Thread Index |
Old Index