about summary refs log tree commit diff
path: root/sysdeps
diff options
context:
space:
mode:
Diffstat (limited to 'sysdeps')
-rw-r--r--sysdeps/alpha/fpu/bits/fenv.h6
-rw-r--r--sysdeps/generic/bits/fenv.h4
-rw-r--r--sysdeps/i386/fpu/bits/fenv.h6
-rw-r--r--sysdeps/i386/fpu/bits/mathinline.h24
-rw-r--r--sysdeps/i386/fpu/fesetenv.c8
-rw-r--r--sysdeps/libm-ieee754/s_log2.c7
-rw-r--r--sysdeps/libm-ieee754/s_log2f.c7
-rw-r--r--sysdeps/m68k/fpu/bits/fenv.h6
-rw-r--r--sysdeps/powerpc/bits/fenv.h6
-rw-r--r--sysdeps/powerpc/test-arith.c4
-rw-r--r--sysdeps/sparc/sparc32/fpu/bits/fenv.h6
-rw-r--r--sysdeps/sparc/sparc64/fpu/bits/fenv.h6
12 files changed, 46 insertions, 44 deletions
diff --git a/sysdeps/alpha/fpu/bits/fenv.h b/sysdeps/alpha/fpu/bits/fenv.h
index 4482f0a2dd..2ccf149f0a 100644
--- a/sysdeps/alpha/fpu/bits/fenv.h
+++ b/sysdeps/alpha/fpu/bits/fenv.h
@@ -1,4 +1,4 @@
-/* Copyright (C) 1997 Free Software Foundation, Inc.
+/* Copyright (C) 1997, 1998 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -68,8 +68,8 @@ enum
 
 enum
   {
-    FE_TOWARDSZERO =	0,
-#define FE_TOWARDSZERO	FE_TOWARDSZERO
+    FE_TOWARDZERO =	0,
+#define FE_TOWARDZERO	FE_TOWARDZERO
 
     FE_DOWNWARD = 	1,
 #define FE_DOWNWARD	FE_DOWNWARD
diff --git a/sysdeps/generic/bits/fenv.h b/sysdeps/generic/bits/fenv.h
index c42540fa61..7552677d84 100644
--- a/sysdeps/generic/bits/fenv.h
+++ b/sysdeps/generic/bits/fenv.h
@@ -1,4 +1,4 @@
-/* Copyright (C) 1997 Free Software Foundation, Inc.
+/* Copyright (C) 1997, 1998 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -36,7 +36,7 @@
     FE_TONEAREST
     FE_DOWNWARD
     FE_UPWARD
-    FE_TOWARDSZERO
+    FE_TOWARDZERO
    We define no macro which signals no rounding mode is selectable.  */
 
 
