ipow.h

Go to the documentation of this file.

/*
  Copyright (c) Mircea Neacsu (2014-2025) Licensed under MIT License.
  This file is part of MLIB project. See LICENSE file for full license terms.
*/

 
#pragma once
 
#if __has_include("defs.h")
#include "defs.h"
#endif
#include <assert.h>
 
namespace mlib {

template <typename T>
T ipow (T base, int exp)
{
  assert (exp >= 0);
  T result = 1;
  while (exp)
  {
    if (exp & 1)
      result *= base;
    exp >>= 1;
    base *= base;
  }
  return result;
}

template <>
inline double ipow (double base, int exp)
{
  if (exp < 0)
  {
    base = 1. / base;
    exp = -exp;
  }
  double result = 1;
  while (exp)
  {
    if (exp & 1)
      result *= base;
    exp >>= 1;
    base *= base;
  }
  return result;
}

template <typename T>
inline T squared (T base)
{
  return base * base;
}

template <typename T>
inline T cubed (T base)
{
  return base * base * base;
}
 
 
} // namespace mlib