diff options
| author | Ulrich Drepper <drepper@redhat.com> | 1998-05-06 14:43:15 +0000 |
|---|---|---|
| committer | Ulrich Drepper <drepper@redhat.com> | 1998-05-06 14:43:15 +0000 |
| commit | 6600049466b586e3decaf24bd70c06b21382cf98 (patch) | |
| tree | 536fc5c8efe78035527223283bcd5a48ac259576 /sysdeps/libm-ieee754/s_exp2.c | |
| parent | b22fc5f5651706304ac09305ac3ee5bf84516378 (diff) | |
| download | glibc-6600049466b586e3decaf24bd70c06b21382cf98.tar glibc-6600049466b586e3decaf24bd70c06b21382cf98.tar.gz glibc-6600049466b586e3decaf24bd70c06b21382cf98.tar.bz2 glibc-6600049466b586e3decaf24bd70c06b21382cf98.zip | |
Update.
1998-04-16 07:42 Geoff Keating <geoffk@ozemail.com.au>
* Makeconfig [!build-static]: Link `static' binaries with libc_pic.a.
Still need *FLAGS-.o because we still sometimes build .o files.
* db2/Makefile: Don't build libndbm.a if !build-static.
1998-04-16 07:42 Geoff Keating <geoffk@ozemail.com.au>
* configure.in: New test for broken gcc on PowerPC.
* sysdeps/powerpc/atomicity.h: Use result of test.
* linuxthreads/sysdeps/powerpc/pt-machine.h: Use result of test.
* math/libm-test.c: Update many of the epsilon to match actual
performance.
* sysdeps/libm-ieee754/e_exp.c: Reduce the number of branches.
* sysdeps/libm-ieee754/e_expf.c: Likewise.
* sysdeps/libm-ieee754/s_exp2.c: Likewise.
* sysdeps/libm-ieee754/s_exp2f.c: Likewise.
* sysdeps/libm-ieee754/e_pow.c: Correct typo.
* sysdeps/powerpc/elf/libc-start.c: New file.
* sysdeps/powerpc/elf/start.S: New file, use libc-start.
* sysdeps/powerpc/elf/start.c: Delete.
* sysdeps/unix/sysv/linux/powerpc/Dist: Remove syscall.h
* sysdeps/unix/sysv/linux/powerpc/syscall.h: Delete. It was unused.
* sysdeps/unix/sysv/linux/powerpc/dl-sysdep.c: Correct previous
few patches.
Diffstat (limited to 'sysdeps/libm-ieee754/s_exp2.c')
| -rw-r--r-- | sysdeps/libm-ieee754/s_exp2.c | 14 |
1 files changed, 3 insertions, 11 deletions
diff --git a/sysdeps/libm-ieee754/s_exp2.c b/sysdeps/libm-ieee754/s_exp2.c index ead1ce89eb..875d4d6f2c 100644 --- a/sysdeps/libm-ieee754/s_exp2.c +++ b/sysdeps/libm-ieee754/s_exp2.c @@ -48,7 +48,7 @@ __ieee754_exp2 (double x) /* Check for usual case. */ if (isless (x, himark) && isgreater (x, lomark)) { - static const double TWO43 = 8796093022208.0; + static const double THREEp42 = 13194139533312.0; int tval, unsafe; double rx, x22, result; union ieee754_double ex2_u, scale_u; @@ -66,16 +66,8 @@ __ieee754_exp2 (double x) x = ex + t/512 + x1. First, calculate rx = ex + t/512. */ - if (x >= 0) - { - rx = x + TWO43; - rx -= TWO43; - } - else - { - rx = x - TWO43; - rx += TWO43; - } + rx = x + THREEp42; + rx -= THREEp42; x -= rx; /* Compute x=x1. */ /* Compute tval = (ex*512 + t)+256. Now, t = (tval mod 512)-256 and ex=tval/512 [that's mod, NOT %; and |