diff options
Diffstat (limited to 'src/complex/casinf.c')
-rw-r--r-- | src/complex/casinf.c | 15 |
1 files changed, 15 insertions, 0 deletions
diff --git a/src/complex/casinf.c b/src/complex/casinf.c new file mode 100644 index 0000000..2cda2f0 --- /dev/null +++ b/src/complex/casinf.c @@ -0,0 +1,15 @@ +#include "complex_impl.h" + +// FIXME + +float complex casinf(float complex z) +{ + float complex w; + float x, y; + + x = crealf(z); + y = cimagf(z); + w = CMPLXF(1.0 - (x - y)*(x + y), -2.0*x*y); + float complex r = clogf(CMPLXF(-y, x) + csqrtf(w)); + return CMPLXF(cimagf(r), -crealf(r)); +} |