make_op_es.h

/* (c) Maarten Keijzer, Marc Schoenauer and GeNeura Team, 2001

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: http://eodev.sourceforge.net
         todos@geneura.ugr.es, http://geneura.ugr.es
         Marc.Schoenauer@polytechnique.fr
         mkeijzer@dhi.dk
 */


#ifndef EO_make_op_h
#define EO_make_op_h

// the operators
#include <eoOp.h>
#include <eoGenOp.h>
#include <eoCloneOps.h>
#include <eoOpContainer.h>
// combinations of simple eoOps (eoMonOp and eoQuadOp)
#include <eoProportionalCombinedOp.h>

// the specialized Real stuff
#include <es/eoReal.h>
#include <es/eoRealAtomXover.h>
#include <es/eoEsChromInit.h>
#include <es/eoEsMutationInit.h>
#include <es/eoEsMutate.h>
#include <es/eoEsGlobalXover.h>
#include <es/eoEsStandardXover.h>
  // also need the parser and param includes
#include <utils/eoParser.h>
#include <utils/eoState.h>


/*
 * This function builds the operators that will be applied to the eoReal
 *
 * It uses a parser (to get user parameters) and a state (to store the memory)
 * the last argument is an individual, needed for 2 reasons
 *     it disambiguates the call after instanciations
 *     some operator might need some private information about the indis
 *
 * This is why the template is the complete EOT even though only the fitness
 * is actually templatized here: the following only applies to bitstrings
 *
 * Note : the last parameter is an eoInit: if some operator needs some info
 *        about the gneotypes, the init has it all (e.g. bounds, ...)
 *        Simply do
 *        EOT myEO;
 *        _init(myEO);
 *        and myEO is then an ACTUAL object
*/

template <class EOT>
eoGenOp<EOT> & do_make_op(eoParser& _parser, eoState& _state, eoRealInitBounded<EOT>& _init)
{
  // get std::vector size
  unsigned vecSize = _init.size();

  // First, decide whether the objective variables are bounded
  eoValueParam<eoRealVectorBounds>& boundsParam
      = _parser.getORcreateParam(eoRealVectorBounds(vecSize,eoDummyRealNoBounds),
                                 "objectBounds", "Bounds for variables",
                                 'B', "Variation Operators");

  std::cerr << boundsParam.value() << std::endl;

    // now we read Pcross and Pmut,
  eoValueParam<std::string>& operatorParam
      = _parser.getORcreateParam(std::string("SGA"), "operator",
                                 "Description of the operator (SGA only now)",
                                 'o', "Variation Operators");

  if (operatorParam.value() != std::string("SGA"))
    throw std::runtime_error("Sorry, only SGA-like operator available right now\n");

    // now we read Pcross and Pmut,
    // and create the eoGenOp that is exactly
    // crossover with pcross + mutation with pmut

  eoValueParam<double>& pCrossParam
      = _parser.getORcreateParam(1.0, "pCross", "Probability of Crossover",
                                 'C', "Variation Operators" );
  // minimum check
  if ( (pCrossParam.value() < 0) || (pCrossParam.value() > 1) )
    throw std::runtime_error("Invalid pCross");

  eoValueParam<double>& pMutParam
      = _parser.getORcreateParam(1.0, "pMut", "Probability of Mutation",
                                 'M', "Variation Operators" );
  // minimum check
  if ( (pMutParam.value() < 0) || (pMutParam.value() > 1) )
    throw std::runtime_error("Invalid pMut");


  // crossover
  // ES crossover
  eoValueParam<std::string>& crossTypeParam
      = _parser.getORcreateParam(std::string("global"), "crossType",
                                 "Type of ES recombination (global or standard)",
                                 'C', "Variation Operators");

  eoValueParam<std::string>& crossObjParam
      = _parser.getORcreateParam(std::string("discrete"), "crossObj",
                                 "Recombination of object variables (discrete, intermediate or none)",
                                 'O', "Variation Operators");
  eoValueParam<std::string>& crossStdevParam
      = _parser.getORcreateParam(std::string("intermediate"), "crossStdev",
                                 "Recombination of mutation strategy parameters "
                                 "(intermediate, discrete or none)",
                                 'S', "Variation Operators");

  // The pointers: first the atom Xover
  eoBinOp<double> *ptObjAtomCross = NULL;
  eoBinOp<double> *ptStdevAtomCross = NULL;
  // then the EOT-level one (need to be an eoGenOp as global Xover is
  eoGenOp<EOT> *ptCross;

  // check for the atom Xovers
  if (crossObjParam.value() == std::string("discrete"))
    ptObjAtomCross = new eoDoubleExchange;
  else if (crossObjParam.value() == std::string("intermediate"))
    ptObjAtomCross = new eoDoubleIntermediate;
  else if (crossObjParam.value() == std::string("none"))
    ptObjAtomCross = new eoBinCloneOp<double>;
  else throw std::runtime_error("Invalid Object variable crossover type");

  if (crossStdevParam.value() == std::string("discrete"))
    ptStdevAtomCross = new eoDoubleExchange;
  else if (crossStdevParam.value() == std::string("intermediate"))
    ptStdevAtomCross = new eoDoubleIntermediate;
  else if (crossStdevParam.value() == std::string("none"))
    ptStdevAtomCross = new eoBinCloneOp<double>;
  else throw std::runtime_error("Invalid mutation strategy parameter crossover type");

  // and build the indi Xover
  if (crossTypeParam.value() == std::string("global"))
    ptCross = new eoEsGlobalXover<EOT>(*ptObjAtomCross, *ptStdevAtomCross);
  else if (crossTypeParam.value() == std::string("standard"))
    {      // using a standard eoBinOp, but wrap it into an eoGenOp
      eoBinOp<EOT> & crossTmp = _state.storeFunctor(
             new eoEsStandardXover<EOT>(*ptObjAtomCross, *ptStdevAtomCross)
             );
      ptCross = new eoBinGenOp<EOT>(crossTmp);
    }
  else throw std::runtime_error("Invalide Object variable crossover type");

  // now that everything is OK, DON'T FORGET TO STORE MEMORY
  _state.storeFunctor(ptObjAtomCross);
  _state.storeFunctor(ptStdevAtomCross);
  _state.storeFunctor(ptCross);

  //  mutation

  // Ok, time to set up the self-adaptive mutation
  // Proxy for the mutation parameters
  eoEsMutationInit mutateInit(_parser, "Variation Operators");

  eoEsMutate<EOT> & mut =  _state.storeFunctor(
      new eoEsMutate<EOT>(mutateInit, boundsParam.value()));

  // now the general op - a sequential application of crossover and mutatation
  // no need to first have crossover combined with a clone as it is an
  // eoBinOp and not an eoQuadOp as in SGA paradigm

  eoSequentialOp<EOT> & op = _state.storeFunctor(new eoSequentialOp<EOT>);
  op.add(*ptCross, pCrossParam.value());
  op.add(mut, pMutParam.value());

  // that's it!
  return op;
}
#endif // EO_make_op_h


// Local Variables:
// coding: iso-8859-1
// mode:C++
// c-file-style: "Stroustrup"
// comment-column: 35
// fill-column: 80
// End: