/home/kevin/Documents/PYLikeC/mutation.c

Go to the documentation of this file.

#include <math.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include "pga.h"

void mutation(struct individual parents[], struct individual children[], int gen, gsl_rng * r)
{
   double tau;
   double tau_prime;
   int i, j;
   
   tau       = 1.0 / sqrt( 2 * sqrt(NUM_PARAMS));
   tau_prime = 1.0 / sqrt( 2 * NUM_PARAMS) * gsl_ran_gaussian(r, 1.0);
   for (i = 0; i < NUM_PARAMS; i++)
     for (j = 0; j < GEN_SIZE; j++)
       {
          
          children[j].std_dev[i] = parents[j].std_dev[i] * exp(tau_prime + tau * gsl_ran_gaussian(r, 1.0));
          children[j].params[i] = parents[j].params[i] + children[j].std_dev[i] * gsl_ran_cauchy(r, 1.0);
          children[j].created = gen;
          /* Copy fitness for recombination choosiness */
          children[j].fitness = parents[j].fitness;
         
          /*
          children[j].std_dev[i] = parents[j].std_dev[i] + parents[j].std_dev[i] * 0.5 * g05ddc(0.0, 1.0);
          children[j].params[i] = parents[j].params[i] + children[j].std_dev[i] * g05ddc(0.0, 1.0);
          */
          if (children[j].params[i] < 0)
            children[j].params[i] = 0;
       }
}

        

---