1#include <limits.h> 2#include <fenv.h> 3#include "libm.h" 4 5 6#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 7long long llrintl(long double x) 8{ 9 return llrint(x); 10} 11#elif defined(FE_INEXACT) 12/* 13see comments in lrint.c 14 15Note that if LLONG_MAX == 0x7fffffffffffffff && LDBL_MANT_DIG == 64 16then x == 2**63 - 0.5 is the only input that overflows and 17raises inexact (with tonearest or upward rounding mode) 18*/ 19long long llrintl(long double x) 20{ 21 #pragma STDC FENV_ACCESS ON 22 int e; 23 24 e = fetestexcept(FE_INEXACT); 25 x = rintl(x); 26 if (!e && (x > LLONG_MAX || x < LLONG_MIN)) 27 feclearexcept(FE_INEXACT); 28 /* conversion */ 29 return x; 30} 31#else 32long long llrintl(long double x) 33{ 34 return rintl(x); 35} 36#endif 37