Mutation.hh

// Copyright QED Software 2023.
 
#ifndef GRAIL_MUTATION_H
#define GRAIL_MUTATION_H
 
#include <memory>
#include <random>
#include <vector>
 
namespace grail
{
namespace evolution
{
    class Individual;
    class BaseEvoParam;
 
    struct Mutation
    {
        friend class EAOptimizer;
 
    public:
        enum class MutationSelectionType
        {
            ITERATE_AND_TEST_P,
            SELECT_P_OF_N
        };
        
        Mutation(double mutationPhaseRate = 1.0,
                 double mutateIndividualRate = 0.09,
                 double mutateGeneRate = 0.5,
                 bool mutateBothParentsAndChildren = true,
                 MutationSelectionType randomResolutionTypeForIndividual = MutationSelectionType::SELECT_P_OF_N,
                 MutationSelectionType randomResolutionTypeForGene = MutationSelectionType::SELECT_P_OF_N,
                 std::mt19937_64::result_type seed = std::random_device{}());
        
        virtual ~Mutation();
        
        virtual void Perform(const std::vector<std::unique_ptr<Individual>>& parentPopulation,
                             const std::vector<std::unique_ptr<Individual>>& childrenPopulation);
 
        double MutationPhaseRate = 1.0;
 
        double MutateIndividualRate = 0.09;
 
        double MutateGeneRate = 0.5;
 
        bool MutateBothParentsAndChildren = false;
 
        MutationSelectionType RandomResolutionTypeForIndividual = MutationSelectionType::SELECT_P_OF_N;
 
        MutationSelectionType RandomResolutionTypeForGene = MutationSelectionType::SELECT_P_OF_N;
 
    protected:
        virtual void MutateIndividual(Individual& individual);
 
        virtual void MutateGene(Individual& individual, BaseEvoParam& gene);
 
    protected:
        std::mt19937_64 randGen;
 
    private:
        std::vector<std::unique_ptr<Individual>> populationAfter{};
    };
}
}
 
#endif