freebsd-dev/lib/libF77/pow_ri.c
L Jonas Olsson 09c656ca68 Library for f2c (part 2 of 2)
Obtained from: netlib.att.com
1994-10-26 18:17:41 +00:00

35 lines
348 B
C

#include "f2c.h"
#ifdef KR_headers
double pow_ri(ap, bp) real *ap; integer *bp;
#else
double pow_ri(real *ap, integer *bp)
#endif
{
double pow, x;
integer n;
pow = 1;
x = *ap;
n = *bp;
if(n != 0)
{
if(n < 0)
{
n = -n;
x = 1/x;
}
for( ; ; )
{
if(n & 01)
pow *= x;
if(n >>= 1)
x *= x;
else
break;
}
}
return(pow);
}