mirror of
git://gcc.gnu.org/git/gcc.git
synced 2024-12-26 07:08:57 +08:00
32 lines
299 B
C
32 lines
299 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);
|
||
|
|
}
|