CartesianCrxUniform.cpp

#include "CartesianCrxUniform.h"
#include "Cartesian_genotype.h"
 
namespace cartesian{
    void CartesianCrxUniform::registerParameters(StateP state)
    {
        myGenotype_->registerParameter(state, "crx.uniform", (voidP) new double(0), ECF::DOUBLE);
    }
 
    bool CartesianCrxUniform::initialize(StateP state)
    {
        voidP sptr = myGenotype_->getParameterValue(state, "crx.uniform");
        probability_ = *((double*)sptr.get());
        return true;
    }
 
    bool CartesianCrxUniform::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();
 
        for(uint i = 0; i < minimumSize; i++) {
            double randomResult = p1->get_random_double(0,1);
            if(randomResult < 0.5) {
                ch->push_back(p1->operator[](i));
            }
            else{
                ch->push_back(p2->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));
            }
        }
        //Resolve outputs at random.
        for(uint i = 0; i < p1->nOutputs; i++) {
            double randomResult = p1->get_random_double(0,1);
            if(randomResult < 0.5) {
                ch->push_back(p1->operator[](parent1Size + i));
            }
            else {
                ch->push_back(p1->operator[](parent2Size + i));
            }
        }
        return true;
    }
}