mirror of
https://github.com/Cantera/cantera.git
synced 2025-02-25 18:55:29 -06:00
40 lines
489 B
C
40 lines
489 B
C
#include "f2c.h"
|
|
#ifdef __cplusplus
|
|
extern "C" {
|
|
#endif
|
|
|
|
#ifdef KR_headers
|
|
shortint pow_hh(ap, bp) shortint *ap, *bp;
|
|
#else
|
|
shortint pow_hh(shortint *ap, shortint *bp)
|
|
#endif
|
|
{
|
|
shortint pow, x, n;
|
|
unsigned u;
|
|
|
|
x = *ap;
|
|
n = *bp;
|
|
|
|
if (n <= 0) {
|
|
if (n == 0 || x == 1)
|
|
return 1;
|
|
if (x != -1)
|
|
return x == 0 ? 1/x : 0;
|
|
n = -n;
|
|
}
|
|
u = n;
|
|
for(pow = 1; ; )
|
|
{
|
|
if(u & 01)
|
|
pow *= x;
|
|
if(u >>= 1)
|
|
x *= x;
|
|
else
|
|
break;
|
|
}
|
|
return(pow);
|
|
}
|
|
#ifdef __cplusplus
|
|
}
|
|
#endif
|