CartesianCrxOnePoint.cpp

#include "CartesianCrxOnePoint.h"
#include "Cartesian_genotype.h"
namespace cartesian{
 
    void CartesianCrxOnePoint::registerParameters(StateP state)
    {
        myGenotype_->registerParameter(state, "crx.onepoint", (voidP) new double(0), ECF::DOUBLE);
    }
 
    bool CartesianCrxOnePoint::initialize(StateP state)
    {
        voidP sptr = myGenotype_->getParameterValue(state, "crx.onepoint");
        probability_ = *((double*)sptr.get());
        return true;
    }
 
    bool CartesianCrxOnePoint::mate(GenotypeP gen1, GenotypeP gen2, GenotypeP child)
    {
        Cartesian* p1 = (Cartesian*) gen1.get();
        Cartesian* p2 = (Cartesian*) gen2.get();
        Cartesian* ch = (Cartesian*) child.get();
 
        uint parent1Size = p1->size() - p1->nOutputs;
        uint parent2Size = p2->size() - p2->nOutputs;
        uint minimumSize = parent1Size > parent2Size ? parent1Size : parent2Size;
        bool parent1IsLonger = parent1Size > parent2Size;
        ch->clear();
        double whichParent = p1->get_random_double(0,1.0);
        uint index = p1->get_random_int(0, minimumSize);
        if(whichParent < 0.5) {
            for(uint i = 0; i < index; i++) {
                ch->push_back(p1->operator[](i));
            }
            for(uint i = index; i < minimumSize; i++) {
                ch->push_back(p2->operator[](i));
            }
        }
        else {
            for(uint i = 0; i < index; i++) {
                ch->push_back(p2->operator[](i));
            }
            for(uint i = index; i < minimumSize; i++) {
                ch->push_back(p1->operator[](i));
            }
        }
        if(parent1IsLonger) {
            for(uint i = parent2Size; i < parent1Size; i++) {
                ch->push_back(p1->operator[](i));
            }
        }
        else{
            for(uint i = parent1Size; i < parent2Size; i++) {
                ch->push_back(p1->operator[](i));
            }
        }
        for(uint i = 0; i < p1->nOutputs; i++) {
            if(whichParent < 0.5) {
                ch->push_back(p1->operator[](parent1Size + i));
            }
            else {
                ch->push_back(p1->operator[](parent2Size + i));
            }
        }
        return true;
    }
}