eoProportionalCombinedOp.h

// -*- mode: c++; c-indent-level: 4; c++-member-init-indent: 8; comment-column: 35; -*-
//-----------------------------------------------------------------------------
// eoCombinedOp.h
// (c) GeNeura Team, 1998, Marc Schoenauer, 2000
/* 
    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
 */
//-----------------------------------------------------------------------------

#ifndef _eoCombinedOp_H
#define _eoCombinedOp_H

#include <eoObject.h>
#include <eoPrintable.h>
#include <eoFunctor.h>
#include <eoOp.h>
#include <utils/eoRNG.h>




template <class EOT>
class eoPropCombinedMonOp: public eoMonOp<EOT>
{
public:
  eoPropCombinedMonOp(eoMonOp<EOT> & _first, const double _rate)
  {
    ops.push_back(&_first);
    rates.push_back(_rate);
  }

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

  virtual void add(eoMonOp<EOT> & _op, const double _rate, bool _verbose=false)
  {
    ops.push_back(&_op);
    rates.push_back(_rate);
    // compute the relative rates in percent - to warn the user!
    if (_verbose)
      printOn(std::cout);
  }

  // outputs the operators and percentages
  virtual void printOn(std::ostream & _os)
  {
    double total = 0;
    unsigned i;
    for (i=0; i<ops.size(); i++)
      total += rates[i];
    _os << "In " << className() << "\n" ;
    for (i=0; i<ops.size(); i++)
      _os << ops[i]->className() << " with rate " << 100*rates[i]/total << " %\n";
  }

  virtual bool operator()(EOT & _indi)
  {
    unsigned what = rng.roulette_wheel(rates); // choose one op
    return (*ops[what])(_indi);            // apply it
  }
protected:
std::vector<eoMonOp<EOT>*> ops;
std::vector<double> rates;
};


template <class EOT>
class eoPropCombinedBinOp: public eoBinOp<EOT>
{
public:
  eoPropCombinedBinOp(eoBinOp<EOT> & _first, const double _rate)
  {
    ops.push_back(&_first);
    rates.push_back(_rate);
  }

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

virtual void add(eoBinOp<EOT> & _op, const double _rate, bool _verbose=false)
  {
    ops.push_back(&_op);
    rates.push_back(_rate);
    // compute the relative rates in percent - to warn the user!
    if (_verbose)
      {
        double total = 0;
        unsigned i;
        for (i=0; i<ops.size(); i++)
          total += rates[i];
        std::cout << "In " << className() << "\n" ;
        for (i=0; i<ops.size(); i++)
          std::cout << ops[i]->className() << " with rate " << 100*rates[i]/total << " %\n";
      }
  }

  virtual void operator()(EOT & _indi1, const EOT & _indi2)
  {
    unsigned what = rng.roulette_wheel(rates); // choose one op index
    return (*ops[what])(_indi1, _indi2);                   // apply it
  }
private:
std::vector<eoBinOp<EOT>*> ops;
std::vector<double> rates;
};



template <class EOT>
class eoPropCombinedQuadOp: public eoQuadOp<EOT>
{
public:
  eoPropCombinedQuadOp(eoQuadOp<EOT> & _first, const double _rate)
  {
    ops.push_back(&_first);
    rates.push_back(_rate);
  }

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

  // addition of a true operator
virtual void add(eoQuadOp<EOT> & _op, const double _rate, bool _verbose=false)
  {
    ops.push_back(&_op);
    rates.push_back(_rate);
    // compute the relative rates in percent - to warn the user!
    if (_verbose)
      printOn(std::cout);
  }

  // outputs the operators and percentages
  virtual void printOn(std::ostream & _os)
  {
    double total = 0;
    unsigned i;
    for (i=0; i<ops.size(); i++)
      total += rates[i];
    _os << "In " << className() << "\n" ;
    for (i=0; i<ops.size(); i++)
      _os << ops[i]->className() << " with rate " << 100*rates[i]/total << " %\n";
  }

  virtual bool operator()(EOT & _indi1, EOT & _indi2)
  {
    unsigned what = rng.roulette_wheel(rates); // choose one op index
    return (*ops[what])(_indi1, _indi2);                   // apply it
  }
private:
std::vector<eoQuadOp<EOT>*> ops;
std::vector<double> rates;
};


// for General Ops, it's another story - see eoOpCOntainer
#endif

---