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