00001 #include "../ECF_base.h" 00002 #include "FloatingPoint.h" 00003 00004 00005 namespace FloatingPoint 00006 { 00007 00008 void FloatingPointCrsDiscrete::registerParameters(StateP state) 00009 { 00010 myGenotype_->registerParameter(state, "crx.discrete", (voidP) new double(0), ECF::DOUBLE); 00011 } 00012 00013 00014 bool FloatingPointCrsDiscrete::initialize(StateP state) 00015 { 00016 voidP sptr = myGenotype_->getParameterValue(state, "crx.discrete"); 00017 probability_ = *((double*)sptr.get()); 00018 return true; 00019 } 00020 00021 00022 bool FloatingPointCrsDiscrete::mate(GenotypeP gen1, GenotypeP gen2, GenotypeP child) 00023 { 00024 FloatingPoint* p1 = (FloatingPoint*) (gen1.get()); 00025 FloatingPoint* p2 = (FloatingPoint*) (gen2.get()); 00026 FloatingPoint* ch = (FloatingPoint*) (child.get()); 00027 int a; 00028 00029 switch (state_->getRandomizer()->getRandomInteger(0, 1)) { 00030 case 0: for (uint i = 0; i < p1->realValue.size(); i++) { 00031 a = state_->getRandomizer()->getRandomInteger(0, 1); 00032 if (a == 0) 00033 ch->realValue[i] = p1->realValue[i]; 00034 else 00035 ch->realValue[i] = p2->realValue[i]; 00036 } 00037 break; 00038 case 1: for (uint i = 0; i < p1->realValue.size(); i++) { 00039 a = state_->getRandomizer()->getRandomInteger(0, 1); 00040 if (a == 0) 00041 ch->realValue[i] = p2->realValue[i]; 00042 else 00043 ch->realValue[i] = p1->realValue[i]; 00044 } 00045 } 00046 00047 return true; 00048 } 00049 00050 }