misc_scalar_functions.h

Go to the documentation of this file.

#ifndef TIL_MISC_SCALAR_FUNCTIONS_H_
#define TIL_MISC_SCALAR_FUNCTIONS_H_


namespace til
{
  template < typename T > INLINE const T square(const T &a) { return a*a;}  
  template < typename T > INLINE const T cube(const T &a) { return a*a*a;}
  template < typename T > INLINE const T max(const T &a, const T &b) { return ((a>b)?a:b);}
  template < typename T > INLINE const T min(const T &a, const T &b) { return ((a<b)?a:b);}
  template < typename T > INLINE const T max(const T &a, const T &b, const T &c) { return max(max(a,b),c);}
  template < typename T > INLINE const T min(const T &a, const T &b, const T &c) { return min(min(a,b),c);}
  template < typename T > INLINE const T min(const T &a, const T &b, const T &c, const T &d) { return min(min(min(a,b),c),d);}
  template < typename T > 
  inline const T crop(const T & x, const T & a, const T & b) { return min(max(x, a), b); }
  template < typename T, typename TInt >
  inline const T rangecrop(const T & x, TInt dim) { return crop<T>(x, 0, dim-1); }
  // NB: I commented it out because std::abs should be used instead, it may be faster because working at the bit level.
  //template < typename T > INLINE const T abs(const T &a) { return max(a, -a); }
  template < typename T > INLINE double norm(const T &a, const T &b) { return std::sqrt(double(a*a+b*b)); }
  template < typename T > INLINE double norm(const T &a, const T &b, const T &c) { return std::sqrt(double(a*a+b*b+c*c)); }


  template < typename T >
  inline
  T integer_pow(T x, unsigned int n)
  {
    T res = (n%2 ? x : 1);
    while (n >>= 1)
    {
      x *= x;
      if (n%2)
      {
        res *= x;
      }
    }
    return res;
  }

} // namespace til

#endif /*MISC_SCALAR_FUNCTIONS_H_*/