1 #include "libm.h"
2 
3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
rintl(long double x)4 long double rintl(long double x)
5 {
6 	return rint(x);
7 }
8 #elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
9 
10 static const long double toint = 1/LDBL_EPSILON;
11 
rintl(long double x)12 long double rintl(long double x)
13 {
14 	union ldshape u = {x};
15 	int e = u.i.se & 0x7fff;
16 	int s = u.i.se >> 15;
17 	long double y;
18 
19 	if (e >= 0x3fff+LDBL_MANT_DIG-1)
20 		return x;
21 	if (s)
22 		y = x - toint + toint;
23 	else
24 		y = x + toint - toint;
25 	if (y == 0)
26 		return 0*x;
27 	return y;
28 }
29 #endif
30