From b02eed9c4841913d690a2d0029737d72615384fe Mon Sep 17 00:00:00 2001 From: Szabolcs Nagy Date: Tue, 21 Apr 2015 23:04:11 +0000 Subject: complex: fix ctanh(+-0+i*nan) and ctanh(+-0+-i*inf) These cases were incorrect in C11 as described by http://www.open-std.org/jtc1/sc22/wg14/www/docs/n1886.htm --- src/complex/ctanh.c | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) (limited to 'src/complex/ctanh.c') diff --git a/src/complex/ctanh.c b/src/complex/ctanh.c index 0461050d..3ba3a899 100644 --- a/src/complex/ctanh.c +++ b/src/complex/ctanh.c @@ -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 -- cgit v1.2.1