4 extern double sqrt(), f__cabs();
6 VOID
c_sqrt(resx
, z
) complex *resx
, *z
;
10 extern double f__cabs(double, double);
12 void c_sqrt(complex *resx
, complex *z
)
18 if( (mag
= f__cabs(z
->r
, z
->i
)) == 0.)
22 res
.r
= t
= sqrt(0.5 * (mag
+ z
->r
) );
28 t
= sqrt(0.5 * (mag
- z
->r
) );