BasicMathFunctions.h

#pragma once

#if defined(BasicMathFunctions_RECURSES)
#error Recursive header files inclusion detected in BasicMathFunctions.h
#else // defined(BasicMathFunctions_RECURSES)

#define BasicMathFunctions_RECURSES

#if !defined BasicMathFunctions_h

#define BasicMathFunctions_h

// Inclusions
#include <algorithm>
#include <functional>
#include "DGtal/base/Common.h"
#include "DGtal/base/Bits.h"
#include "DGtal/kernel/NumberTraits.h"

namespace DGtal 
{
  namespace functions
  {

    template<typename T>
    T power(const T&aVal, const unsigned int exponent)
    {
      unsigned int q=exponent;
      T p(aVal);
      
      if (exponent == 0) { return 1;    }
      
      T result = NumberTraits<T>::ONE;
      while (q != 0) 
        {
          if (q % 2 == 1) {    // q is odd
            result *= p;
            q--;
          }
          p *= p;
          q /= 2;
        }
      return result;
    }


    template<typename T> 
    T roundToUpperPowerOfTwo(const  T &n){
      return (T)  1 << (1+DGtal::Bits::mostSignificantBit( (unsigned int) n-1 ) );
    }

    template<typename T>
    T abs(const T & a)
    {
      BOOST_CONCEPT_ASSERT((boost::LessThanComparable<T>));
      if (a<0) 
        return -a;
      else
        return a;
    }

    template <typename T>
    inline
    T square( T x )
    { return x * x; }

    template <typename T>
    inline
    T cube( T x )
    { return x * x * x; }

  }  // namespace functions
} // namespace DGTal




#endif // !defined BasicMathFunctions_h

#undef BasicMathFunctions_RECURSES
#endif // else defined(BasicMathFunctions_RECURSES)