diff options
Diffstat (limited to 'math/k_casinh.c')
| -rw-r--r-- | math/k_casinh.c | 20 |
1 files changed, 20 insertions, 0 deletions
diff --git a/math/k_casinh.c b/math/k_casinh.c index 7f98f24a80..41cd5ec47d 100644 --- a/math/k_casinh.c +++ b/math/k_casinh.c @@ -57,6 +57,26 @@ __kernel_casinh (__complex__ double x, int adj) res = __clog (y); __real__ res += M_LN2; } + else if (rx >= 0.5 && ix < DBL_EPSILON / 8.0) + { + double s = __ieee754_hypot (1.0, rx); + + __real__ res = __ieee754_log (rx + s); + if (adj) + __imag__ res = __ieee754_atan2 (s, __imag__ x); + else + __imag__ res = __ieee754_atan2 (ix, s); + } + else if (rx < DBL_EPSILON / 8.0 && ix >= 1.5) + { + double s = __ieee754_sqrt ((ix + 1.0) * (ix - 1.0)); + + __real__ res = __ieee754_log (ix + s); + if (adj) + __imag__ res = __ieee754_atan2 (rx, __copysign (s, __imag__ x)); + else + __imag__ res = __ieee754_atan2 (s, rx); + } else { __real__ y = (rx - ix) * (rx + ix) + 1.0; |