1 #include "complex_impl.h"
2 
3 // FIXME
4 
casinf(float complex z)5 float complex casinf(float complex z)
6 {
7 	float complex w;
8 	float x, y;
9 
10 	x = crealf(z);
11 	y = cimagf(z);
12 	w = CMPLXF(1.0 - (x - y)*(x + y), -2.0*x*y);
13 	float complex r = clogf(CMPLXF(-y, x) + csqrtf(w));
14 	return CMPLXF(cimagf(r), -crealf(r));
15 }
16