#include "libm.h" | |
// FIXME | |
/* asin(z) = -i log(i z + sqrt(1 - z*z)) */ | |
double complex casin(double complex z) | |
{ | |
double complex w; | |
double x, y; | |
x = creal(z); | |
y = cimag(z); | |
w = CMPLX(1.0 - (x - y)*(x + y), -2.0*x*y); | |
return clog(CMPLX(-y, x) + csqrt(w)); | |
} |