BlendCrossover.h

Go to the documentation of this file.

#include "boost/smart_ptr.hpp"
#include "Util/CmdArguments.h"
#include <cmath>

template <class T> struct BlendCrossover
{
    void crossover(boost::shared_ptr<T> ind1, boost::shared_ptr<T> ind2,
                   CmdArguments_t args) {

        // BLX-0.5 performs better than BLX operators with any other \alpha
        // value
        const double alpha = 0.5;

        for(size_t i = 0; i < ind1->genes_m.size(); i++) {

            double ming = std::min(ind1->genes_m[i], ind2->genes_m[i]);
            double maxg = std::max(ind1->genes_m[i], ind2->genes_m[i]);
            double gamma1 = (1 + 2 * alpha) *
                static_cast<double>(rand() / (RAND_MAX + 1.0)) - alpha;
            double gamma2 = (1 + 2 * alpha) *
                static_cast<double>(rand() / (RAND_MAX + 1.0)) - alpha;
            ind1->genes_m[i] = (1 - gamma1) * ming + gamma1 * maxg;
            ind2->genes_m[i] = (1 - gamma2) * ming + gamma2 * maxg;
        }
    }
};