From 05e166c887612250d461b5fe7c0f0064cd1a0c41 Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Wed, 3 Apr 2013 16:51:46 +0000 Subject: Fix missing underflow from cexp (bug 14478). --- math/libm-test.inc | 3 +++ math/s_cexp.c | 12 ++++++++++++ math/s_cexpf.c | 12 ++++++++++++ math/s_cexpl.c | 12 ++++++++++++ 4 files changed, 39 insertions(+) (limited to 'math') diff --git a/math/libm-test.inc b/math/libm-test.inc index c9ed719d18..08c80fad42 100644 --- a/math/libm-test.inc +++ b/math/libm-test.inc @@ -4794,6 +4794,9 @@ cexp_test (void) TEST_c_c (cexp, 1e6, min_value, plus_infty, plus_infty, OVERFLOW_EXCEPTION); TEST_c_c (cexp, 1e6, -min_value, plus_infty, minus_infty, OVERFLOW_EXCEPTION); + TEST_c_c (cexp, min_value, min_subnorm_value, 1.0, min_subnorm_value, UNDERFLOW_EXCEPTION); + TEST_c_c (cexp, min_value, -min_subnorm_value, 1.0, -min_subnorm_value, UNDERFLOW_EXCEPTION); + END (cexp, complex); } diff --git a/math/s_cexp.c b/math/s_cexp.c index 36157ff61b..655e4e8dee 100644 --- a/math/s_cexp.c +++ b/math/s_cexp.c @@ -74,6 +74,18 @@ __cexp (__complex__ double x) __real__ retval = exp_val * cosix; __imag__ retval = exp_val * sinix; } + if (fabs (__real__ retval) < DBL_MIN) + { + volatile double force_underflow + = __real__ retval * __real__ retval; + (void) force_underflow; + } + if (fabs (__imag__ retval) < DBL_MIN) + { + volatile double force_underflow + = __imag__ retval * __imag__ retval; + (void) force_underflow; + } } else { diff --git a/math/s_cexpf.c b/math/s_cexpf.c index 364be8ac31..fa942d34f7 100644 --- a/math/s_cexpf.c +++ b/math/s_cexpf.c @@ -74,6 +74,18 @@ __cexpf (__complex__ float x) __real__ retval = exp_val * cosix; __imag__ retval = exp_val * sinix; } + if (fabsf (__real__ retval) < FLT_MIN) + { + volatile float force_underflow + = __real__ retval * __real__ retval; + (void) force_underflow; + } + if (fabsf (__imag__ retval) < FLT_MIN) + { + volatile float force_underflow + = __imag__ retval * __imag__ retval; + (void) force_underflow; + } } else { diff --git a/math/s_cexpl.c b/math/s_cexpl.c index 1bfce78aeb..d827bc3de4 100644 --- a/math/s_cexpl.c +++ b/math/s_cexpl.c @@ -74,6 +74,18 @@ __cexpl (__complex__ long double x) __real__ retval = exp_val * cosix; __imag__ retval = exp_val * sinix; } + if (fabsl (__real__ retval) < LDBL_MIN) + { + volatile long double force_underflow + = __real__ retval * __real__ retval; + (void) force_underflow; + } + if (fabsl (__imag__ retval) < LDBL_MIN) + { + volatile long double force_underflow + = __imag__ retval * __imag__ retval; + (void) force_underflow; + } } else { -- cgit v1.2.3-70-g09d2