freebsd-dev/contrib/libf2c/libF77/pow_di.c
2003-07-11 03:42:19 +00:00

33 lines
347 B
C

#include "f2c.h"
double
pow_di (doublereal * ap, integer * bp)
{
double pow, x;
integer n;
unsigned long u;
pow = 1;
x = *ap;
n = *bp;
if (n != 0)
{
if (n < 0)
{
n = -n;
x = 1 / x;
}
for (u = n;;)
{
if (u & 01)
pow *= x;
if (u >>= 1)
x *= x;
else
break;
}
}
return (pow);
}