IntGenotypeCrxTwoPoint.cpp

#include <ecf/ECF.h>
#include <ecf/ECF_base.h>
#include <cmath>
#include <sstream>
#include "IntGenotype.h"
 
namespace IntGenotype
{
    
        bool IntGenotypeCrxTwoPoint::mate(GenotypeP gen1, GenotypeP gen2, GenotypeP child) {
 
            IntGenotype* p1 = (IntGenotype*)(gen1.get());
            IntGenotype* p2 = (IntGenotype*)(gen2.get());
            IntGenotype* ch = (IntGenotype*)(child.get());
 
            //int prva = (int) p1->getUBound();
            //int druga = (int) p2->getUBound();
            //
            //for (uint i = 0; i < p1->intValues.size(); i++) {
            //  if (p1->intValues[i] < prva)
            //      prva = p1->intValues[i];
            //  if (p2->intValues[i] < druga)
            //      druga = p2->intValues[i];
            //}
 
            //for (uint i = 0; i < p1->intValues.size(); i++) {
            //  p1->intValues[i] = p1->intValues[i] - prva;
            //  p2->intValues[i] = p2->intValues[i] - druga;
            //}
 
            uint dimensionCrs1 = state_->getRandomizer()->getRandomInteger((int)p1->intValues.size());
            uint dimensionCrs2, tmp;
 
            do
            {
                dimensionCrs2 = state_->getRandomizer()->getRandomInteger(p1->intValues.size());
            } while (dimensionCrs1 == dimensionCrs2);
 
            if (dimensionCrs1 > dimensionCrs2)
            {
                tmp = dimensionCrs1;
                dimensionCrs1 = dimensionCrs2;
                dimensionCrs2 = tmp;
            }
 
            switch (state_->getRandomizer()->getRandomInteger(0, 1)) {
            case 0: for (uint i = 0; i < dimensionCrs1; i++) {
                        ch->intValues[i] = p1->intValues[i];
            }
                    for (uint i = dimensionCrs1; i < dimensionCrs2; i++) {
                        ch->intValues[i] = p2->intValues[i];
                    }
                    for (uint i = dimensionCrs2; i < p1->intValues.size(); i++) {
                        ch->intValues[i] = p1->intValues[i];
                    }
                    break;
            case 1: for (uint i = 0; i < dimensionCrs1; i++) {
                        ch->intValues[i] = p2->intValues[i];
            }
                    for (uint i = dimensionCrs1; i < dimensionCrs2; i++) {
                        ch->intValues[i] = p1->intValues[i];
                    }
                    for (uint i = dimensionCrs2; i < p2->intValues.size(); i++) {
                        ch->intValues[i] = p2->intValues[i];
                    }
            }
    
            return true;
        }
    
}