eoEsGlobalXover.h

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


#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>

template<class EOT>
class eoEsGlobalXover: public eoGenOp<EOT>
{
public:
  typedef typename EOT::Fitness FitT;

  eoEsGlobalXover(eoBinOp<double> & _crossObj, eoBinOp<double> & _crossMut) :
    crossObj(_crossObj), crossMut(_crossMut) {}

  virtual std::string className() const { return "eoEsGlobalXover"; }

  unsigned max_production(void) { return 1; }

  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:

  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