diff options
author | law <law@138bc75d-0d04-0410-961f-82ee72b054a4> | 1998-02-01 01:37:08 +0000 |
---|---|---|
committer | law <law@138bc75d-0d04-0410-961f-82ee72b054a4> | 1998-02-01 01:37:08 +0000 |
commit | a9fe81af3b8d02e1529b3f1f2a679b93cd9c36e1 (patch) | |
tree | 4b30cdf598ee6a8b24903fba52bb3e0c06bf459b /libf2c/libF77/pow_di.c | |
parent | 9682c2b47acc059592d168a09c7162f2bf333ee2 (diff) | |
download | gcc-a9fe81af3b8d02e1529b3f1f2a679b93cd9c36e1.tar.gz |
* Previous contents of gcc/f/runtime moved into toplevel
"libf2c" directory.
git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/trunk@17568 138bc75d-0d04-0410-961f-82ee72b054a4
Diffstat (limited to 'libf2c/libF77/pow_di.c')
-rw-r--r-- | libf2c/libF77/pow_di.c | 35 |
1 files changed, 35 insertions, 0 deletions
diff --git a/libf2c/libF77/pow_di.c b/libf2c/libF77/pow_di.c new file mode 100644 index 00000000000..affed625a91 --- /dev/null +++ b/libf2c/libF77/pow_di.c @@ -0,0 +1,35 @@ +#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; +unsigned long u; + +pow = 1; +x = *ap; +n = *bp; + +if(n != 0) + { + if(n < 0) + { + n = -n; + x = 1/x; + } + for(u = n; ; ) + { + if(u & 01) + pow *= x; + if(u >>= 1) + x *= x; + else + break; + } + } +return(pow); +} |