diff options
Diffstat (limited to 'contrib/libf2c/libF77/z_div.c')
| -rw-r--r-- | contrib/libf2c/libF77/z_div.c | 53 |
1 files changed, 28 insertions, 25 deletions
diff --git a/contrib/libf2c/libF77/z_div.c b/contrib/libf2c/libF77/z_div.c index 22153fa4514a..4a987ab255a4 100644 --- a/contrib/libf2c/libF77/z_div.c +++ b/contrib/libf2c/libF77/z_div.c @@ -2,35 +2,38 @@ #ifdef KR_headers extern VOID sig_die(); -VOID z_div(c, a, b) doublecomplex *a, *b, *c; +VOID z_div(resx, a, b) doublecomplex *a, *b, *resx; #else extern void sig_die(char*, int); -void z_div(doublecomplex *c, doublecomplex *a, doublecomplex *b) +void z_div(doublecomplex *resx, doublecomplex *a, doublecomplex *b) #endif { - double ratio, den; - double abr, abi, cr; +double ratio, den; +double abr, abi; +doublecomplex res; - 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); - cr = (a->r*ratio + a->i) / den; - c->i = (a->i*ratio - a->r) / den; - } +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); + res.r = (a->r*ratio + a->i) / den; + res.i = (a->i*ratio - a->r) / den; + } - else - { - ratio = b->i / b->r ; - den = b->r * (1 + ratio*ratio); - cr = (a->r + a->i*ratio) / den; - c->i = (a->i - a->r*ratio) / den; - } - c->r = cr; +else + { + ratio = b->i / b->r ; + den = b->r * (1 + ratio*ratio); + res.r = (a->r + a->i*ratio) / den; + res.i = (a->i - a->r*ratio) / den; } + +resx->r = res.r; +resx->i = res.i; +} |