00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026 #include <algorithm>
00027 #include <utils/eoParser.h>
00028 #include <utils/eoRNG.h>
00029 #include <es/eoReal.h>
00030 #include <utils/eoRealBounds.h>
00031 #include <utils/eoRealVectorBounds.h>
00032
00033
00034
00035 template<class EOT> class eoSBXCrossover: public eoQuadOp<EOT>
00036 {
00037 public:
00038
00039
00040
00041
00042
00043
00044 eoSBXCrossover(const double& _eta = 1.0) :
00045 bounds(eoDummyVectorNoBounds), eta(_eta), range(1) {}
00046
00047
00049
00061 eoSBXCrossover(eoRealVectorBounds & _bounds,
00062 const double& _eta = 1.0) :
00063 bounds(_bounds), eta(_eta), range(1) {}
00064
00066
00068
00075 eoSBXCrossover(eoParser & _parser) :
00076
00077
00078 bounds (_parser.getORcreateParam(eoDummyVectorNoBounds, "objectBounds", "Bounds for variables", 'B', "Variation Operators").value()) ,
00079
00080 eta (_parser.getORcreateParam(1.0, "eta", "SBX eta parameter", '\0', "Variation Operators").value()) ,
00081 range(1) {}
00082
00083
00085 virtual std::string className() const { return "eoSBXCrossover"; }
00086
00087
00088
00089
00090
00091
00092 bool operator()(EOT& _eo1, EOT& _eo2)
00093 {
00094 unsigned i;
00095 double r1, r2, beta;
00096
00097 for (i=0; i<_eo1.size(); i++)
00098 {
00099 double u = rng.uniform(range) ;
00100
00101 if ( u <= 0.5 )
00102 beta = exp( (1/(eta+1))*log(2*u));
00103 else
00104 beta = exp((1/(eta+1))*log(1/(2*(1-u))));
00105
00106
00107
00108 r1=_eo1[i];
00109 r2=_eo2[i];
00110 _eo1[i] =0.5*((1+beta)*r1+(1-beta)*r2);
00111 _eo2[i] =0.5*((1-beta)*r1+(1+beta)*r2);
00112
00113
00114 if(!(bounds.isInBounds(i,_eo1[i])))
00115 bounds.foldsInBounds(i,_eo1[i]);
00116 if(!(bounds.isInBounds(i,_eo2[i])))
00117 bounds.foldsInBounds(i,_eo2[i]);
00118
00119
00120
00121 }
00122 return true;
00123 }
00124
00125
00126
00127 protected:
00128 eoRealVectorBounds & bounds;
00129 double eta;
00130 double range;
00131 };
00132
00133