diff options
Diffstat (limited to 'src/complex/catanh.c')
| -rw-r--r-- | src/complex/catanh.c | 57 |
1 files changed, 56 insertions, 1 deletions
diff --git a/src/complex/catanh.c b/src/complex/catanh.c index 2be96dfc..198fc5ea 100644 --- a/src/complex/catanh.c +++ b/src/complex/catanh.c @@ -1,10 +1,65 @@ # define TGSOURCE "../complex/catanh.c" #include "_tgmath.h" - #include <complex.h> +#include "math.h" +#include "fenv.h" + +#include "../math/M_PI_2.c" TYPE complex TGFN(catanh)(TYPE complex z) { + int classr = fpclassify(TGFN(creal)(z)); + int classi = fpclassify(TGFN(cimag)(z)); + int signr = signbit(TGFN(creal)(z)); + int signi = signbit(TGFN(cimag)(z)); + + if (classr == FP_ZERO && classi == FP_ZERO) { + return TGCMPLX(0.0, 0.0); + } + + if (classr == FP_ZERO && classi == FP_NAN) { + return TGCMPLX(0.0, NAN); + } + + if (TGFN(creal)(z) == 1.0 && classi == FP_ZERO) { + feraiseexcept(FE_DIVBYZERO); + return TGCMPLX(INFINITY, 0.0); + } + + if (classr != FP_INFINITE && !signr && classi == FP_INFINITE) { + return TGCMPLX(0.0, M_PI_2); + } + + if (classr != FP_INFINITE && !signr && classi == FP_NAN) { + feraiseexcept(FE_INVALID); + return TGCMPLX(NAN, NAN); + } + + if (classr == FP_INFINITE && classi != FP_INFINITE && !signi) { + return TGCMPLX(0.0, M_PI_2); + } + + if (classr == FP_INFINITE && classi == FP_INFINITE) { + return TGCMPLX(0.0, M_PI_2); + } + + if (classr == FP_INFINITE && classi == FP_NAN) { + return TGCMPLX(0.0, NAN); + } + + if (classr == FP_NAN && classi != FP_INFINITE) { + feraiseexcept(FE_INVALID); + return TGCMPLX(NAN, NAN); + } + + if (classr == FP_NAN && classi == FP_INFINITE) { + return TGCMPLX(0.0, M_PI_2); + } + + if (classr == FP_NAN && classi == FP_NAN) { + return TGCMPLX(NAN, NAN); + } + return z; } |