qiurui
/
glibc
mirror of git://sourceware.org/git/glibc.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

* sysdeps/ieee754/ldbl-128ibm/s_nextafterl.c (__nextafterl): Use 

math_opt_barrier and math_force_eval macros.
This commit is contained in:
Ulrich Drepper 2007-04-16 20:42:23 +00:00
parent 3e336a8754
commit 7a82b74f17
5 changed files with 43 additions and 38 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
ChangeLog
View File

@ -1,3 +1,8 @@
2007-04-01 Jakub Jelinek <jakub@redhat.com>
* sysdeps/ieee754/ldbl-128ibm/s_nextafterl.c (__nextafterl): Use
math_opt_barrier and math_force_eval macros.
2007-03-27 Jakub Jelinek <jakub@redhat.com> 2007-03-27 Jakub Jelinek <jakub@redhat.com>
[BZ #3306] [BZ #3306]

18
sysdeps/i386/fpu/s_nextafterl.c
View File

@ -27,7 +27,7 @@ static char rcsid[] = "$NetBSD: $";
*/ */
#include "math.h" #include "math.h"
#include "math_private.h" #include <math_private.h>
#ifdef __STDC__ #ifdef __STDC__
long double __nextafterl(long double x, long double y) long double __nextafterl(long double x, long double y)
@ -52,9 +52,12 @@ static char rcsid[] = "$NetBSD: $";
return x+y; return x+y;
if(x==y) return y; /* x=y, return y */ if(x==y) return y; /* x=y, return y */
if((ix|hx|lx)==0) { /* x == 0 */ if((ix|hx|lx)==0) { /* x == 0 */
long double u;
SET_LDOUBLE_WORDS(x,esy&0x8000,0,1);/* return +-minsubnormal */ SET_LDOUBLE_WORDS(x,esy&0x8000,0,1);/* return +-minsubnormal */
y = x*x; u = math_opt_barrier (x);
if(y==x) return y; else return x; /* raise underflow flag */ u = u * u;
math_force_eval (u); /* raise underflow flag */
return x;
} }
if(esx>=0) { /* x > 0 */ if(esx>=0) { /* x > 0 */
if(esx>esy||((esx==esy) && (hx>hy||((hx==hy)&&(lx>ly))))) { if(esx>esy||((esx==esy) && (hx>hy||((hx==hy)&&(lx>ly))))) {
@ -109,12 +112,9 @@ static char rcsid[] = "$NetBSD: $";
} }
esy = esx&0x7fff; esy = esx&0x7fff;
if(esy==0x7fff) return x+x; /* overflow */ if(esy==0x7fff) return x+x; /* overflow */
if(esy==0) { /* underflow */ if(esy==0) {
y = x*x; long double u = x*x; /* underflow */
if(y!=x) { /* raise underflow flag */ math_force_eval (u); /* raise underflow flag */
SET_LDOUBLE_WORDS(y,esx,hx,lx);
return y;
}
} }
SET_LDOUBLE_WORDS(x,esx,hx,lx); SET_LDOUBLE_WORDS(x,esx,hx,lx);
return x; return x;

18
sysdeps/ieee754/ldbl-128/s_nextafterl.c
View File

@ -25,7 +25,7 @@ static char rcsid[] = "$NetBSD: $";
*/ */
#include "math.h" #include "math.h"
#include "math_private.h" #include <math_private.h>
#ifdef __STDC__ #ifdef __STDC__
long double __nextafterl(long double x, long double y) long double __nextafterl(long double x, long double y)
@ -47,9 +47,12 @@ static char rcsid[] = "$NetBSD: $";
return x+y; return x+y;
if(x==y) return y; /* x=y, return y */ if(x==y) return y; /* x=y, return y */
if((ix|lx)==0) { /* x == 0 */ if((ix|lx)==0) { /* x == 0 */
long double u;
SET_LDOUBLE_WORDS64(x,hy&0x8000000000000000ULL,1);/* return +-minsubnormal */ SET_LDOUBLE_WORDS64(x,hy&0x8000000000000000ULL,1);/* return +-minsubnormal */
y = x*x; u = math_opt_barrier (x);
if(y==x) return y; else return x; /* raise underflow flag */ u = u * u;
math_force_eval (u); /* raise underflow flag */
return x;
} }
if(hx>=0) { /* x > 0 */ if(hx>=0) { /* x > 0 */
if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */ if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */
@ -70,12 +73,9 @@ static char rcsid[] = "$NetBSD: $";
} }
hy = hx&0x7fff000000000000LL; hy = hx&0x7fff000000000000LL;
if(hy==0x7fff000000000000LL) return x+x;/* overflow */ if(hy==0x7fff000000000000LL) return x+x;/* overflow */
if(hy==0) { /* underflow */ if(hy==0) {
y = x*x; long double u = x*x; /* underflow */
if(y!=x) { /* raise underflow flag */ math_force_eval (u); /* raise underflow flag */
SET_LDOUBLE_WORDS64(y,hx,lx);
return y;
}
} }
SET_LDOUBLE_WORDS64(x,hx,lx); SET_LDOUBLE_WORDS64(x,hx,lx);
return x; return x;

20
sysdeps/ieee754/ldbl-128ibm/s_nextafterl.c
View File

@ -24,8 +24,8 @@ static char rcsid[] = "$NetBSD: $";
* Special cases: * Special cases:
*/ */
#include "math.h" #include <math.h>
#include "math_private.h" #include <math_private.h>
#include <math_ldbl_opt.h> #include <math_ldbl_opt.h>
#ifdef __STDC__ #ifdef __STDC__
@ -53,9 +53,12 @@ static char rcsid[] = "$NetBSD: $";
if(x==y) if(x==y)
return y; /* x=y, return y */ return y; /* x=y, return y */
if(ihx == 0 && ilx == 0) { /* x == 0 */ if(ihx == 0 && ilx == 0) { /* x == 0 */
long double u;
SET_LDOUBLE_WORDS64(x,hy&0x8000000000000000ULL,1);/* return +-minsubnormal */ SET_LDOUBLE_WORDS64(x,hy&0x8000000000000000ULL,1);/* return +-minsubnormal */
y = x*x; u = math_opt_barrier (u);
if(y==x) return y; else return x; /* raise underflow flag */ u = u * u;
math_force_eval (u); /* raise underflow flag */
return x;
} }
if(ihx>=0) { /* x > 0 */ if(ihx>=0) { /* x > 0 */
if(ihx>ihy||((ihx==ihy)&&(ilx>ily))) { /* x > y, x -= ulp */ if(ihx>ihy||((ihx==ihy)&&(ilx>ily))) { /* x > y, x -= ulp */
@ -93,12 +96,9 @@ static char rcsid[] = "$NetBSD: $";
} }
hy = hx&0x7ff0000000000000LL; hy = hx&0x7ff0000000000000LL;
if(hy==0x7ff0000000000000LL) return x+x;/* overflow */ if(hy==0x7ff0000000000000LL) return x+x;/* overflow */
if(hy==0) { /* underflow */ if(hy==0) {
y = x*x; long double u = x * x; /* underflow */
if(y!=x) { /* raise underflow flag */ math_force_eval (u); /* raise underflow flag */
SET_LDOUBLE_WORDS64(y,hx,lx);
return y;
}
} }
SET_LDOUBLE_WORDS64(x,hx,lx); SET_LDOUBLE_WORDS64(x,hx,lx);
return x; return x;

18
sysdeps/ieee754/ldbl-96/s_nextafterl.c
View File

@ -26,7 +26,7 @@ static char rcsid[] = "$NetBSD: $";
*/ */
#include "math.h" #include "math.h"
#include "math_private.h" #include <math_private.h>
#ifdef __STDC__ #ifdef __STDC__
long double __nextafterl(long double x, long double y) long double __nextafterl(long double x, long double y)
@ -48,9 +48,12 @@ static char rcsid[] = "$NetBSD: $";
return x+y; return x+y;
if(x==y) return y; /* x=y, return y */ if(x==y) return y; /* x=y, return y */
if((ix|hx|lx)==0) { /* x == 0 */ if((ix|hx|lx)==0) { /* x == 0 */
long double u;
SET_LDOUBLE_WORDS(x,esy&0x8000,0,1);/* return +-minsubnormal */ SET_LDOUBLE_WORDS(x,esy&0x8000,0,1);/* return +-minsubnormal */
y = x*x; u = math_opt_barrier (x);
if(y==x) return y; else return x; /* raise underflow flag */ u = u * u;
math_force_eval (u); /* raise underflow flag */
return x;
} }
if(esx<0x8000) { /* x > 0 */ if(esx<0x8000) { /* x > 0 */
if(ix>iy||((ix==iy) && (hx>hy||((hx==hy)&&(lx>ly))))) { if(ix>iy||((ix==iy) && (hx>hy||((hx==hy)&&(lx>ly))))) {
@ -86,12 +89,9 @@ static char rcsid[] = "$NetBSD: $";
} }
esy = esx&0x7fff; esy = esx&0x7fff;
if(esy==0x7fff) return x+x; /* overflow */ if(esy==0x7fff) return x+x; /* overflow */
if(esy==0) { /* underflow */ if(esy==0) {
y = x*x; long double u = x*x; /* underflow */
if(y!=x) { /* raise underflow flag */ math_force_eval (u); /* raise underflow flag */
SET_LDOUBLE_WORDS(y,esx,hx,lx);
return y;
}
} }
SET_LDOUBLE_WORDS(x,esx,hx,lx); SET_LDOUBLE_WORDS(x,esx,hx,lx);
return x; return x;