3 float nexttowardf(float x, long double y)
5 union {float f; uint32_t i;} ux = {x};
8 if (isnan(x) || isnan(y))
27 e = ux.i & 0x7f800000;
28 /* raise overflow if ux.f is infinite and x is finite */
31 /* raise underflow if ux.f is subnormal or zero */
33 FORCE_EVAL(x*x + ux.f*ux.f);