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.

1 files changed, 15 insertions, 33 deletions

diff --git a/sysdeps/libm-ieee754/e_exp.c b/sysdeps/libm-ieee754/e_exp.c
index 660c5bc88d..5eae12a19c 100644
--- a/sysdeps/libm-ieee754/e_exp.c
+++ b/sysdeps/libm-ieee754/e_exp.c
@@ -76,8 +76,8 @@ __ieee754_exp (double x)
   /* Check for usual case.  */
   if (isless (x, himark) && isgreater (x, lomark))
     {
-      static const double TWO43 = 8796093022208.0;
-      static const double TWO52 = 4503599627370496.0;
+      static const double THREEp42 = 13194139533312.0;
+      static const double THREEp51 = 6755399441055744.0;
       /* 1/ln(2).  */
       static const double M_1_LN2 = 1.442695040888963387;
       /* ln(2), part 1 */
@@ -94,40 +94,22 @@ __ieee754_exp (double x)
       fesetround (FE_TONEAREST);
 
       /* Calculate n.  */
-      if (x >= 0)
-	{
-	  n = x * M_1_LN2 + TWO52;
-	  n -= TWO52;
-	}
-      else
-	{
-	  n = x * M_1_LN2 - TWO52;
-	  n += TWO52;
-	}
+      n = x * M_1_LN2 + THREEp51;
+      n -= THREEp51;
       x = x - n*M_LN2_0;
-      if (x >= 0)
-	{
-	  /* Calculate t/512.  */
-	  t = x + TWO43;
-	  t -= TWO43;
-	  x -= t;
-
-	  /* Compute tval = t.  */
-	  tval = (int) (t * 512.0);
-
-	  x -= __exp_deltatable[tval];
-	}
-      else
-	{
-	  /* As above, but x is negative.  */
-	  t = x - TWO43;
-	  t += TWO43;
-	  x -= t;
 
-	  tval = (int) (t * 512.0);
+      /* Calculate t/512.  */
+      t = x + THREEp42;
+      t -= THREEp42;
+      x -= t;
 
-	  x += __exp_deltatable[-tval];
-	}
+      /* Compute tval = t.  */
+      tval = (int) (t * 512.0);
+
+      if (t >= 0)
+	x -= __exp_deltatable[tval];
+      else
+	x += __exp_deltatable[-tval];
 
       /* Now, the variable x contains x + n*ln(2)_1.  */
       dely = n*M_LN2_1;