* Copyright (c) 1980 Regents of the University of California.
* All rights reserved. The Berkeley software License Agreement
* specifies the terms and conditions for redistribution.
double log(), exp(), cos(), sin(), atan2(), cabs();
logr
= log( cabs(a
->dreal
, a
->dimag
) );
logi
= atan2(a
->dimag
, a
->dreal
);
x
= exp( logr
* b
->dreal
- logi
* b
->dimag
);
y
= logr
* b
->dimag
+ logi
* b
->dreal
;