if (fabsl (__imag__ x) < 2.0e-4 * fabsl (__real__ x)) t = 0.25 * __imag__ x * (__imag__ x / __real__ x);
d = __ieee754_hypotl (__real__ x, __imag__ x);  if (__real__ x > 0) { r = __ieee754_sqrtl (0.5L * d + 0.5L * __real__ x); s = (0.5L * __imag__ x) / r; }
__csqrtl (__complex__ long double x){  __complex__ long double res;  int rcls = fpclassify (__real__ x);  int icls = fpclassify (__imag__ x);  if (rcls <= FP_INFINITE || icls <= FP_INFINITE)    {      if (icls == FP_INFINITE)	{	  __real__ res = HUGE_VALL;	  __imag__ res = __imag__ x;	}      else if (rcls == FP_INFINITE)	{	  if (__real__ x < 0.0)	    {	      __real__ res = icls == FP_NAN ? __nanl ("") : 0;	      __imag__ res = __copysignl (HUGE_VALL, __imag__ x);	    }	  else	    {	      __real__ res = __real__ x;	      __imag__ res = (icls == FP_NAN			      ? __nanl ("") : __copysignl (0.0, __imag__ x));	    }	}      else	{	  __real__ res = __nanl ("");	  __imag__ res = __nanl ("");	}    }  else    {      if (icls == FP_ZERO)	{	  if (__real__ x < 0.0)	    {	      __real__ res = 0.0;	      __imag__ res = __copysignl (__ieee754_sqrtl (-__real__ x),					  __imag__ x);	    }	  else	    {	      __real__ res = fabsl (__ieee754_sqrtl (__real__ x));	      __imag__ res = __copysignl (0.0, __imag__ x);	    }	}      else if (rcls == FP_ZERO)	{	  long double r = __ieee754_sqrtl (0.5 * fabsl (__imag__ x));	  __real__ res = __copysignl (r, __imag__ x);	  __imag__ res = r;	}      else	{#if 0 /* FIXME: this is broken. */	  __complex__ long double q;	  long double t, r;	  if (fabsl (__imag__ x) < 2.0e-4 * fabsl (__real__ x))	    t = 0.25 * __imag__ x * (__imag__ x / __real__ x);	  else	    t = 0.5 * (__ieee754_hypotl (__real__ x, __imag__ x) - __real__ x);	  r = __ieee754_sqrtl (t);	  __real__ q = __imag__ x / (2.0 * r);	  __imag__ q = r;	  /* Heron iteration in complex arithmetic.  */	  res = 0.5 * (q + q / 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	}    }  return res;}