find_bin.h

Go to the documentation of this file.

 
 
#ifndef IDTPM2_BINS_FIND_BIN
#define IDTPM2_BINS_FIND_BIN
 
#include <vector>
#include <stdexcept>
 
 
namespace IDTPM2::bins {
 
template<typename T,typename U=T> 
size_t find_bin( const std::vector<T>& v, const U& x ) {
  if ( v.size()<2) throw std::runtime_error("cannot define a bin range with fewer than two limits");              
  if ( x>=v.back() ) return v.size()-2;
  if ( x<v.front() ) return 0;
  typename std::vector<T>::const_iterator itr = std::lower_bound( v.begin(), v.end(), x );
  return (itr-v.begin())-1;
}
 
}
 
#endif // IDTPM2_BINS_FIND_BIN