00001 #include "../ECF_base.h" 00002 #include "Binary.h" 00003 00004 namespace Binary 00005 { 00006 00007 void BinaryCrsRandomRespectful::registerParameters(StateP state) 00008 { 00009 myGenotype_->registerParameter(state, "crx.randomrespectful", (voidP) new double(0), ECF::DOUBLE); 00010 } 00011 00012 00013 bool BinaryCrsRandomRespectful::initialize(StateP state) 00014 { 00015 voidP sptr = myGenotype_->getParameterValue(state, "crx.randomrespectful"); 00016 probability_ = *((double*)sptr.get()); 00017 return true; 00018 } 00019 00020 00021 bool BinaryCrsRandomRespectful::mate(GenotypeP gen1, GenotypeP gen2, GenotypeP child) 00022 { 00023 Binary* p1 = (Binary*) (gen1.get()); 00024 Binary* p2 = (Binary*) (gen2.get()); 00025 Binary* ch = (Binary*) (child.get()); 00026 00027 double rand; 00028 typedef std::vector<uint> v_uint; 00029 00030 v_uint similarity; 00031 00032 for (uint dimension = 0; dimension < p1->variables.size(); dimension++) { 00033 for(uint i = 0; i < p1->getNumBits(); i++) { 00034 if (p1->variables[dimension][i] == p2->variables[dimension][i]) 00035 similarity.push_back(p2->variables[dimension][i]); 00036 else 00037 similarity.push_back(2); 00038 } 00039 } 00040 00041 for (uint dimension = 0; dimension < p1->variables.size(); dimension++) { 00042 for(uint i = 0; i < p1->getNumBits(); i++) { 00043 if (similarity[dimension*p1->getNumBits()+i]==1) 00044 ch->variables[dimension][i]= 1; 00045 else if (similarity[dimension*p1->getNumBits()+i]==0) 00046 ch->variables[dimension][i]= 0; 00047 else { 00048 rand = state_->getRandomizer()->getRandomDouble(); 00049 if (rand <= 0.5) 00050 ch->variables[dimension][i]= 1; 00051 else 00052 ch->variables[dimension][i]= 0; 00053 } 00054 } 00055 } 00056 00057 ch->update(); 00058 00059 return true; 00060 } 00061 00062 }