VOID
z_div(c
, a
, b
) doublecomplex
*a
, *b
, *c
;
extern void sig_die(char*, int);
void z_div(doublecomplex
*c
, doublecomplex
*a
, doublecomplex
*b
)
sig_die("complex division by zero", 1);
den
= b
->i
* (1 + ratio
*ratio
);
c
->r
= (a
->r
*ratio
+ a
->i
) / den
;
c
->i
= (a
->i
*ratio
- a
->r
) / den
;
den
= b
->r
* (1 + ratio
*ratio
);
c
->r
= (a
->r
+ a
->i
*ratio
) / den
;
c
->i
= (a
->i
- a
->r
*ratio
) / den
;