diff options
Diffstat (limited to 'Build/source/utils/pmx/pmx-2.6.18/libf2c/pow_di.c')
-rw-r--r-- | Build/source/utils/pmx/pmx-2.6.18/libf2c/pow_di.c | 41 |
1 files changed, 41 insertions, 0 deletions
diff --git a/Build/source/utils/pmx/pmx-2.6.18/libf2c/pow_di.c b/Build/source/utils/pmx/pmx-2.6.18/libf2c/pow_di.c new file mode 100644 index 00000000000..abf36cb74c0 --- /dev/null +++ b/Build/source/utils/pmx/pmx-2.6.18/libf2c/pow_di.c @@ -0,0 +1,41 @@ +#include "f2c.h" +#ifdef __cplusplus +extern "C" { +#endif + +#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); +} +#ifdef __cplusplus +} +#endif |