4 double log(), exp(), cos(), sin(), atan2(), f__cabs();
5 VOID
pow_zz(r
,a
,b
) doublecomplex
*r
, *a
, *b
;
9 extern double f__cabs(double,double);
10 void pow_zz(doublecomplex
*r
, doublecomplex
*a
, doublecomplex
*b
)
13 double logr
, logi
, x
, y
;
15 logr
= log( f__cabs(a
->r
, a
->i
) );
16 logi
= atan2(a
->i
, a
->r
);
18 x
= exp( logr
* b
->r
- logi
* b
->i
);
19 y
= logr
* b
->i
+ logi
* b
->r
;