#ifndef __NO_MATH_INLINES # define __NO_MATH_INLINES 1 /* for glibc */ #endif #include #include /* Override the values of , like done in float.in.h. */ #if defined __i386__ && (defined __BEOS__ || defined __OpenBSD__) # undef LDBL_MANT_DIG # define LDBL_MANT_DIG 64 # undef LDBL_MIN_EXP # define LDBL_MIN_EXP (-16381) # undef LDBL_MAX_EXP # define LDBL_MAX_EXP 16384 #endif #if defined __i386__ && (defined __FreeBSD__ || defined __DragonFly__) # undef LDBL_MANT_DIG # define LDBL_MANT_DIG 64 # undef LDBL_MIN_EXP # define LDBL_MIN_EXP (-16381) # undef LDBL_MAX_EXP # define LDBL_MAX_EXP 16384 #endif #if (defined _ARCH_PPC || defined _POWER) && defined _AIX && (LDBL_MANT_DIG == 106) && defined __GNUC__ # undef LDBL_MIN_EXP # define LDBL_MIN_EXP DBL_MIN_EXP #endif #if defined __sgi && (LDBL_MANT_DIG >= 106) # undef LDBL_MANT_DIG # define LDBL_MANT_DIG 106 # if defined __GNUC__ # undef LDBL_MIN_EXP # define LDBL_MIN_EXP DBL_MIN_EXP # endif #endif #undef exp2l extern #ifdef __cplusplus "C" #endif long double exp2l (long double); static long double dummy (long double x) { return 0; } static long double zero; int main (int argc, char *argv[]) { long double (*my_exp2l) (long double) = argc ? exp2l : dummy; int result = 0; /* This test fails on NetBSD 8.0. */ { const long double TWO_LDBL_MANT_DIG = /* 2^LDBL_MANT_DIG */ (long double) (1U << ((LDBL_MANT_DIG - 1) / 5)) * (long double) (1U << ((LDBL_MANT_DIG - 1 + 1) / 5)) * (long double) (1U << ((LDBL_MANT_DIG - 1 + 2) / 5)) * (long double) (1U << ((LDBL_MANT_DIG - 1 + 3) / 5)) * (long double) (1U << ((LDBL_MANT_DIG - 1 + 4) / 5)); long double x = 11.358L; long double err = (my_exp2l (x) * my_exp2l (- x) - 1.0L) * TWO_LDBL_MANT_DIG; if (!(err >= -100.0L && err <= 100.0L)) result |= 1; } return result; }