Very simple classes for vectors of (unbounded) real numbers.
A few straightforward operators only (BLX-0, uniform mutation). Mostly for tutorial reasons - though it might improve in usefull classes later
This commit is contained in:
evomarc
parent
2b8c049283
commit
d0e566fdf4
2 changed files with 330 additions and 0 deletions
58
eo/src/es/eoReal.h
Normal file
58
eo/src/es/eoReal.h
Normal file
|
|
@ -0,0 +1,58 @@
|
|||
/*
|
||||
eoReal.h
|
||||
// (c) Marc Schoenauer, Maarten Keijzer and GeNeura Team, 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: Marc.Schoenauer@polytechnique.fr
|
||||
todos@geneura.ugr.es, http://geneura.ugr.es
|
||||
mkeijzer@dhi.dk
|
||||
*/
|
||||
|
||||
#ifndef eoReal_h
|
||||
#define eoReal_h
|
||||
|
||||
//-----------------------------------------------------------------------------
|
||||
|
||||
#include <iostream> // ostream, istream
|
||||
#include <string> // string
|
||||
|
||||
#include <eoFixedLength.h>
|
||||
|
||||
/** eoReal: implementation of simple real-valued chromosome.
|
||||
* based on eoFixedLength class
|
||||
*/
|
||||
template <class F> class eoReal: public eoFixedLength<F, double>
|
||||
{
|
||||
public:
|
||||
|
||||
/**
|
||||
* (Default) Constructor.
|
||||
* @param size Size of the vector
|
||||
*/
|
||||
eoReal(unsigned size = 0, double value = 0.0):
|
||||
eoFixedLength<F, double>(size, value) {}
|
||||
|
||||
/// My class name.
|
||||
string className() const
|
||||
{
|
||||
return "eoReal";
|
||||
}
|
||||
|
||||
};
|
||||
|
||||
//-----------------------------------------------------------------------------
|
||||
|
||||
#endif //eoReal_h
|
||||
272
eo/src/es/eoRealOp.h
Normal file
272
eo/src/es/eoRealOp.h
Normal file
|
|
@ -0,0 +1,272 @@
|
|||
// -*- mode: c++; c-indent-level: 4; c++-member-init-indent: 8; comment-column: 35; -*-
|
||||
|
||||
//-----------------------------------------------------------------------------
|
||||
// eoRealOp.h
|
||||
// (c) EEAAX 2000 - Maarten Keijzer 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: Marc.Schoenauer@polytechnique.fr
|
||||
mak@dhi.dk
|
||||
*/
|
||||
//-----------------------------------------------------------------------------
|
||||
|
||||
#ifndef eoRealOp_h
|
||||
#define eoRealOp_h
|
||||
|
||||
//-----------------------------------------------------------------------------
|
||||
|
||||
#include <algorithm> // swap_ranges
|
||||
#include <utils/eoRNG.h>
|
||||
#include <es/eoReal.h>
|
||||
|
||||
//-----------------------------------------------------------------------------
|
||||
|
||||
/** eoUniformMutation --> changes all values of the vector
|
||||
by uniform choice with range epsilon
|
||||
with probability p_change per variable
|
||||
\class eoUniformMutation eoRealOp.h Tutorial/eoRealOp.h
|
||||
\ingroup parameteric
|
||||
*/
|
||||
|
||||
template<class Chrom> class eoUniformMutation: public eoMonOp<Chrom>
|
||||
{
|
||||
public:
|
||||
/**
|
||||
* (Default) Constructor.
|
||||
* @param _epsilon the range for uniform nutation
|
||||
* @param _p_change the probability to change a given coordinate
|
||||
*/
|
||||
eoUniformMutation(const double& _epsilon, const double& _p_change = 1.0):
|
||||
epsilon(_epsilon), p_change(_p_change) {}
|
||||
|
||||
/// The class name.
|
||||
string className() const { return "eoUniformMutation"; }
|
||||
|
||||
/**
|
||||
* Do it!
|
||||
* @param chrom The cromosome undergoing the mutation
|
||||
*/
|
||||
void operator()(Chrom& chrom)
|
||||
{
|
||||
chrom.invalidate();
|
||||
for (unsigned lieu=0; lieu<chrom.size(); lieu++)
|
||||
{
|
||||
if (rng.flip(p_change))
|
||||
{
|
||||
chrom[lieu] += 2*epsilon*rng.uniform()-epsilon;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
private:
|
||||
double epsilon;
|
||||
double p_change;
|
||||
};
|
||||
|
||||
/** eoDetUniformMutation --> changes exactly k values of the vector
|
||||
by uniform choice with range epsilon
|
||||
\class eoDetUniformMutation eoRealOp.h Tutorial/eoRealOp.h
|
||||
\ingroup parameteric
|
||||
*/
|
||||
|
||||
template<class Chrom> class eoDetUniformMutation: public eoMonOp<Chrom>
|
||||
{
|
||||
public:
|
||||
/**
|
||||
* (Default) Constructor.
|
||||
* @param _epsilon the range for uniform nutation
|
||||
* @param number of coordinate to modify
|
||||
*/
|
||||
eoDetUniformMutation(const double& _epsilon, const unsigned& _no = 1):
|
||||
epsilon(_epsilon), no(_no) {}
|
||||
|
||||
/// The class name.
|
||||
string className() const { return "eoDetUniformMutation"; }
|
||||
|
||||
/**
|
||||
* Do it!
|
||||
* @param chrom The cromosome undergoing the mutation
|
||||
*/
|
||||
void operator()(Chrom& chrom)
|
||||
{
|
||||
chrom.invalidate();
|
||||
for (unsigned i=0; i<no; i++)
|
||||
{
|
||||
unsigned lieu = rng.random(chrom.size());
|
||||
// actually, we should test that we don't re-modify same variable!
|
||||
chrom[lieu] += 2*epsilon*rng.uniform()-epsilon;
|
||||
}
|
||||
}
|
||||
|
||||
private:
|
||||
double epsilon;
|
||||
unsigned no;
|
||||
};
|
||||
|
||||
|
||||
// two arithmetical crossovers
|
||||
|
||||
/** eoSegmentCrossover --> uniform choice in segment
|
||||
== arithmetical with same value along all coordinates
|
||||
\class eoSegmentCrossover eoRealOp.h Tutorial/eoRealOp.h
|
||||
\ingroup parameteric
|
||||
*/
|
||||
|
||||
template<class Chrom> class eoSegmentCrossover: public eoQuadraticOp<Chrom>
|
||||
{
|
||||
public:
|
||||
/**
|
||||
* (Default) Constructor.
|
||||
* @param _alpha the amount of exploration OUTSIDE the parents
|
||||
* as in BLX-alpha notation (Eshelman and Schaffer)
|
||||
* 0 == contractive application
|
||||
*/
|
||||
eoSegmentCrossover(const double& _alpha = 0.0) :
|
||||
alpha(_alpha), range(1+2*alpha) {}
|
||||
|
||||
/// The class name.
|
||||
string className() const { return "eoSegmentCrossover"; }
|
||||
|
||||
/**
|
||||
* segment crossover - modifies both parents
|
||||
* @param chrom1 The first parent
|
||||
* @param chrom2 The first parent
|
||||
*/
|
||||
void operator()(Chrom& chrom1, Chrom& chrom2)
|
||||
{
|
||||
unsigned i;
|
||||
double r1, r2, fact;
|
||||
fact = rng.uniform(range); // in [0,range)
|
||||
for (i=0; i<chrom1.size(); i++)
|
||||
{
|
||||
r1=chrom1[i];
|
||||
r2=chrom2[i];
|
||||
chrom1[i] = fact * r1 + (1-fact) * r2;
|
||||
chrom2[i] = (1-fact) * r1 + fact * r2;
|
||||
}
|
||||
chrom1.invalidate(); // shoudl test if fact was 0 or 1 :-)))
|
||||
chrom2.invalidate();
|
||||
}
|
||||
|
||||
protected:
|
||||
double alpha;
|
||||
double range; // == 1+2*alpha
|
||||
};
|
||||
|
||||
/** eoArithmeticCrossover --> uniform choice in hypercube
|
||||
== arithmetical with different values for each coordinate
|
||||
\class eoArithmeticCrossover eoRealOp.h Tutorial/eoRealOp.h
|
||||
\ingroup parameteric
|
||||
*/
|
||||
|
||||
template<class Chrom> class eoArithmeticCrossover: public eoQuadraticOp<Chrom>
|
||||
{
|
||||
public:
|
||||
/**
|
||||
* (Default) Constructor.
|
||||
* @param _alpha the amount of exploration OUTSIDE the parents
|
||||
* as in BLX-alpha notation (Eshelman and Schaffer)
|
||||
* 0 == contractive application
|
||||
*/
|
||||
eoArithmeticCrossover(const double& _alpha = 0.0):
|
||||
alpha(_alpha), range(1+2*alpha) {}
|
||||
|
||||
/// The class name.
|
||||
string className() const { return "eoArithmeticCrossover"; }
|
||||
|
||||
/**
|
||||
* arithmetical crossover - modifies both parents
|
||||
* @param chrom1 The first parent
|
||||
* @param chrom2 The first parent
|
||||
*/
|
||||
void operator()(Chrom& chrom1, Chrom& chrom2)
|
||||
{
|
||||
unsigned i;
|
||||
double r1, r2, fact;
|
||||
for (i=0; i<chrom1.size(); i++)
|
||||
{
|
||||
r1=chrom1[i];
|
||||
r2=chrom2[i];
|
||||
fact = rng.uniform(range); // in [0,range)
|
||||
chrom1[i] = fact * r1 + (1-fact) * r2;
|
||||
chrom2[i] = (1-fact) * r1 + fact * r2;
|
||||
}
|
||||
}
|
||||
|
||||
protected:
|
||||
double alpha;
|
||||
double range; // == 1+2*alpha
|
||||
};
|
||||
|
||||
|
||||
/** eoRealUxOver --> Uniform crossover, also termed intermediate crossover
|
||||
\class eoRealUxOver eoRealOp.h Tutorial/eoRealOp.h
|
||||
\ingroup parameteric
|
||||
*/
|
||||
|
||||
template<class Chrom> class eoRealUxOver: public eoQuadraticOp<Chrom>
|
||||
{
|
||||
public:
|
||||
/**
|
||||
* (Default) Constructor.
|
||||
* @param _preference bias in the choice (usually, no bias == 0.5)
|
||||
*/
|
||||
eoRealUxOver(const float& _preference = 0.5): preference(_preference)
|
||||
{
|
||||
if ( (_preference <= 0.0) || (_preference >= 1.0) )
|
||||
runtime_error("UxOver --> invalid preference");
|
||||
}
|
||||
|
||||
/// The class name.
|
||||
string className() const { return "eoRealUxOver"; }
|
||||
|
||||
/**
|
||||
* Uniform crossover for real vectors
|
||||
* @param chrom1 The first parent
|
||||
* @param chrom2 The second parent
|
||||
* @runtime_error if sizes don't match
|
||||
*/
|
||||
void operator()(Chrom& chrom1, Chrom& chrom2)
|
||||
{
|
||||
if ( chrom1.size() != chrom2.size())
|
||||
runtime_error("UxOver --> chromosomes sizes don't match" );
|
||||
bool changed = false;
|
||||
for (unsigned int i=0; i<chrom1.size(); i++)
|
||||
{
|
||||
if (rng.flip(preference))
|
||||
if (chrom1[i] == chrom2[i])
|
||||
{
|
||||
double tmp = chrom1[i];
|
||||
chrom1[i]=chrom2[i];
|
||||
chrom2[i] = tmp;
|
||||
changed = true;
|
||||
}
|
||||
}
|
||||
if (changed)
|
||||
{
|
||||
chrom1.invalidate();
|
||||
chrom2.invalidate();
|
||||
}
|
||||
}
|
||||
private:
|
||||
float preference;
|
||||
};
|
||||
|
||||
|
||||
//-----------------------------------------------------------------------------
|
||||
//@}
|
||||
#endif eoRealOp_h
|
||||
|
||||
Reference in a new issue