summaryrefslogtreecommitdiff
path: root/sysdeps/generic/s_ctanf.c
diff options
context:
space:
mode:
Diffstat (limited to 'sysdeps/generic/s_ctanf.c')
-rw-r--r--sysdeps/generic/s_ctanf.c17
1 files changed, 14 insertions, 3 deletions
diff --git a/sysdeps/generic/s_ctanf.c b/sysdeps/generic/s_ctanf.c
index e02971dd44..58d9d13298 100644
--- a/sysdeps/generic/s_ctanf.c
+++ b/sysdeps/generic/s_ctanf.c
@@ -1,5 +1,5 @@
/* Complex tangent function for float.
- Copyright (C) 1997 Free Software Foundation, Inc.
+ Copyright (C) 1997, 2005 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@@ -61,8 +61,19 @@ __ctanf (__complex__ float x)
den = cos2rx + __ieee754_coshf (2.0 * __imag__ x);
- __real__ res = sin2rx / den;
- __imag__ res = __ieee754_sinhf (2.0 * __imag__ x) / den;
+
+ if (den == 0.0)
+ {
+ __complex__ float ez = __cexpf (1.0i * x);
+ __complex__ float emz = __cexpf (-1.0i * x);
+
+ res = (ez - emz) / (ez + emz) * -1.0i;
+ }
+ else
+ {
+ __real__ res = sin2rx / den;
+ __imag__ res = __ieee754_sinhf (2.0 * __imag__ x) / den;
+ }
}
return res;