summaryrefslogtreecommitdiff
path: root/libf2c/libF77/pow_di.c
diff options
context:
space:
mode:
authorlaw <law@138bc75d-0d04-0410-961f-82ee72b054a4>1998-02-01 01:37:08 +0000
committerlaw <law@138bc75d-0d04-0410-961f-82ee72b054a4>1998-02-01 01:37:08 +0000
commita9fe81af3b8d02e1529b3f1f2a679b93cd9c36e1 (patch)
tree4b30cdf598ee6a8b24903fba52bb3e0c06bf459b /libf2c/libF77/pow_di.c
parent9682c2b47acc059592d168a09c7162f2bf333ee2 (diff)
downloadgcc-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.c35
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);
+}