00001 #include "../ECF_base.h"
00002 #include "Binary.h"
00003
00004 namespace Binary
00005 {
00006
00007 void BinaryCrsMasked::registerParameters(StateP state)
00008 {
00009 myGenotype_->registerParameter(state, "crx.masked", (voidP) new double(0), ECF::DOUBLE);
00010 }
00011
00012
00013 bool BinaryCrsMasked::initialize(StateP state)
00014 {
00015 voidP sptr = myGenotype_->getParameterValue(state, "crx.masked");
00016 probability_ = *((double*)sptr.get());
00017 return true;
00018 }
00019
00020
00021 bool BinaryCrsMasked::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 uint rand;
00028
00029 rand = state_->getRandomizer()->getRandomInteger(0, 1);
00030
00031 for (uint dimension = 0; dimension < p1->variables.size(); dimension++) {
00032 uint i;
00033 switch (rand) {
00034 case 0: for (i = 0; i < p1->getNumBits(); i++) {
00035 ch->variables[dimension][i] = p1->variables[dimension][i];
00036 }
00037 break;
00038 case 1: for (i = 0; i < p2->getNumBits(); i++) {
00039 ch->variables[dimension][i] = p2->variables[dimension][i];
00040 }
00041 }
00042 }
00043
00044 for (uint dimension = 0; dimension < p1->variables.size(); dimension++) {
00045 uint i;
00046 switch (rand) {
00047 case 0: for (i = 0; i < p1->getNumBits(); i++) {
00048 if (p2->variables[dimension][i] == 1 && p1->variables[dimension][i] == 0)
00049 ch->variables[dimension][i] = p2->variables[dimension][i];
00050 }
00051 break;
00052 case 1: for (i = 0; i < p2->getNumBits(); i++) {
00053 if (p2->variables[dimension][i] == 0 && p1->variables[dimension][i] == 1)
00054 ch->variables[dimension][i] = p1->variables[dimension][i];
00055 }
00056 }
00057 }
00058
00059
00060 ch->update();
00061
00062 return true;
00063 }
00064
00065 }