1 #include "complex_impl.h"
5 /* asin(z) = -i log(i z + sqrt(1 - z*z)) */
7 double complex casin(double complex z)
14 w = CMPLX(1.0 - (x - y)*(x + y), -2.0*x*y);
15 double complex r = clog(CMPLX(-y, x) + csqrt(w));
16 return CMPLX(cimag(r), -creal(r));