1 /* s_nextafterl.c -- long double version of s_nextafter.c.
2 * Conversion to long double by Ulrich Drepper,
3 * Cygnus Support, drepper@cygnus.com.
4 * Fixed for m68k by Andreas Schwab <schwab@suse.de>.
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>
32 long double __nextafterl(long double x
, long double y
)
34 int32_t ix
,iy
,esx
,esy
;
35 u_int32_t hx
,hy
,lx
,ly
;
37 GET_LDOUBLE_WORDS(esx
,hx
,lx
,x
);
38 GET_LDOUBLE_WORDS(esy
,hy
,ly
,y
);
39 ix
= esx
&0x7fff; /* |x| */
40 iy
= esy
&0x7fff; /* |y| */
42 if(((ix
==0x7fff)&&((hx
&0x7fffffff)|lx
)!=0) || /* x is nan */
43 ((iy
==0x7fff)&&((hy
&0x7fffffff)|ly
)!=0)) /* y is nan */
45 if(x
==y
) return y
; /* x=y, return y */
46 if((ix
|hx
|lx
)==0) { /* x == 0 */
47 SET_LDOUBLE_WORDS(x
,esy
&0x8000,0,1);/* return +-minsubnormal */
49 if(y
==x
) return y
; else return x
; /* raise underflow flag */
51 if(esx
>=0) { /* x > 0 */
52 if(esx
>esy
||((esx
==esy
) && (hx
>hy
||((hx
==hy
)&&(lx
>ly
))))) {
55 if (ix
!= 0 && hx
== 0x80000000) hx
= 0;
60 } else { /* x < y, x += ulp */
71 if(esy
>=0||esx
>esy
||((esx
==esy
) && (hx
>hy
||((hx
==hy
)&&(lx
>ly
))))){
74 if (ix
!= 0 && hx
== 0x80000000) hx
= 0;
79 } else { /* x > y, x += ulp */
91 if(esy
==0x7fff) return x
+x
; /* overflow */
92 if(esy
==0 && (hx
& 0x80000000) == 0) { /* underflow */
94 if(y
!=x
) { /* raise underflow flag */
95 SET_LDOUBLE_WORDS(y
,esx
,hx
,lx
);
99 SET_LDOUBLE_WORDS(x
,esx
,hx
,lx
);
102 weak_alias (__nextafterl
, nextafterl
)
103 strong_alias (__nextafterl
, __nexttowardl
)
104 weak_alias (__nextafterl
, nexttowardl
)