From: dannysmith <dannysmith>
Date: Wed, 12 Oct 2005 06:46:17 +0000 (+0000)
Subject: 	* mingwex/complex/csqrt.c (csqrt): The sign of real part
X-Git-Url: http://git.osdn.net/view?a=commitdiff_plain;h=3acd35cfe2b1ea5449a123623885efdc302b6390;p=pf3gnuchains%2Fpf3gnuchains3x.git

	* 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.
---

diff --git a/winsup/mingw/ChangeLog b/winsup/mingw/ChangeLog
index 2cdd5c4d29..896c5b1d0b 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 3717939f40..b5f8868e94 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 b39cbb4689..7c37e99ce6 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 f058b216a0..1b2ebbe5be 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