src/complex/casinf.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23

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


#ifdef TEST
int
main(void) {
	return 0;
}
#endif