fix incorrect results for catan with some inputs

[musl] / src / complex / catan.c

diff --git a/src/complex/catan.c b/src/complex/catan.c
index 39ce6cf..7dc2afe 100644 (file)
--- a/src/complex/catan.c
+++ b/src/complex/catan.c
@@ -91,29 +91,17 @@ double complex catan(double complex z)
        x = creal(z);
        y = cimag(z);
 
-       if (x == 0.0 && y > 1.0)
-               goto ovrf;
-
        x2 = x * x;
        a = 1.0 - x2 - (y * y);
-       if (a == 0.0)
-               goto ovrf;
 
        t = 0.5 * atan2(2.0 * x, a);
        w = _redupi(t);
 
        t = y - 1.0;
        a = x2 + (t * t);
-       if (a == 0.0)
-               goto ovrf;
 
        t = y + 1.0;
        a = (x2 + t * t)/a;
-       w = w + (0.25 * log(a)) * I;
-       return w;
-
-ovrf:
-       // FIXME
-       w = MAXNUM + MAXNUM * I;
+       w = CMPLX(w, 0.25 * log(a));
        return w;
 }