diff --git a/sysdeps/i386/fpu/bits/fenv.h b/sysdeps/i386/fpu/bits/fenv.h
index 4d06b01be3..cc2aa41634 100644
--- a/sysdeps/i386/fpu/bits/fenv.h
+++ b/sysdeps/i386/fpu/bits/fenv.h
@@ -1,4 +1,4 @@
-/* Copyright (C) 1997 Free Software Foundation, Inc.
+/* Copyright (C) 1997, 1998 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -52,8 +52,8 @@ enum
 #define FE_DOWNWARD	FE_DOWNWARD
     FE_UPWARD = 0x800,
 #define FE_UPWARD	FE_UPWARD
-    FE_TOWARDSZERO = 0xc00
-#define FE_TOWARDSZERO	FE_TOWARDSZERO
+    FE_TOWARDZERO = 0xc00
+#define FE_TOWARDZERO	FE_TOWARDZERO
   };
 
 
diff --git a/sysdeps/i386/fpu/bits/mathinline.h b/sysdeps/i386/fpu/bits/mathinline.h
index 9917316f7a..77b6cd5e51 100644
--- a/sysdeps/i386/fpu/bits/mathinline.h
+++ b/sysdeps/i386/fpu/bits/mathinline.h
@@ -39,36 +39,36 @@
    better code.  */
 #  define isgreater(x, y) \
      ({ register char __result;						      \
-	__asm__ ("fucomip; seta %%al"					      \
+	__asm__ ("fucomip %%st(1), %%st; seta %%al"			      \
 		 : "=a" (__result) : "u" (y), "t" (x) : "cc", "st");	      \
 	__result; })
 #  define isgreaterequal(x, y) \
      ({ register char __result;						      \
-	__asm__ ("fucomip; setae %%al"					      \
+	__asm__ ("fucomip %%st(1), %%st; setae %%al"			      \
 		 : "=a" (__result) : "u" (y), "t" (x) : "cc", "st");	      \
 	__result; })
 
 #  define isless(x, y) \
      ({ register char __result;						      \
-	__asm__ ("fucomip; setb %%al"					      \
-		 : "=a" (__result) : "u" (y), "t" (x) : "cc", "st");	      \
+	__asm__ ("fucomip %%st(1), %%st; seta %%al"			      \
+		 : "=a" (__result) : "u" (x), "t" (y) : "cc", "st");	      \
 	__result; })
 
 #  define islessequal(x, y) \
      ({ register char __result;						      \
-	__asm__ ("fucomip; setbe %%al"					      \
-		 : "=a" (__result) : "u" (y), "t" (x) : "cc", "st");	      \
+	__asm__ ("fucomip %%st(1), %%st; setae %%al"			      \
+		 : "=a" (__result) : "u" (x), "t" (y) : "cc", "st");	      \
 	__result; })
 
 #  define islessgreater(x, y) \
      ({ register char __result;						      \
-	__asm__ ("fucomip; setne %%al"					      \
+	__asm__ ("fucomip %%st(1), %%st; setne %%al"			      \
 		 : "=a" (__result) : "u" (y), "t" (x) : "cc", "st");	      \
 	__result; })
 
 #  define isunordered(x, y) \
      ({ register char __result;						      \
-	__asm__ ("fucomip; setp %%al"					      \
+	__asm__ ("fucomip %%st(1), %%st; setp %%al"			      \
 		 : "=a" (__result) : "u" (y), "t" (x) : "cc", "st");	      \
 	__result; })
 # else
@@ -87,14 +87,14 @@
 
 #  define isless(x, y) \
      ({ register char __result;						      \
-	__asm__ ("fucompp; fnstsw; sahf; setb %%al"			      \
-		 : "=a" (__result) : "u" (y), "t" (x) : "cc", "st", "st(1)"); \
+	__asm__ ("fucompp; fnstsw; testb $0x45, %%ah; setz %%al"	      \
+		 : "=a" (__result) : "u" (x), "t" (y) : "cc", "st", "st(1)"); \
 	__result; })
 
 #  define islessequal(x, y) \
      ({ register char __result;						      \
-	__asm__ ("fucompp; fnstsw; sahf; setbe %%al"			      \
-		 : "=a" (__result) : "u" (y), "t" (x) : "cc", "st", "st(1)"); \
+	__asm__ ("fucompp; fnstsw; testb $0x05, %%ah; setz %%al"	      \
+		 : "=a" (__result) : "u" (x), "t" (y) : "cc", "st", "st(1)"); \
 	__result; })
 
 #  define islessgreater(x, y) \
diff --git a/sysdeps/i386/fpu/fesetenv.c b/sysdeps/i386/fpu/fesetenv.c
index e7300ea1e6..3faf01a568 100644
--- a/sysdeps/i386/fpu/fesetenv.c
+++ b/sysdeps/i386/fpu/fesetenv.c
@@ -1,5 +1,5 @@
 /* Install given floating-point environment.
-   Copyright (C) 1997 Free Software Foundation, Inc.
+   Copyright (C) 1997, 1998 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
@@ -50,7 +50,7 @@ fesetenv (const fenv_t *envp)
     }
   else if (envp == FE_NOMASK_ENV)
     {
-      temp.control_word &= ~(FE_ALL_EXCEPT | FE_TOWARDSZERO);
+      temp.control_word &= ~(FE_ALL_EXCEPT | FE_TOWARDZERO);
       temp.status_word &= ~FE_ALL_EXCEPT;
       temp.eip = 0;
       temp.cs_selector = 0;
@@ -60,9 +60,9 @@ fesetenv (const fenv_t *envp)
     }
   else
     {
-      temp.control_word &= ~(FE_ALL_EXCEPT | FE_TOWARDSZERO);
+      temp.control_word &= ~(FE_ALL_EXCEPT | FE_TOWARDZERO);
       temp.control_word |= (envp->control_word
-			    & (FE_ALL_EXCEPT | FE_TOWARDSZERO));
+			    & (FE_ALL_EXCEPT | FE_TOWARDZERO));
       temp.status_word &= ~FE_ALL_EXCEPT;
       temp.status_word |= envp->status_word & FE_ALL_EXCEPT;
       temp.eip = envp->eip;
diff --git a/sysdeps/libm-ieee754/s_log2.c b/sysdeps/libm-ieee754/s_log2.c
index 93c20b3db1..7379ce85e7 100644
--- a/sysdeps/libm-ieee754/s_log2.c
+++ b/sysdeps/libm-ieee754/s_log2.c
@@ -62,6 +62,7 @@ static const double
 #else
 static double
 #endif
+ln2 = 0.69314718055994530942,
 two54   =  1.80143985094819840000e+16,  /* 43500000 00000000 */
 Lg1 = 6.666666666666735130e-01,  /* 3FE55555 55555593 */
 Lg2 = 3.999999999940941908e-01,  /* 3FD99999 9997FA04 */
@@ -109,7 +110,7 @@ static double zero   =  0.0;
 	if((0x000fffff&(2+hx))<3) {	/* |f| < 2**-20 */
 	    if(f==zero) return dk;
 	    R = f*f*(0.5-0.33333333333333333*f);
-	    return dk-(R-f);
+	    return dk-(R-f)/ln2;
 	}
  	s = f/(2.0+f);
 	z = s*s;
@@ -122,9 +123,9 @@ static double zero   =  0.0;
 	R = t2+t1;
 	if(i>0) {
 	    hfsq=0.5*f*f;
-	    return dk-((hfsq-(s*(hfsq+R)))-f);
+	    return dk-((hfsq-(s*(hfsq+R)))-f)/ln2;
 	} else {
-	    return dk-((s*(f-R))-f);
+	    return dk-((s*(f-R))-f)/ln2;
 	}
 }
 
diff --git a/sysdeps/libm-ieee754/s_log2f.c b/sysdeps/libm-ieee754/s_log2f.c
index da9dfb8324..2377acd36a 100644
--- a/sysdeps/libm-ieee754/s_log2f.c
+++ b/sysdeps/libm-ieee754/s_log2f.c
@@ -23,6 +23,7 @@ static const float
 #else
 static float
 #endif
+ln2 = 0.69314718055994530942,
 two25 =    3.355443200e+07,	/* 0x4c000000 */
 Lg1 = 6.6666668653e-01,	/* 3F2AAAAB */
 Lg2 = 4.0000000596e-01,	/* 3ECCCCCD */
@@ -69,7 +70,7 @@ static float zero   =  0.0;
 	if((0x007fffff&(15+ix))<16) {	/* |f| < 2**-20 */
 	    if(f==zero) return dk;
 	    R = f*f*((float)0.5-(float)0.33333333333333333*f);
-	    return dk-(R-f);
+	    return dk-(R-f)/ln2;
 	}
  	s = f/((float)2.0+f);
 	z = s*s;
@@ -82,9 +83,9 @@ static float zero   =  0.0;
 	R = t2+t1;
 	if(i>0) {
 	    hfsq=(float)0.5*f*f;
-	    return dk-((hfsq-(s*(hfsq+R)))-f);
+	    return dk-((hfsq-(s*(hfsq+R)))-f)/ln2;
 	} else {
-	    return dk-((s*(f-R))-f);
+	    return dk-((s*(f-R))-f)/ln2;
 	}
 }
 weak_alias (__log2f, log2f)
diff --git a/sysdeps/m68k/fpu/bits/fenv.h b/sysdeps/m68k/fpu/bits/fenv.h
index b1608b9ddc..e1a278d333 100644
--- a/sysdeps/m68k/fpu/bits/fenv.h
+++ b/sysdeps/m68k/fpu/bits/fenv.h
@@ -1,4 +1,4 @@
-/* Copyright (C) 1997 Free Software Foundation, Inc.
+/* Copyright (C) 1997, 1998 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -47,8 +47,8 @@ enum
   {
     FE_TONEAREST = 0,
 #define FE_TONEAREST	FE_TONEAREST
-    FE_TOWARDSZERO = 1 << 4,
-#define FE_TOWARDSZERO	FE_TOWARDSZERO
+    FE_TOWARDZERO = 1 << 4,
+#define FE_TOWARDZERO	FE_TOWARDZERO
     FE_DOWNWARD = 2 << 4,
 #define FE_DOWNWARD	FE_DOWNWARD
     FE_UPWARD = 3 << 4
diff --git a/sysdeps/powerpc/bits/fenv.h b/sysdeps/powerpc/bits/fenv.h
index 14fbda8db0..4f394fd64f 100644
--- a/sysdeps/powerpc/bits/fenv.h
+++ b/sysdeps/powerpc/bits/fenv.h
@@ -1,4 +1,4 @@
-/* Copyright (C) 1997 Free Software Foundation, Inc.
+/* Copyright (C) 1997, 1998 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -106,8 +106,8 @@ enum
   {
     FE_TONEAREST = 0,
 #define FE_TONEAREST	FE_TONEAREST
-    FE_TOWARDSZERO = 1,
-#define FE_TOWARDSZERO	FE_TOWARDSZERO
+    FE_TOWARDZERO = 1,
+#define FE_TOWARDZERO	FE_TOWARDZERO
     FE_UPWARD = 2,
 #define FE_UPWARD	FE_UPWARD
     FE_DOWNWARD = 3,
diff --git a/sysdeps/powerpc/test-arith.c b/sysdeps/powerpc/test-arith.c
index 9e1be88098..3cefc3b8cc 100644
--- a/sysdeps/powerpc/test-arith.c
+++ b/sysdeps/powerpc/test-arith.c
@@ -1,5 +1,5 @@
 /* Test floating-point arithmetic operations.
-   Copyright (C) 1997 Free Software Foundation, Inc.
+   Copyright (C) 1997, 1998 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -574,7 +574,7 @@ int main(int argc, char **argv)
   fesetenv(FE_DFL_ENV);
   fesetround(FE_TONEAREST);
   fegetenv(rmodes+0);
-  fesetround(FE_TOWARDSZERO);
+  fesetround(FE_TOWARDZERO);
   fegetenv(rmodes+1);
   fesetround(FE_UPWARD);
   fegetenv(rmodes+2);
diff --git a/sysdeps/sparc/sparc32/fpu/bits/fenv.h b/sysdeps/sparc/sparc32/fpu/bits/fenv.h
index 9fb33cef28..6858caac97 100644
--- a/sysdeps/sparc/sparc32/fpu/bits/fenv.h
+++ b/sysdeps/sparc/sparc32/fpu/bits/fenv.h
@@ -1,4 +1,4 @@
-/* Copyright (C) 1997 Free Software Foundation, Inc.
+/* Copyright (C) 1997, 1998 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -47,8 +47,8 @@ enum
   {
     FE_TONEAREST = 	(0U << 30),
 #define FE_TONEAREST	FE_TONEAREST
-    FE_TOWARDSZERO = 	(1U << 30),
-#define FE_TOWARDSZERO	FE_TOWARDSZERO
+    FE_TOWARDZERO = 	(1U << 30),
+#define FE_TOWARDZERO	FE_TOWARDZERO
     FE_UPWARD = 	(2U << 30),
 #define FE_UPWARD	FE_UPWARD
     FE_DOWNWARD = 	(3U << 30)
diff --git a/sysdeps/sparc/sparc64/fpu/bits/fenv.h b/sysdeps/sparc/sparc64/fpu/bits/fenv.h
index 5728e2b588..b864abdedf 100644
--- a/sysdeps/sparc/sparc64/fpu/bits/fenv.h
+++ b/sysdeps/sparc/sparc64/fpu/bits/fenv.h
@@ -1,4 +1,4 @@
-/* Copyright (C) 1997 Free Software Foundation, Inc.
+/* Copyright (C) 1997, 1998 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -47,8 +47,8 @@ enum
   {
     FE_TONEAREST = 	(0U << 30),
 #define FE_TONEAREST	FE_TONEAREST
-    FE_TOWARDSZERO = 	(1U << 30),
-#define FE_TOWARDSZERO	FE_TOWARDSZERO
+    FE_TOWARDZERO = 	(1U << 30),
+#define FE_TOWARDZERO	FE_TOWARDZERO
     FE_UPWARD = 	(2U << 30),
 #define FE_UPWARD	FE_UPWARD
     FE_DOWNWARD = 	(3U << 30)