diff options
Diffstat (limited to 'sysdeps/ieee754/ldbl-128/s_atanl.c')
| -rw-r--r-- | sysdeps/ieee754/ldbl-128/s_atanl.c | 7 |
1 files changed, 7 insertions, 0 deletions
diff --git a/sysdeps/ieee754/ldbl-128/s_atanl.c b/sysdeps/ieee754/ldbl-128/s_atanl.c index dc5e7ef462..1367b6b15d 100644 --- a/sysdeps/ieee754/ldbl-128/s_atanl.c +++ b/sysdeps/ieee754/ldbl-128/s_atanl.c @@ -59,6 +59,8 @@ <http://www.gnu.org/licenses/>. */ +#include <float.h> +#include <math.h> #include <math_private.h> /* arctan(k/8), k = 0, ..., 82 */ @@ -200,6 +202,11 @@ __atanl (long double x) if (k <= 0x3fc50000) /* |x| < 2**-58 */ { + if (fabsl (x) < LDBL_MIN) + { + long double force_underflow = x * x; + math_force_eval (force_underflow); + } /* Raise inexact. */ if (huge + x > 0.0) return x; |