8 extern void sig_die(char*,int);
9 void c_div(complex *c
, complex *a
, complex *b
)
15 if( (abr
= b
->r
) < 0.)
17 if( (abi
= b
->i
) < 0.)
22 sig_die("complex division by zero", 1);
23 ratio
= (double)b
->r
/ b
->i
;
24 den
= b
->i
* (1 + ratio
*ratio
);
25 cr
= (a
->r
*ratio
+ a
->i
) / den
;
26 c
->i
= (a
->i
*ratio
- a
->r
) / den
;
31 ratio
= (double)b
->i
/ b
->r
;
32 den
= b
->r
* (1 + ratio
*ratio
);
33 cr
= (a
->r
+ a
->i
*ratio
) / den
;
34 c
->i
= (a
->i
- a
->r
*ratio
) / den
;