8 extern void sig_die(char*,int);
9 void c_div(complex *resx
, complex *a
, complex *b
)
16 if( (abr
= b
->r
) < 0.)
18 if( (abi
= b
->i
) < 0.)
23 sig_die("complex division by zero", 1);
24 ratio
= (double)b
->r
/ b
->i
;
25 den
= b
->i
* (1 + ratio
*ratio
);
26 res
.r
= (a
->r
*ratio
+ a
->i
) / den
;
27 res
.i
= (a
->i
*ratio
- a
->r
) / den
;
32 ratio
= (double)b
->i
/ b
->r
;
33 den
= b
->r
* (1 + ratio
*ratio
);
34 res
.r
= (a
->r
+ a
->i
*ratio
) / den
;
35 res
.i
= (a
->i
- a
->r
*ratio
) / den
;