| #include "libm.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)); | |
| } |