VOID
pow_zi(p
, a
, b
) /* p = a**b */
doublecomplex
*p
, *a
; integer
*b
;
extern void z_div(doublecomplex
*, doublecomplex
*, doublecomplex
*);
void pow_zi(doublecomplex
*p
, doublecomplex
*a
, integer
*b
) /* p = a**b */
static doublecomplex one
= {1.0, 0.0};
t
= p
->r
* x
.r
- p
->i
* x
.i
;
p
->i
= p
->r
* x
.i
+ p
->i
* x
.r
;
t
= x
.r
* x
.r
- x
.i
* x
.i
;