10064ad967
src/es/eoEsGlobalXover.h and src/es/eoEsLocalXover.h for crossover src/es/make_XXX_es.h for user-input test/t-eoEsAll.cpp to test However, an old bug appeared: className was not const in eoGenOp (and derived classes) so I had to change that throughtout the hierarchy
139 lines
4.3 KiB
C++
139 lines
4.3 KiB
C++
/** -*- mode: c++; c-indent-level: 4; c++-member-init-indent: 8; comment-column: 35; -*-
|
|
|
|
//-----------------------------------------------------------------------------
|
|
// eoEsGlobalXover.h : ES global crossover
|
|
// (c) Marc Schoenauer 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: marc.schoenauer@polytechnique.fr http://eeaax.cmap.polytchnique.fr/
|
|
*/
|
|
//-----------------------------------------------------------------------------
|
|
|
|
|
|
#ifndef _eoEsGlobalXover_H
|
|
#define _eoEsGlobalXover_H
|
|
|
|
#include <utils/eoRNG.h>
|
|
|
|
#include <es/eoEsSimple.h>
|
|
#include <es/eoEsStdev.h>
|
|
#include <es/eoEsFull.h>
|
|
|
|
#include <eoGenOp.h>
|
|
// needs a selector - here random
|
|
#include <eoRandomSelect.h>
|
|
|
|
/** Gloabl crossover operator for ES genotypes.
|
|
* Uses some Atom crossovers to handle both the object variables
|
|
* and the mutation strategy parameters
|
|
*/
|
|
template<class EOT>
|
|
class eoEsGlobalXover: public eoGenOp<EOT>
|
|
{
|
|
public:
|
|
typedef typename EOT::Fitness FitT;
|
|
|
|
/**
|
|
* (Default) Constructor.
|
|
*/
|
|
eoEsGlobalXover(eoBinOp<double> & _crossObj, eoBinOp<double> & _crossMut) :
|
|
crossObj(_crossObj), crossMut(_crossMut) {}
|
|
|
|
/// The class name. Used to display statistics
|
|
virtual string className() const { return "eoEsGlobalXover"; }
|
|
|
|
/// The TOTAL number of offspring (here = nb of parents modified in place)
|
|
unsigned max_production(void) { return 1; }
|
|
|
|
/**
|
|
* modifies one parents in the populator
|
|
* using 2 new parents for each component!
|
|
*
|
|
* @param _pop a POPULATOR (not a simple population)
|
|
*/
|
|
void apply(eoPopulator<EOT>& _plop)
|
|
{
|
|
// First, select as many parents as you will have offspring
|
|
EOT& parent = *_plop; // select the first parent
|
|
|
|
// first, the object variables
|
|
for (unsigned i=0; i<parent.size(); i++)
|
|
{
|
|
// get extra parents - use private selector
|
|
// _plop.source() is the eoPop<EOT> used by _plop to get parents
|
|
const EOT& realParent1 = sel(_plop.source());
|
|
const EOT& realParent2 = sel(_plop.source());
|
|
parent[i] = realParent1[i];
|
|
crossObj(parent[i], realParent2[i]); // apply eoBinOp
|
|
}
|
|
// then the self-adaptation parameters
|
|
cross_self_adapt(parent, _plop.source());
|
|
// dont' forget to invalidate
|
|
parent.invalidate();
|
|
}
|
|
|
|
private:
|
|
|
|
// the method to cross slef-adaptation parameters: need to specialize
|
|
|
|
void cross_self_adapt(eoEsSimple<FitT> & _parent, const eoPop<eoEsSimple<FitT> >& _pop)
|
|
{
|
|
const EOT& realParent1 = sel(_pop);
|
|
const EOT& realParent2 = sel(_pop);
|
|
_parent.stdev = realParent1.stdev;
|
|
crossMut(_parent.stdev, realParent2.stdev); // apply eoBinOp
|
|
}
|
|
|
|
void cross_self_adapt(eoEsStdev<FitT> & _parent, const eoPop<eoEsStdev<FitT> >& _pop)
|
|
{
|
|
for (unsigned i=0; i<_parent.size(); i++)
|
|
{
|
|
const EOT& realParent1 = sel(_pop);
|
|
const EOT& realParent2 = sel(_pop);
|
|
_parent.stdevs[i] = realParent1.stdevs[i];
|
|
crossMut(_parent.stdevs[i], realParent2.stdevs[i]); // apply eoBinOp
|
|
}
|
|
}
|
|
|
|
void cross_self_adapt(eoEsFull<FitT> & _parent, const eoPop<eoEsFull<FitT> >& _pop)
|
|
{
|
|
unsigned i;
|
|
// the StDev
|
|
for (i=0; i<_parent.size(); i++)
|
|
{
|
|
const EOT& realParent1 = sel(_pop);
|
|
const EOT& realParent2 = sel(_pop);
|
|
_parent.stdevs[i] = realParent1.stdevs[i];
|
|
crossMut(_parent.stdevs[i], realParent2.stdevs[i]); // apply eoBinOp
|
|
}
|
|
// the roataion angles
|
|
for (i=0; i<_parent.correlations.size(); i++)
|
|
{
|
|
const EOT& realParent1 = sel(_pop);
|
|
const EOT& realParent2 = sel(_pop);
|
|
_parent.correlations[i] = realParent1.correlations[i];
|
|
crossMut(_parent.correlations[i], realParent2.correlations[i]); // apply eoBinOp
|
|
}
|
|
|
|
}
|
|
|
|
// the data
|
|
eoRandomSelect<EOT> sel;
|
|
eoBinOp<double> & crossObj;
|
|
eoBinOp<double> & crossMut;
|
|
};
|
|
|
|
#endif
|