summaryrefslogtreecommitdiff
path: root/libc/sysdeps/ieee754/ldbl-96/s_nextafterl.c
diff options
context:
space:
mode:
authorjoseph <joseph@7b3dc134-2b1b-0410-93df-9e9f96275f8d>2007-04-17 11:15:02 +0000
committerjoseph <joseph@7b3dc134-2b1b-0410-93df-9e9f96275f8d>2007-04-17 11:15:02 +0000
commit1aab636c1bf39a768891da3abcb72fac07e25329 (patch)
tree12dee84a1a058b1790287d152efa91bc4f9c6049 /libc/sysdeps/ieee754/ldbl-96/s_nextafterl.c
parent26fd795bec9789d6226b8643d8ad4a6998ce56f8 (diff)
downloadeglibc2-1aab636c1bf39a768891da3abcb72fac07e25329.tar.gz
Merge changes between r1968 and r2007 from /fsf/trunk.
git-svn-id: svn://svn.eglibc.org/trunk@2008 7b3dc134-2b1b-0410-93df-9e9f96275f8d
Diffstat (limited to 'libc/sysdeps/ieee754/ldbl-96/s_nextafterl.c')
-rw-r--r--libc/sysdeps/ieee754/ldbl-96/s_nextafterl.c20
1 files changed, 10 insertions, 10 deletions
diff --git a/libc/sysdeps/ieee754/ldbl-96/s_nextafterl.c b/libc/sysdeps/ieee754/ldbl-96/s_nextafterl.c
index 1798261b2..ed0877d6d 100644
--- a/libc/sysdeps/ieee754/ldbl-96/s_nextafterl.c
+++ b/libc/sysdeps/ieee754/ldbl-96/s_nextafterl.c
@@ -26,7 +26,7 @@ static char rcsid[] = "$NetBSD: $";
*/
#include "math.h"
-#include "math_private.h"
+#include <math_private.h>
#ifdef __STDC__
long double __nextafterl(long double x, long double y)
@@ -48,9 +48,12 @@ static char rcsid[] = "$NetBSD: $";
return x+y;
if(x==y) return y; /* x=y, return y */
if((ix|hx|lx)==0) { /* x == 0 */
+ long double u;
SET_LDOUBLE_WORDS(x,esy&0x8000,0,1);/* return +-minsubnormal */
- y = x*x;
- if(y==x) return y; else return x; /* raise underflow flag */
+ u = math_opt_barrier (x);
+ u = u * u;
+ math_force_eval (u); /* raise underflow flag */
+ return x;
}
if(esx<0x8000) { /* x > 0 */
if(ix>iy||((ix==iy) && (hx>hy||((hx==hy)&&(lx>ly))))) {
@@ -85,13 +88,10 @@ static char rcsid[] = "$NetBSD: $";
}
}
esy = esx&0x7fff;
- if(esy==0x7fff) return x+x; /* overflow */
- if(esy==0) { /* underflow */
- y = x*x;
- if(y!=x) { /* raise underflow flag */
- SET_LDOUBLE_WORDS(y,esx,hx,lx);
- return y;
- }
+ if(esy==0x7fff) return x+x; /* overflow */
+ if(esy==0) {
+ long double u = x*x; /* underflow */
+ math_force_eval (u); /* raise underflow flag */
}
SET_LDOUBLE_WORDS(x,esx,hx,lx);
return x;