summaryrefslogtreecommitdiff
path: root/gnu/egcs/libf2c/libF77/pow_ii.c
diff options
context:
space:
mode:
Diffstat (limited to 'gnu/egcs/libf2c/libF77/pow_ii.c')
-rw-r--r--gnu/egcs/libf2c/libF77/pow_ii.c33
1 files changed, 33 insertions, 0 deletions
diff --git a/gnu/egcs/libf2c/libF77/pow_ii.c b/gnu/egcs/libf2c/libF77/pow_ii.c
new file mode 100644
index 00000000000..84d1c7e0b5e
--- /dev/null
+++ b/gnu/egcs/libf2c/libF77/pow_ii.c
@@ -0,0 +1,33 @@
+#include "f2c.h"
+
+#ifdef KR_headers
+integer pow_ii(ap, bp) integer *ap, *bp;
+#else
+integer pow_ii(integer *ap, integer *bp)
+#endif
+{
+ integer pow, x, n;
+ unsigned long u;
+
+ 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;
+ }
+ u = n;
+ for(pow = 1; ; )
+ {
+ if(u & 01)
+ pow *= x;
+ if(u >>= 1)
+ x *= x;
+ else
+ break;
+ }
+ return(pow);
+ }