* mingwex/complex/(cabsf.c cacosf.c cacoshf.c cargf.c casinf.c

	casinhf.c catanf.c catanhf.c ccosf.c ccoshf.c cexpf.c cimagf.c
	clogf.c cpowf.c cprojf.c crealf.c csinf.c csinhf.c csqrtf.c
	ctanf.c ctanhf.c): New files.
	* mingwex/Makefile.in (COMPLEX_DISTFILES): Adjust.
	(COMPLEX_OBJS(: Adjust.
	* include/complex.h (cabsf, cacosf, cacoshf, cargf, casinf.
	casinhf, catanf, catanhf, ccosf, ccoshf, cexpf, cimagf, clogf,
	cpowf, cprojf, crealf, csinf, csinhf, csqrtf, ctanf, ctanhf):
	Declare.

1 files changed, 43 insertions, 0 deletions

diff --git a/winsup/mingw/mingwex/complex/cpowf.c b/winsup/mingw/mingwex/complex/cpowf.c
new file mode 100755
index 000000000..eeeed1b7f
--- /dev/null
+++ b/winsup/mingw/mingwex/complex/cpowf.c
@@ -0,0 +1,43 @@
+/*  cpowf.c */
+/*
+   Contributed by Danny Smith
+   2004-12-24
+*/
+
+/* cpow(X, Y) = cexp(X * clog(Y)) */
+
+#include <math.h>
+#include <complex.h>
+
+float complex cpowf (float complex X, float complex Y)
+{
+  float complex Res;
+  float i;
+  float r = _hypot (__real__ X, __imag__ X);
+  if (r == 0.0f)
+    {
+       __real__ Res = __imag__ Res = 0.0;
+    }
+  else
+    {
+      float rho;
+      float theta;
+      i = cargf (X);
+      theta = i * __real__ Y;
+ 
+      if (__imag__ Y == 0.0f)
+	/* This gives slightly more accurate results in these cases. */
+   	rho = powf (r, __real__ Y);
+      else
+	{
+          r = logf (r);
+	  /* rearrangement of cexp(X * clog(Y)) */
+	  theta += r * __imag__ Y;
+	  rho = expf (r * __real__ Y - i * __imag__ Y);
+	}
+
+      __real__ Res = rho * cosf (theta);
+      __imag__ Res = rho * sinf (theta);
+    }
+  return  Res;
+}