3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
4 double nexttoward(double x
, long double y
)
6 return nextafter(x
, y
);
9 double nexttoward(double x
, long double y
)
11 union {double f
; uint64_t i
;} ux
= {x
};
14 if (isnan(x
) || isnan(y
))
34 /* raise overflow if ux.f is infinite and x is finite */
37 /* raise underflow if ux.f is subnormal or zero */
39 FORCE_EVAL(x
*x
+ ux
.f
*ux
.f
);