5 #if FLT_EVAL_METHOD==0 || FLT_EVAL_METHOD==1
6 #define EPS DBL_EPSILON
7 #elif FLT_EVAL_METHOD==2
8 #define EPS LDBL_EPSILON
10 static const double_t toint = 1/EPS;
14 union {double f; uint64_t i;} u = {x};
15 int e = u.i>>52 & 0x7ff;
22 y = x - toint + toint;
24 y = x + toint - toint;