#ifndef ICL_CLIPPED_CAST_H
#define ICL_CLIPPED_CAST_H

#include <limits>

namespace icl{
  /// clips a value into the range [tMin,tMax] \ingroup GENERAL
  template <class T>
  inline T clip(T tX, T tMin, T tMax){ 
    return tX < tMin ? tMin : tX > tMax ? tMax : tX; 
  }
  
  /// utility cast function wrapping the standard lib's numerical_limits template
  template<class S, class D> 
  inline D clipped_cast(S src){
    return src < std::numeric_limits<D>::min() ? std::numeric_limits<D>::min() : 
           src > std::numeric_limits<D>::max() ? std::numeric_limits<D>::max() : 
           static_cast<D>(src);
  }
  
  /** \cond */
  /// specializations for all buildin data types
#define SPECIALISE_CLIPPED_CAST(T) template<> inline T clipped_cast<T,T>(T t) { return t; }
  SPECIALISE_CLIPPED_CAST(int)
  SPECIALISE_CLIPPED_CAST(unsigned int)
  SPECIALISE_CLIPPED_CAST(char)
  SPECIALISE_CLIPPED_CAST(unsigned char)
  SPECIALISE_CLIPPED_CAST(short)
  SPECIALISE_CLIPPED_CAST(unsigned short)
  SPECIALISE_CLIPPED_CAST(long)
  SPECIALISE_CLIPPED_CAST(unsigned long)
  SPECIALISE_CLIPPED_CAST(bool)
  SPECIALISE_CLIPPED_CAST(float)
  SPECIALISE_CLIPPED_CAST(double)
#undef SPECIALISE_CLIPPED_CAST
  /** \endcond */

}


#endif