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);