aboutsummaryrefslogtreecommitdiff
path: root/gcc/f/runtime/libF77/pow_hh.c
diff options
context:
space:
mode:
Diffstat (limited to 'gcc/f/runtime/libF77/pow_hh.c')
-rw-r--r--gcc/f/runtime/libF77/pow_hh.c33
1 files changed, 33 insertions, 0 deletions
diff --git a/gcc/f/runtime/libF77/pow_hh.c b/gcc/f/runtime/libF77/pow_hh.c
new file mode 100644
index 00000000000..24a019734da
--- /dev/null
+++ b/gcc/f/runtime/libF77/pow_hh.c
@@ -0,0 +1,33 @@
+#include "f2c.h"
+
+#ifdef KR_headers
+shortint pow_hh(ap, bp) shortint *ap, *bp;
+#else
+shortint pow_hh(shortint *ap, shortint *bp)
+#endif
+{
+ shortint pow, x, n;
+ unsigned 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);
+ }