diff options
author | jmz <jmz@FreeBSD.org> | 1997-04-13 01:16:58 +0000 |
---|---|---|
committer | jmz <jmz@FreeBSD.org> | 1997-04-13 01:16:58 +0000 |
commit | 438bdd8a9c9d390f13a4b05980c19b18e974d16c (patch) | |
tree | e809113557ae09d5ce13e9239ae81534708bd955 /lib/libF77/z_div.c | |
parent | 5d93c9d5fb208d10eacf608b44ee02d3cd5b4a16 (diff) | |
download | FreeBSD-src-438bdd8a9c9d390f13a4b05980c19b18e974d16c.zip FreeBSD-src-438bdd8a9c9d390f13a4b05980c19b18e974d16c.tar.gz |
Upgrade to the 1997/02/26 version.
Diffstat (limited to 'lib/libF77/z_div.c')
-rw-r--r-- | lib/libF77/z_div.c | 46 |
1 files changed, 23 insertions, 23 deletions
diff --git a/lib/libF77/z_div.c b/lib/libF77/z_div.c index fd53733..ed7ee66 100644 --- a/lib/libF77/z_div.c +++ b/lib/libF77/z_div.c @@ -8,29 +8,29 @@ extern void sig_die(char*, int); void z_div(doublecomplex *c, doublecomplex *a, doublecomplex *b) #endif { -double ratio, den; -double abr, abi; + double ratio, den; + double abr, abi; + double ai = a->i, ar = a->r, bi = b->i, br = b->r; -if( (abr = b->r) < 0.) - abr = - abr; -if( (abi = b->i) < 0.) - abi = - abi; -if( abr <= abi ) - { - if(abi == 0) - sig_die("complex division by zero", 1); - ratio = b->r / b->i ; - den = b->i * (1 + ratio*ratio); - c->r = (a->r*ratio + a->i) / den; - c->i = (a->i*ratio - a->r) / den; - } + if( (abr = br) < 0.) + abr = - abr; + if( (abi = bi) < 0.) + abi = - abi; + if( abr <= abi ) + { + if(abi == 0) + sig_die("complex division by zero", 1); + ratio = br / bi ; + den = bi * (1 + ratio*ratio); + c->r = (ar*ratio + ai) / den; + c->i = (ai*ratio - ar) / den; + } -else - { - ratio = b->i / b->r ; - den = b->r * (1 + ratio*ratio); - c->r = (a->r + a->i*ratio) / den; - c->i = (a->i - a->r*ratio) / den; + else + { + ratio = bi / br ; + den = br * (1 + ratio*ratio); + c->r = (ar + ai*ratio) / den; + c->i = (ai - ar*ratio) / den; + } } - -} |