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)