eoG3Replacement.h

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

#ifndef _eoG3Replacement_h
#define _eoG3Replacement_h


//-----------------------------------------------------------------------------
#include <eoPop.h>
#include <eoFunctor.h>
#include <eoMerge.h>
#include <eoReduce.h>
#include <utils/eoHowMany.h>
#include <eoReduceSplit.h>
//-----------------------------------------------------------------------------

template <class EOT>
class eoG3Replacement : public eoReplacement<EOT>
{
public:
  eoG3Replacement(eoHowMany _howManyEliminatedParents = eoHowMany(2, false)) :
    // split truncates the parents and returns eliminated parents
    split(_howManyEliminatedParents, true),
    // reduce truncates the offpsring and does not return eliminated guys
    reduce(-_howManyEliminatedParents, false)
  {}

    void operator()(eoPop<EOT> & _parents, eoPop<EOT> & _offspring)
    {
      eoPop<EOT> temp;
      split(_parents, temp);
      unsigned toKeep = temp.size(); // how many to keep from merged populations
      // merge temp into offspring
      plus(temp, _offspring);      // add temp to _offspring (a little inconsistent!)

      // reduce merged
      reduce(_offspring, temp);    // temp dummy arg. will not be modified
      // minimla check:
      if (_offspring.size() != toKeep)
        {
          std::cerr << "Les tailles " << _offspring.size() << " " << toKeep << std::endl;
        throw std::runtime_error("eoG3Replacement: wrong number of remaining offspring");
        }
      // and put back into _parents
      plus(_offspring, _parents);
    }

private:
  eoLinearTruncateSplit<EOT> split; // few parents to truncate -> linear
  eoTruncateSplit<EOT> reduce;     // supposedly many offspring to truncate
  eoPlus<EOT> plus;
};

#endif