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