make_general_replacement.h

// -*- mode: c++; c-indent-level: 4; c++-member-init-indent: 8; comment-column: 35; -*-

//-----------------------------------------------------------------------------
// make_general_replacement.h
// (c) Marc Schoenauer and Pierre Collet, 2002
/*
    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Lesser General Public
    License as published by the Free Software Foundation; either
    version 2 of the License, or (at your option) any later version.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public
    License along with this library; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

    Contact: todos@geneura.ugr.es, http://geneura.ugr.es
             Marc.Schoenauer@polytechnique.fr
             mkeijzer@dhi.dk
 */
//-----------------------------------------------------------------------------

#ifndef _make_general_replacement_h
#define _make_general_replacement_h

#include <utils/eoData.h>     // for eo_is_a_rate

// Replacement
#include <eoReduceMergeReduce.h>

// also need the parser and param includes
#include <utils/eoParser.h>
#include <utils/eoState.h>


template <class EOT>
eoReduce<EOT> & decode_reduce(eoParamParamType & _ppReduce, eoState & _state)
{
  unsigned int detSize;
  eoReduce<EOT> * ptReduce;

  // ---------- Deterministic
  if ( (_ppReduce.first == std::string("Deterministic")) ||
       (_ppReduce.first == std::string("Sequential"))
       )
  {
    ptReduce = new eoTruncate<EOT>;
  }
  // ---------- EP
  else if (_ppReduce.first == std::string("EP"))
  {
    if (!_ppReduce.second.size())   // no parameter added
    {
      std::cerr << "WARNING, no parameter passed to EP, using 6" << std::endl;
      detSize = 6;
      // put back 6 in parameter for consistency (and status file)
      _ppReduce.second.push_back(std::string("6"));
    }
    else          // parameter passed by user as EP(T)
      detSize = atoi(_ppReduce.second[0].c_str());
    ptReduce = new eoEPReduce<EOT>(detSize);
  }
  // ---------- DetTour
  else if (_ppReduce.first == std::string("DetTour"))
  {
    if (!_ppReduce.second.size())   // no parameter added
      {
        std::cerr << "WARNING, no parameter passed to DetTour, using 2" << std::endl;
        detSize = 2;
        // put back 2 in parameter for consistency (and status file)
        _ppReduce.second.push_back(std::string("2"));
      }
    else          // parameter passed by user as DetTour(T)
      detSize = atoi(_ppReduce.second[0].c_str());
    ptReduce = new eoDetTournamentTruncate<EOT>(detSize);
  }
  else if (_ppReduce.first == std::string("StochTour"))
    {
      double p;
      if (!_ppReduce.second.size())   // no parameter added
        {
          std::cerr << "WARNING, no parameter passed to StochTour, using 1" << std::endl;
          p = 1;
          // put back p in parameter for consistency (and status file)
          _ppReduce.second.push_back(std::string("1"));
        }
      else        // parameter passed by user as DetTour(T)
        {
          p = atof(_ppReduce.second[0].c_str());
          if ( (p<=0.5) || (p>1) )
            throw std::runtime_error("Stochastic tournament size should be in [0.5,1]");
        }

      ptReduce = new eoStochTournamentTruncate<EOT>(p);
    }
  else if ( (_ppReduce.first == std::string("Uniform")) ||
            (_ppReduce.first == std::string("Random"))
            )
    {
      ptReduce = new eoRandomReduce<EOT>;
    }
  else // no known reduction entered
    {
      throw std::runtime_error("Unknown reducer: " + _ppReduce.first);
    }
  // all done, stores and return a reference
  _state.storeFunctor(ptReduce);
  return (*ptReduce);
}

template <class EOT>
eoReplacement<EOT> & make_general_replacement(
    eoParser& _parser, eoState& _state,
    eoHowMany _elite = eoHowMany(0),
    bool _strongElitism = false,
    eoHowMany _surviveParents = eoHowMany(0.0),
    eoParamParamType & _reduceParentType = eoParamParamType("Deterministic"),
    eoHowMany _surviveOffspring = eoHowMany(1.0),
    eoParamParamType & _reduceOffspringType = eoParamParamType("Deterministic"),
    eoParamParamType & _reduceFinalType  = eoParamParamType("Deterministic")
        )
{
  // the replacement

  // Elitism
    eoHowMany elite =  _parser.createParam(_elite, "elite", "Nb of elite parents (percentage or absolute)", '\0', "Evolution Engine / Replacement").value();

    bool strongElitism = _parser.createParam(_strongElitism,"eliteType", "Strong (true) or weak (false) elitism (set elite to 0 for none)", '\0', "Evolution Engine / Replacement").value();

  // reduce the parents
    eoHowMany surviveParents =  _parser.createParam(_surviveParents, "surviveParents", "Nb of surviving parents (percentage or absolute)", '\0', "Evolution Engine / Replacement").value();

  eoParamParamType & reduceParentType = _parser.createParam(_reduceParentType, "reduceParents", "Parents reducer: Deterministic, EP(T), DetTour(T), StochTour(t), Uniform", '\0', "Evolution Engine / Replacement").value();

  eoReduce<EOT> & reduceParent = decode_reduce<EOT>(reduceParentType, _state);

  // reduce the offspring
    eoHowMany surviveOffspring =  _parser.createParam(_surviveOffspring, "surviveOffspring", "Nb of surviving offspring (percentage or absolute)", '\0', "Evolution Engine / Replacement").value();

  eoParamParamType & reduceOffspringType = _parser.createParam(_reduceOffspringType, "reduceOffspring", "Offspring reducer: Deterministic, EP(T), DetTour(T), StochTour(t), Uniform", '\0', "Evolution Engine / Replacement").value();

  eoReduce<EOT> & reduceOffspring = decode_reduce<EOT>(reduceOffspringType, _state);

  eoParamParamType & reduceFinalType = _parser.createParam(_reduceFinalType, "reduceFinal", "Final reducer: Deterministic, EP(T), DetTour(T), StochTour(t), Uniform", '\0', "Evolution Engine / Replacement").value();

  eoReduce<EOT> & reduceFinal = decode_reduce<EOT>(reduceFinalType, _state);

  // now the replacement itself
  eoReduceMergeReduce<EOT> *ptReplace = new eoReduceMergeReduce<EOT>(elite, strongElitism, surviveParents, reduceParent, surviveOffspring, reduceOffspring, reduceFinal);
  _state.storeFunctor(ptReplace);

  // that's it!
  return *ptReplace;
}

#endif