#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)); | |
} |