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