4 VOID pow_zi(p, a, b) /* p = a**b */
5 doublecomplex *p, *a; integer *b;
7 extern void z_div(doublecomplex*, doublecomplex*, doublecomplex*);
8 void pow_zi(doublecomplex *p, doublecomplex *a, integer *b) /* p = a**b */
15 static doublecomplex one = {1.0, 0.0};
38 t = q.r * x.r - q.i * x.i;
39 q.i = q.r * x.i + q.i * x.r;
44 t = x.r * x.r - x.i * x.i;