* mingwex/complex/csqrt.c (csqrt): The sign of real part

	of result is positive when real part of arg == 0;
	* mingwex/complex/csqrtf.c (csqrtf): Ditto.
	* mingwex/complex/csqrtl.c (csqrtl): Ditto.

4 files changed, 16 insertions, 8 deletions

diff --git a/winsup/mingw/ChangeLog b/winsup/mingw/ChangeLog
index 2cdd5c4d2..896c5b1d0 100644
--- a/winsup/mingw/ChangeLog
+++ b/winsup/mingw/ChangeLog
@@ -1,5 +1,12 @@
 2005-10-12  Danny Smith  <dannysmith@users.sourceforge.net>
 
+	* mingwex/complex/csqrt.c (csqrt): The sign of real part
+	of result is positive when real part of arg == 0;
+	* mingwex/complex/csqrtf.c (csqrtf): Ditto.
+	* mingwex/complex/csqrtl.c (csqrtl): Ditto.
+
+2005-10-12  Danny Smith  <dannysmith@users.sourceforge.net>
+
 	* include/time.h (_time64): Correct prototype.
 
 2005-10-08  Danny Smith  <dannysmith@users.sourceforge.net>
diff --git a/winsup/mingw/mingwex/complex/csqrt.c b/winsup/mingw/mingwex/complex/csqrt.c
index 3717939f4..b5f8868e9 100644
--- a/winsup/mingw/mingwex/complex/csqrt.c
+++ b/winsup/mingw/mingwex/complex/csqrt.c
@@ -31,22 +31,23 @@ double complex  csqrt (double complex Z)
   else if (x == 0.0)
     {
       t = sqrt(0.5 * fabs (y));
-      __real__ Res = y > 0 ? t : -t;
-      __imag__ Res = t;
+      __real__ Res = t;
+      __imag__ Res = y > 0 ? t : -t;
     }
 
   else
     {
       t = sqrt (2.0  * (_hypot (x, y) + fabs (x)));
+      double u = t / 2.0;
       if ( x > 0.0)
         {	
-          __real__ Res = 0.5 * t;
+          __real__ Res = u;
           __imag__ Res = y / t;
         }
       else
         {
 	  __real__ Res = fabs ( y / t);
-	  __imag__ Res = (y < 0.0 ? -0.5 : 0.5) * t;
+	  __imag__ Res = y < 0.0 ? -u : u;
         }
     }
 
diff --git a/winsup/mingw/mingwex/complex/csqrtf.c b/winsup/mingw/mingwex/complex/csqrtf.c
index b39cbb468..7c37e99ce 100755
--- a/winsup/mingw/mingwex/complex/csqrtf.c
+++ b/winsup/mingw/mingwex/complex/csqrtf.c
@@ -25,8 +25,8 @@ float complex  csqrtf (float complex Z)
   else if (x == 0.0f)
     {
       r = sqrtf(0.5f * fabsf (y));
-      __real__ Res = y > 0 ? r : -r;
-      __imag__ Res = r;
+      __real__ Res = r;
+      __imag__ Res = y > 0 ? r : -r;
     }
 
   else
diff --git a/winsup/mingw/mingwex/complex/csqrtl.c b/winsup/mingw/mingwex/complex/csqrtl.c
index f058b216a..1b2ebbe5b 100755
--- a/winsup/mingw/mingwex/complex/csqrtl.c
+++ b/winsup/mingw/mingwex/complex/csqrtl.c
@@ -31,8 +31,8 @@ long double complex  csqrtl (long double complex Z)
   else if (x == 0.0L)
     {
       r = sqrtl(0.5L * fabsl (y));
-      __real__ Res = y > 0 ? r : -r;
-      __imag__ Res = r;
+      __real__ Res = r;
+      __imag__ Res = y > 0 ? r : -r;
     }
 
   else