+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.
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;
}
}
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
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