freebsd-dev/contrib/libf2c/libF77/pow_hh.c

#include "f2c.h"

shortint
pow_hh (shortint * ap, shortint * bp)
{
  shortint pow, x, n;
  unsigned u;

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