update ctype data to unicode 12.1.0
[musl] / src / complex / ctanh.c
index 0461050..54004cd 100644 (file)
@@ -63,7 +63,7 @@
  *   precision.  I also handle large x differently.
  */
 
-#include "libm.h"
+#include "complex_impl.h"
 
 double complex ctanh(double complex z)
 {
@@ -101,11 +101,13 @@ double complex ctanh(double complex z)
        }
 
        /*
+        * ctanh(+-0 + i NAN) = +-0 + i NaN
+        * ctanh(+-0 +- i Inf) = +-0 + i NaN
         * ctanh(x + i NAN) = NaN + i NaN
         * ctanh(x +- i Inf) = NaN + i NaN
         */
        if (!isfinite(y))
-               return CMPLX(y - y, y - y);
+               return CMPLX(x ? y - y : x, y - y);
 
        /*
         * ctanh(+-huge + i +-y) ~= +-1 +- i 2sin(2y)/exp(2x), using the