5 VOID
z_div(resx
, a
, b
) doublecomplex
*a
, *b
, *resx
;
7 extern void sig_die(char*, int);
8 void z_div(doublecomplex
*resx
, doublecomplex
*a
, doublecomplex
*b
)
15 if( (abr
= b
->r
) < 0.)
17 if( (abi
= b
->i
) < 0.)
22 sig_die("complex division by zero", 1);
24 den
= b
->i
* (1 + ratio
*ratio
);
25 res
.r
= (a
->r
*ratio
+ a
->i
) / den
;
26 res
.i
= (a
->i
*ratio
- a
->r
) / den
;
32 den
= b
->r
* (1 + ratio
*ratio
);
33 res
.r
= (a
->r
+ a
->i
*ratio
) / den
;
34 res
.i
= (a
->i
- a
->r
*ratio
) / den
;