diff options
Diffstat (limited to 'gnu/egcs/libf2c/libF77/pow_ii.c')
-rw-r--r-- | gnu/egcs/libf2c/libF77/pow_ii.c | 33 |
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); + } |