summaryrefslogtreecommitdiffstats
path: root/lib/libF77/pow_qq.c
diff options
context:
space:
mode:
authorljo <ljo@FreeBSD.org>1994-10-26 18:17:41 +0000
committerljo <ljo@FreeBSD.org>1994-10-26 18:17:41 +0000
commit2195d627fb793fec6c00bdddf1c4d46982b9f27f (patch)
treec2961f38c1507b80baefaca1d568747c113f28c0 /lib/libF77/pow_qq.c
parent80b2e1963004e4ae0173b4f9bc965518ba81d3b0 (diff)
downloadFreeBSD-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_qq.c')
-rw-r--r--lib/libF77/pow_qq.c31
1 files changed, 31 insertions, 0 deletions
diff --git a/lib/libF77/pow_qq.c b/lib/libF77/pow_qq.c
new file mode 100644
index 0000000..d80c40a9
--- /dev/null
+++ b/lib/libF77/pow_qq.c
@@ -0,0 +1,31 @@
+#include "f2c.h"
+
+#ifdef KR_headers
+longint pow_qq(ap, bp) longint *ap, *bp;
+#else
+longint pow_qq(longint *ap, longint *bp)
+#endif
+{
+ longint pow, x, n;
+
+ x = *ap;
+ n = *bp;
+
+ if (n <= 0) {
+ if (n == 0 || x == 1)
+ return 1;
+ if (x != -1)
+ return x == 0 ? 1/x : 0;
+ n = -n;
+ }
+ for(pow = 1; ; )
+ {
+ if(n & 01)
+ pow *= x;
+ if(n >>= 1)
+ x *= x;
+ else
+ break;
+ }
+ return(pow);
+ }
OpenPOWER on IntegriCloud