3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
4 long double nextafterl(long double x, long double y)
6 return nextafter(x, y);
8 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
9 #define MSB ((uint64_t)1<<63)
10 long double nextafterl(long double x, long double y)
14 if (isnan(x) || isnan(y))
22 ux.bits.sign = uy.bits.sign;
23 } else if (x < y ^ ux.bits.sign) {
25 if ((ux.bits.m & ~MSB) == 0) {
30 if ((ux.bits.m & ~MSB) == 0) {
37 /* raise overflow if ux.value is infinite and x is finite */
38 if (ux.bits.exp == 0x7fff)
40 /* raise underflow if ux.value is subnormal or zero */
42 FORCE_EVAL(x*x + ux.value*ux.value);
45 #elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
46 long double nextafterl(long double x, long double y)
50 if (isnan(x) || isnan(y))
58 ux.bits.sign = uy.bits.sign;
59 } else if (x < y ^ ux.bits.sign) {
61 if (ux.bits.mlo == 0) {
67 if (ux.bits.mlo == 0) {
74 /* raise overflow if ux.value is infinite and x is finite */
75 if (ux.bits.exp == 0x7fff)
77 /* raise underflow if ux.value is subnormal or zero */
79 FORCE_EVAL(x*x + ux.value*ux.value);