2
0
mirror of git://gcc.gnu.org/git/gcc.git synced 2024-12-22 00:59:55 +08:00
gcc/libf2c/libF77/pow_di.c
Jeff Law 81fea2b1d1 * Previous contents of gcc/f/runtime moved into toplevel
"libf2c" directory.

From-SVN: r17568
1998-01-31 18:37:08 -07:00

36 lines
381 B
C

#include "f2c.h"
#ifdef KR_headers
double pow_di(ap, bp) doublereal *ap; integer *bp;
#else
double pow_di(doublereal *ap, integer *bp)
#endif
{
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);
}