eoDetSelect.h

//-----------------------------------------------------------------------------

#ifndef _eoDetSelect_h
#define _eoDetSelect_h


//-----------------------------------------------------------------------------
#include <eoSelect.h>
#include <utils/eoHowMany.h>
#include <math.h>
//-----------------------------------------------------------------------------

template<class EOT>
class eoDetSelect : public eoSelect<EOT>
{
 public:
     eoDetSelect(double _rate = 1.0, bool _interpret_as_rate = true)
         : howMany(_rate, _interpret_as_rate)  {}

  virtual void operator()(const eoPop<EOT>& _source, eoPop<EOT>& _dest)
  {
    unsigned int pSize = _source.size();
    size_t target = howMany(pSize);

    if ( target == 0 )
        {
            eo::log << eo::warnings << "Call to a eoHowMany instance returns 0 (target=" << target << ") it will be replaced by 1 to continue." << std::endl;
            target = 1;
        }

    _dest.resize(target);

    unsigned remain = target % pSize;
    unsigned entireCopy = target / pSize;
    typename eoPop<EOT>::iterator it = _dest.begin();

    if (target >= pSize)
      {
        for (unsigned i=0; i<entireCopy; i++)
          {
              std::copy(_source.begin(), _source.end(), it);
            it += pSize;
          }
      }
    // the last ones
    if (remain)
      {
          std::copy(_source.begin(), _source.begin()+remain, it);
      }
  }

private :
  eoHowMany howMany;
};

#endif