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

#include "../ECF_base.h"
#include "Binary.h"
#include <cmath>

namespace Binary
{

void BinaryMutSimple::registerParameters(StateP state)
{
    myGenotype_->registerParameter(state, "mut.simple", (voidP) new double(0), ECF::DOUBLE);
    myGenotype_->registerParameter(state, "mut.simple.bitprob", (voidP) new double(0.001), ECF::DOUBLE);
}


bool BinaryMutSimple::initialize(StateP state)
{
    voidP sptr = myGenotype_->getParameterValue(state, "mut.simple");
    probability_ = *((double*)sptr.get());

    sptr = myGenotype_->getParameterValue(state, "mut.simple.bitprob");
    bitProb_ = *((double*)sptr.get());

    bUseBitProb_ = false;
    if(myGenotype_->isParameterDefined(state, "mut.simple.bitprob"))
        bUseBitProb_ = true;

    return true;
}


bool BinaryMutSimple::mutate(GenotypeP gene)
{
    Binary* bin = (Binary*) (gene.get());

    // invert all bits with 'bitProb_' probability
    if(bUseBitProb_) {
        for(uint i = 0; i < bin->variables.size(); i++)
            for(uint j = 0; j < bin->getNumBits(); j++)
                if(state_->getRandomizer()->getRandomDouble() < bitProb_)
                    bin->variables[i][j] = !(bin->variables[i][j]);
    }
    // invert a single random bit in the genotype
    else {
        uint iBit = state_->getRandomizer()->getRandomInteger((uint) bin->getNumBits());
        uint dimension = state_->getRandomizer()->getRandomInteger((uint) bin->variables.size());
        bin->variables[dimension][iBit] = !(bin->variables[dimension][iBit]);
    }

    bin->update();

    return true;
}

}