f2c libF77 from netlib.att.com Jan 4 1994

[unix-history] / lib / libF77 / pow_hh.c

#include "f2c.h"

#ifdef KR_headers
shortint pow_hh(ap, bp) shortint *ap, *bp;
#else
shortint pow_hh(shortint *ap, shortint *bp)
#endif
{
       shortint pow, x, n;

       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;
               }
       for(pow = 1; ; )
               {
               if(n & 01)
                       pow *= x;
               if(n >>= 1)
                       x *= x;
               else
                       break;
               }
       return(pow);
       }