1 files changed, 31 insertions, 0 deletions
diff --git a/contrib/libf2c/libF77/pow_hh.c b/contrib/libf2c/libF77/pow_hh.c
new file mode 100644
index 0000000..3379d8a
--- /dev/null
+++ b/contrib/libf2c/libF77/pow_hh.c
@@ -0,0 +1,31 @@
+#include "f2c.h"
+
+shortint
+pow_hh (shortint * ap, shortint * bp)
+{
+  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);
+}