4 VOID
pow_zi(resx
, a
, b
) /* p = a**b */
5 doublecomplex
*resx
, *a
; integer
*b
;
7 extern void z_div(doublecomplex
*, doublecomplex
*, doublecomplex
*);
8 void pow_zi(doublecomplex
*resx
, doublecomplex
*a
, integer
*b
) /* p = a**b */
16 static doublecomplex one
= {1.0, 0.0};
45 t
= res
.r
* x
.r
- res
.i
* x
.i
;
46 res
.i
= res
.r
* x
.i
+ res
.i
* x
.r
;
51 t
= x
.r
* x
.r
- x
.i
* x
.i
;