home *** CD-ROM | disk | FTP | other *** search
/ Linux Cubed Series 3: Developer Tools / Linux Cubed Series 3 - Developer Tools.iso / devel / lang / fortran / libf2c-0.000 / libf2c-0 / libf2c / libF77 / pow_hh.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-03-20  |  401 b   |  32 lines
  1. #include "f2c.h"
  2.  
  3. #ifdef KR_headers
  4. shortint pow_hh(ap, bp) shortint *ap, *bp;
  5. #else
  6. shortint pow_hh(shortint *ap, shortint *bp)
  7. #endif
  8. {
  9.     shortint pow, x, n;
  10.  
  11.     x = *ap;
  12.     n = *bp;
  13.  
  14.     if (n <= 0) {
  15.         if (n == 0 || x == 1)
  16.             return 1;
  17.         if (x != -1)
  18.             return x == 0 ? 1/x : 0;
  19.         n = -n;
  20.         }
  21.     for(pow = 1; ; )
  22.         {
  23.         if(n & 01)
  24.             pow *= x;
  25.         if(n >>= 1)
  26.             x *= x;
  27.         else
  28.             break;
  29.         }
  30.     return(pow);
  31.     }
  32.