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