aboutsummaryrefslogtreecommitdiff
path: root/math/s_cexpf.c
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2013-08-23 19:45:38 +0000
committerJoseph Myers <joseph@codesourcery.com>2013-08-23 19:45:38 +0000
commit8fe89494e6516048759425ec30d8878a6233e00f (patch)
treeb3f65b97d02bac4f183da5fe51e01cad8b8a3780 /math/s_cexpf.c
parent936241e4b2ec90bbb97d1b37bc78726372ec966f (diff)
downloadglibc-8fe89494e6516048759425ec30d8878a6233e00f.zip
glibc-8fe89494e6516048759425ec30d8878a6233e00f.tar.gz
glibc-8fe89494e6516048759425ec30d8878a6233e00f.tar.bz2
Fix cexp (NaN + i0) (bug 15532).
Diffstat (limited to 'math/s_cexpf.c')
-rw-r--r--math/s_cexpf.c14
1 files changed, 10 insertions, 4 deletions
diff --git a/math/s_cexpf.c b/math/s_cexpf.c
index fa942d3..7c42205 100644
--- a/math/s_cexpf.c
+++ b/math/s_cexpf.c
@@ -145,12 +145,18 @@ __cexpf (__complex__ float x)
}
else
{
- /* If the real part is NaN the result is NaN + iNaN. */
+ /* If the real part is NaN the result is NaN + iNaN unless the
+ imaginary part is zero. */
__real__ retval = __nanf ("");
- __imag__ retval = __nanf ("");
+ if (icls == FP_ZERO)
+ __imag__ retval = __imag__ x;
+ else
+ {
+ __imag__ retval = __nanf ("");
- if (rcls != FP_NAN || icls != FP_NAN)
- feraiseexcept (FE_INVALID);
+ if (rcls != FP_NAN || icls != FP_NAN)
+ feraiseexcept (FE_INVALID);
+ }
}
return retval;