diff options
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 c9ed719..08c80fa 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 36157ff..655e4e8 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 364be8a..fa942d3 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 1bfce78..d827bc3 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 { |