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

#include "../ECF_base.h"
#include "Binary.h"

namespace Binary
{

void BinaryCrsReducedSurrogate::registerParameters(StateP state)
{
    myGenotype_->registerParameter(state, "crx.reducedsurrogate", (voidP) new double(0), ECF::DOUBLE);
}


bool BinaryCrsReducedSurrogate::initialize(StateP state)
{
    voidP sptr = myGenotype_->getParameterValue(state, "crx.reducedsurrogate");
    probability_ = *((double*)sptr.get());
    return true;
}


bool BinaryCrsReducedSurrogate::mate(GenotypeP gen1, GenotypeP gen2, GenotypeP child)
{
    Binary* p1 = (Binary*) (gen1.get());
    Binary* p2 = (Binary*) (gen2.get());
    Binary* ch = (Binary*) (child.get());
    
    for (uint dimension = 0; dimension < p1->variables.size(); dimension++) {
        
        uint bitCrs = state_->getRandomizer()->getRandomInteger(p1->getNumBits());
        
        //
        // DJ: ovdje ce se zablokirati ako nema razlicitih bitova (npr. za dvije jednake jedinke!)
        //
        while(p1->variables[dimension][bitCrs] != p2->variables[dimension][bitCrs])
        {
            bitCrs = state_->getRandomizer()->getRandomInteger(p1->getNumBits());
        }

        uint parent = state_->getRandomizer()->getRandomInteger(0, 1);
        
        for(uint i = 0; i < p1->getNumBits(); i++) {
            if (p1->variables[dimension][i] == p2->variables[dimension][i]) 
                ch->variables[dimension][i] = p1->variables[dimension][i];
            else    {
                if (parent==0)  {
                    if (i < bitCrs)
                        ch->variables[dimension][i] = p1->variables[dimension][i];
                    else
                        ch->variables[dimension][i] = p2->variables[dimension][i];
                }
                else    {
                    if (i < bitCrs)
                        ch->variables[dimension][i] = p2->variables[dimension][i];
                    else
                        ch->variables[dimension][i] = p1->variables[dimension][i];
                }   
            }
        }
    }

    ch->update();

    return true;
}

}