freebsd-dev/contrib/libf2c/libF77/pow_hh.c

32 lines
386 B
C
Raw Normal View History

1999-09-18 10:51:31 +00:00
#include "f2c.h"
2003-07-11 03:42:19 +00:00
shortint
pow_hh (shortint * ap, shortint * bp)
1999-09-18 10:51:31 +00:00
{
2003-07-11 03:42:19 +00:00
shortint pow, x, n;
unsigned u;
1999-09-18 10:51:31 +00:00
2003-07-11 03:42:19 +00:00
x = *ap;
n = *bp;
1999-09-18 10:51:31 +00:00
2003-07-11 03:42:19 +00:00
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);
}