Fix little oversight.
authordrepper <drepper>
Wed, 29 Jul 1998 14:51:35 +0000 (14:51 +0000)
committerdrepper <drepper>
Wed, 29 Jul 1998 14:51:35 +0000 (14:51 +0000)
sysdeps/libm-ieee754/s_csqrt.c
sysdeps/libm-ieee754/s_csqrtf.c
sysdeps/libm-ieee754/s_csqrtl.c

index ae106bd..a4dedf8 100644 (file)
@@ -84,7 +84,6 @@ __csqrt (__complex__ double x)
        }
       else
        {
-#if 0
          double d, r, s;
 
          d = __ieee754_hypot (__real__ x, __imag__ x);
@@ -98,20 +97,11 @@ __csqrt (__complex__ double x)
          else
            {
              s = __ieee754_sqrt (0.5 * d - 0.5 * __real__ x);
-             r = (0.5 * __imag__ x) / s;
+             r = fabs ((0.5 * __imag__ x) / s);
            }
 
          __real__ res = r;
          __imag__ res = __copysign (s, __imag__ x);
-#else
-         double d, imag;
-
-         d = __ieee754_hypot (__real__ x, __imag__ x);
-         imag = __ieee754_sqrt (0.5 * (d - __real__ x));
-
-         __real__ res = __ieee754_sqrt (0.5 * (d + __real__ x));
-         __imag__ res = __copysign (imag, __imag__ x);
-#endif
        }
     }
 
index 015b0cd..3bc979f 100644 (file)
@@ -84,7 +84,6 @@ __csqrtf (__complex__ float x)
        }
       else
        {
-#if 0
          float d, r, s;
 
          d = __ieee754_hypotf (__real__ x, __imag__ x);
@@ -98,20 +97,11 @@ __csqrtf (__complex__ float x)
          else
            {
              s = __ieee754_sqrtf (0.5f * d - 0.5f * __real__ x);
-             r = (0.5f * __imag__ x) / s;
+             r = fabsf ((0.5f * __imag__ x) / s);
            }
 
          __real__ res = r;
          __imag__ res = __copysignf (s, __imag__ x);
-#else
-         float d, imag;
-
-         d = __ieee754_hypotf (__real__ x, __imag__ x);
-         imag = __ieee754_sqrtf (0.5 * (d - __real__ x));
-
-         __real__ res = __ieee754_sqrtf (0.5 * (d + __real__ x));
-         __imag__ res = __copysignf (imag, __imag__ x);
-#endif
        }
     }
 
index 4b7ed98..b806bc6 100644 (file)
@@ -84,7 +84,6 @@ __csqrtl (__complex__ long double x)
        }
       else
        {
-#if 0
          long double d, r, s;
 
          d = __ieee754_hypotl (__real__ x, __imag__ x);
@@ -98,20 +97,11 @@ __csqrtl (__complex__ long double x)
          else
            {
              s = __ieee754_sqrtl (0.5L * d - 0.5L * __real__ x);
-             r = (0.5L * __imag__ x) / s;
+             r = fabsl ((0.5L * __imag__ x) / s);
            }
 
          __real__ res = r;
          __imag__ res = __copysignl (s, __imag__ x);
-#else
-         long double d, imag;
-
-         d = __ieee754_hypotl (__real__ x, __imag__ x);
-         imag = __ieee754_sqrtl (0.5 * (d - __real__ x));
-
-         __real__ res = __ieee754_sqrtl (0.5 * (d + __real__ x));
-         __imag__ res = __copysignl (imag, __imag__ x);
-#endif
        }
     }