summaryrefslogtreecommitdiffstats
path: root/contrib/libf2c/libF77/pow_di.c
diff options
context:
space:
mode:
Diffstat (limited to 'contrib/libf2c/libF77/pow_di.c')
-rw-r--r--contrib/libf2c/libF77/pow_di.c51
1 files changed, 24 insertions, 27 deletions
diff --git a/contrib/libf2c/libF77/pow_di.c b/contrib/libf2c/libF77/pow_di.c
index affed62..d2298a0 100644
--- a/contrib/libf2c/libF77/pow_di.c
+++ b/contrib/libf2c/libF77/pow_di.c
@@ -1,35 +1,32 @@
#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_di (doublereal * ap, integer * bp)
{
-double pow, x;
-integer n;
-unsigned long u;
+ double pow, x;
+ integer n;
+ unsigned long u;
-pow = 1;
-x = *ap;
-n = *bp;
+ pow = 1;
+ x = *ap;
+ n = *bp;
-if(n != 0)
+ if (n != 0)
+ {
+ 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;
- }
+ n = -n;
+ x = 1 / x;
}
-return(pow);
+ for (u = n;;)
+ {
+ if (u & 01)
+ pow *= x;
+ if (u >>= 1)
+ x *= x;
+ else
+ break;
+ }
+ }
+ return (pow);
}
OpenPOWER on IntegriCloud