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_di.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_di.c')
-rw-r--r-- | lib/libF77/pow_di.c | 34 |
1 files changed, 34 insertions, 0 deletions
diff --git a/lib/libF77/pow_di.c b/lib/libF77/pow_di.c new file mode 100644 index 0000000..7af69a7 --- /dev/null +++ b/lib/libF77/pow_di.c @@ -0,0 +1,34 @@ +#include "f2c.h" + +#ifdef KR_headers +double pow_di(ap, bp) doublereal *ap; integer *bp; +#else +double pow_di(doublereal *ap, integer *bp) +#endif +{ +double pow, x; +integer n; + +pow = 1; +x = *ap; +n = *bp; + +if(n != 0) + { + if(n < 0) + { + n = -n; + x = 1/x; + } + for( ; ; ) + { + if(n & 01) + pow *= x; + if(n >>= 1) + x *= x; + else + break; + } + } +return(pow); +} |