diff options
| author | Joseph Myers <joseph@codesourcery.com> | 2013-04-03 16:51:46 +0000 |
|---|---|---|
| committer | Joseph Myers <joseph@codesourcery.com> | 2013-04-03 16:51:46 +0000 |
| commit | 05e166c887612250d461b5fe7c0f0064cd1a0c41 (patch) | |
| tree | 9cb5d95b1af53bfc90b7bf1bf0a6772336df00d3 /math | |
| parent | 5b535ac41945fa369a5f641acdeb0d7eadf4668e (diff) | |
| download | glibc-05e166c887612250d461b5fe7c0f0064cd1a0c41.tar.gz glibc-05e166c887612250d461b5fe7c0f0064cd1a0c41.tar.xz glibc-05e166c887612250d461b5fe7c0f0064cd1a0c41.zip | |
Fix missing underflow from cexp (bug 14478).
Diffstat (limited to 'math')
| -rw-r--r-- | math/libm-test.inc | 3 | ||||
| -rw-r--r-- | math/s_cexp.c | 12 | ||||
| -rw-r--r-- | math/s_cexpf.c | 12 | ||||
| -rw-r--r-- | math/s_cexpl.c | 12 |
4 files changed, 39 insertions, 0 deletions
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 { |