| #include "libm.h" | |
| #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 | |
| long double rintl(long double x) | |
| { | |
| return rint(x); | |
| } | |
| #elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 | |
| static const long double toint = 1/LDBL_EPSILON; | |
| long double rintl(long double x) | |
| { | |
| union ldshape u = {x}; | |
| int e = u.i.se & 0x7fff; | |
| int s = u.i.se >> 15; | |
| long double y; | |
| if (e >= 0x3fff+LDBL_MANT_DIG-1) | |
| return x; | |
| if (s) | |
| y = x - toint + toint; | |
| else | |
| y = x + toint - toint; | |
| if (y == 0) | |
| return 0*x; | |
| return y; | |
| } | |
| #endif |