1 /* s_nextafterl.c -- long double version of s_nextafter.c.
2 * Special version for i387.
3 * Conversion to long double by Ulrich Drepper,
4 * Cygnus Support, drepper@cygnus.com.
8 * ====================================================
9 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
11 * Developed at SunPro, a Sun Microsystems, Inc. business.
12 * Permission to use, copy, modify, and distribute this
13 * software is freely granted, provided that this notice
15 * ====================================================
18 #if defined(LIBM_SCCS) && !defined(lint)
19 static char rcsid
[] = "$NetBSD: $";
24 * return the next machine floating-point number of x in the
30 #include <math_private.h>
33 long double __nextafterl(long double x
, long double y
)
35 long double __nextafterl(x
,y
)
39 u_int32_t hx
,hy
,ix
,iy
;
43 GET_LDOUBLE_WORDS(esx
,hx
,lx
,x
);
44 GET_LDOUBLE_WORDS(esy
,hy
,ly
,y
);
45 ix
= esx
&0x7fff; /* |x| */
46 iy
= esy
&0x7fff; /* |y| */
48 /* Intel's extended format has the normally implicit 1 explicit
50 if(((ix
==0x7fff)&&(((hx
&0x7fffffff)|lx
)!=0)) || /* x is nan */
51 ((iy
==0x7fff)&&(((hy
&0x7fffffff)|ly
)!=0))) /* y is nan */
53 if(x
==y
) return y
; /* x=y, return y */
54 if((ix
|hx
|lx
)==0) { /* x == 0 */
56 SET_LDOUBLE_WORDS(x
,esy
&0x8000,0,1);/* return +-minsubnormal */
57 u
= math_opt_barrier (x
);
59 math_force_eval (u
); /* raise underflow flag */
62 if(esx
>=0) { /* x > 0 */
63 if(esx
>esy
||((esx
==esy
) && (hx
>hy
||((hx
==hy
)&&(lx
>ly
))))) {
66 if (hx
<= 0x80000000) {
79 } else { /* x < y, x += ulp */
83 if (hx
==0 || (esx
== 0 && hx
== 0x80000000)) {
90 if(esy
>=0||(esx
>esy
||((esx
==esy
)&&(hx
>hy
||((hx
==hy
)&&(lx
>ly
)))))){
93 if (hx
<= 0x80000000) {
102 } else { /* x > y, x += ulp */
106 if (hx
==0 || (esx
== 0xffff8000 && hx
== 0x80000000)) {
114 if(esy
==0x7fff) return x
+x
; /* overflow */
116 long double u
= x
*x
; /* underflow */
117 math_force_eval (u
); /* raise underflow flag */
119 SET_LDOUBLE_WORDS(x
,esx
,hx
,lx
);
122 weak_alias (__nextafterl
, nextafterl
)
123 strong_alias (__nextafterl
, __nexttowardl
)
124 weak_alias (__nextafterl
, nexttowardl
)