• Main Page
  • Modules
  • Classes
  • Files
  • File List

D:/Projekt/ECF_trunk/ECF/binary/BinaryCrsTwoPoint.cpp

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     // update integer and real domain representation
00066     ch->update();
00067 
00068     return true;
00069 }
00070 
00071 }

Generated on Tue Nov 4 2014 13:04:30 for ECF by  doxygen 1.7.1