BSD 4_2 release

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

/*
*      "@(#)pow_hh.c   1.1"
*/

short pow_hh(ap, bp)
short *ap, *bp;
{
short pow, x, n;

pow = 1;
x = *ap;
n = *bp;

if(n < 0)
       { }
else if(n > 0)
       for( ; ; )
               {
               if(n & 01)
                       pow *= x;
               if(n >>= 1)
                       x *= x;
               else
                       break;
               }
return(pow);
}