diff options
author | ljo <ljo@FreeBSD.org> | 1994-10-26 18:17:41 +0000 |
---|---|---|
committer | ljo <ljo@FreeBSD.org> | 1994-10-26 18:17:41 +0000 |
commit | 2195d627fb793fec6c00bdddf1c4d46982b9f27f (patch) | |
tree | c2961f38c1507b80baefaca1d568747c113f28c0 /lib/libF77/pow_zi.c | |
parent | 80b2e1963004e4ae0173b4f9bc965518ba81d3b0 (diff) | |
download | FreeBSD-src-2195d627fb793fec6c00bdddf1c4d46982b9f27f.zip FreeBSD-src-2195d627fb793fec6c00bdddf1c4d46982b9f27f.tar.gz |
Library for f2c (part 2 of 2)
Obtained from: netlib.att.com
Diffstat (limited to 'lib/libF77/pow_zi.c')
-rw-r--r-- | lib/libF77/pow_zi.c | 50 |
1 files changed, 50 insertions, 0 deletions
diff --git a/lib/libF77/pow_zi.c b/lib/libF77/pow_zi.c new file mode 100644 index 0000000..8dd6006 --- /dev/null +++ b/lib/libF77/pow_zi.c @@ -0,0 +1,50 @@ +#include "f2c.h" + +#ifdef KR_headers +VOID pow_zi(p, a, b) /* p = a**b */ + doublecomplex *p, *a; integer *b; +#else +extern void z_div(doublecomplex*, doublecomplex*, doublecomplex*); +void pow_zi(doublecomplex *p, doublecomplex *a, integer *b) /* p = a**b */ +#endif +{ +integer n; +double t; +doublecomplex x; +static doublecomplex one = {1.0, 0.0}; + +n = *b; +p->r = 1; +p->i = 0; + +if(n == 0) + return; +if(n < 0) + { + n = -n; + z_div(&x, &one, a); + } +else + { + x.r = a->r; + x.i = a->i; + } + +for( ; ; ) + { + if(n & 01) + { + t = p->r * x.r - p->i * x.i; + p->i = p->r * x.i + p->i * x.r; + p->r = t; + } + if(n >>= 1) + { + t = x.r * x.r - x.i * x.i; + x.i = 2 * x.r * x.i; + x.r = t; + } + else + break; + } +} |