* mingwex/complex/(cabsl.c cacosl.c cacoshl.c cargl.c casinl.c

	casinhl.c catanl.c catanhl.c ccosl.c ccoshl.c cexpl.c cimagl.c
	clogl.c cpowl.c cprojl.c creall.c csinl.c csinhl.c csqrtl.c
	ctanl.c ctanhl.c): New files.
	* mingwex/Makefile.in (COMPLEX_DISTFILES): Adjust.
	(COMPLEX_OBJS): Adjust.
	* include/complex.h (cabsl, cacosl, cacoshl, cargl, casinf.
	casinhl, catanl, catanhl, ccosl, ccoshl, cexpl, cimagl, clogl,
	cpowl, cprojl, creall, csinl, csinhl, csqrtl, ctanl, ctanhl):
	Declare.

1 files changed, 48 insertions, 0 deletions

diff --git a/winsup/mingw/mingwex/complex/casinl.c b/winsup/mingw/mingwex/complex/casinl.c
new file mode 100755
index 000000000..3093fa542
--- /dev/null
+++ b/winsup/mingw/mingwex/complex/casinl.c
@@ -0,0 +1,48 @@
+/*
+   casinl.c
+   Contributed by Danny Smith
+   2005-01-04
+*/
+
+#include <math.h>
+#include <complex.h>
+
+/* casin (Z ) = -I * clog(I * Z + csqrt (1.0 - Z * Z))) */
+
+long double complex casinl (long double complex Z)
+{
+  long double complex Res;
+  long double x, y;
+
+  x = __real__ Z;
+  y = __imag__ Z;
+
+  if (y == 0.0L)
+    {
+      __real__ Res = asinl (x);
+      __imag__ Res = 0.0L;
+    }
+  else  /* -I * clog(I * Z + csqrt(1.0 - Z * Z))) */
+    {
+      long double complex ZZ;
+                                 
+      /* Z * Z = ((x - y) * (x + y)) + (2.0 * x * y) * I */
+      /* calculate 1 - Z * Z */
+      __real__ ZZ = 1.0L - (x - y) * (x + y);
+      __imag__ ZZ = -2.0L * x * y;
+      ZZ = csqrtl (ZZ);
+
+
+      /* add  I * Z  to ZZ */
+
+      __real__ ZZ -= y;
+      __imag__ ZZ += x;
+
+      ZZ = clogl (ZZ);
+
+      /* mult by -I */
+      __real__ Res = __imag__ ZZ;
+      __imag__ Res = - __real__ ZZ;  
+    }
+  return (Res);
+}