1#include "complex_impl.h"
2
3// FIXME
4
5float 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}