Initial version

1 files changed, 15 insertions, 0 deletions

diff --git a/src/complex/casinf.c b/src/complex/casinf.c
new file mode 100644
index 0000000..2cda2f0
--- /dev/null
+++ b/src/complex/casinf.c
@@ -0,0 +1,15 @@
+#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));
+}