1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
#include <math.h> #include <float.h> #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 long double nearbyintl(long double x) { return nearbyint(x); } #else #include <fenv.h> long double nearbyintl(long double x) { #ifdef FE_INEXACT int e; e = fetestexcept(FE_INEXACT); #endif x = rintl(x); #ifdef FE_INEXACT if (!e) feclearexcept(FE_INEXACT); #endif return x; } #endif