---

mutator.h

Go to the documentation of this file.

#ifndef __MUTATOR_H__
#define __MUTATOR_H__

#include <magic/mmath.h>

class FloatMutator {
  public:
    virtual double  mutate  (double x, double min, double max, double variance) const{FORBIDDEN}
};

class ConstantMutator : public FloatMutator {
    double mVariance;
  public:
                    ConstantMutator (double var=0.22) : mVariance(var){}

    virtual double  mutate  (double p, double min, double max, double variance) const {
        //TRACE2 ("p=%f, v=%f", p, variance);
        double pp=p;
        if (max>min) { // Can be 0 -> gene is immutable
            // Default mutation
            double delta;
            do {
                delta = gaussrnd (mVariance);
            } while (pp+delta<min || pp+delta>max);
            pp += delta;
        }
        return pp;
    }
};

class MutationMutator : public FloatMutator {
    double mPhi;
  public:
                    MutationMutator (double phi=0.22) : mPhi(phi){}
    
    virtual double  mutate  (double p, double min, double max, double variance) const {
        //TRACE2 ("p=%f, v=%f", p, variance);
        double pp = 1/(1+(1-p)/p*exp(-gaussrnd(mPhi)));
        if (pp<min)
            pp=min;
        return pp;
    }
};

class MutationSimpleMutator : public FloatMutator {
    double mTau;
  public:
                    MutationSimpleMutator   (int n=100) : mTau(1/sqrt(double(n))){}
    
    virtual double  mutate  (double p, double min, double max, double variance) const {
        //TRACE2 ("p=%f, v=%f", p, variance);
        double pp = p*exp(mTau*gaussrnd(1));
        if (pp<min)
            pp=min;
        return pp;
    }
};


class CircularMutator : public FloatMutator {
    double mPhi;
  public:
    virtual double  mutate  (double p, double min, double max, double variance) const {
        return 1/(1+(1-p)/p*exp(-mPhi*gaussrnd(1)));
    }
};

#endif

---