Gitea Mirror Bot
28 lines
430 B
C
28 lines
430 B
C
#include "f2c.h"
|
|
|
|
double pow_di(double *ap, int *bp)
|
|
{
|
|
double p = 1;
|
|
double x = *ap;
|
|
int n = *bp;
|
|
|
|
if(n != 0)
|
|
{
|
|
if(n < 0)
|
|
{
|
|
n = -n;
|
|
x = 1/x;
|
|
}
|
|
unsigned u = (unsigned)n;
|
|
for(;;)
|
|
{
|
|
if((u & 1) != 0)
|
|
p *= x;
|
|
if((u >>= 1) == 0)
|
|
break;
|
|
x *= x;
|
|
}
|
|
}
|
|
return p;
|
|
}
|