/home/kevin/Documents/PYLikeC/pga.h

Go to the documentation of this file.

#ifndef PGA_H
#define PGA_H

#include <gsl/gsl_rng.h>

#define NUM_PARAMS 13
#define POP_SIZE 200 /* fits parents and children*/
#define GEN_SIZE 100
#define MAX_PROCS 32 /* has no effect on spga */

#define GIMME_CHILD 1
#define GOT_FIT 2
#define CHILD_UP 3
#define CHILD_DOWN 4



struct individual
{
   double params[NUM_PARAMS];
   double std_dev[NUM_PARAMS];
   double fitness;
   int wins;
   int created;
   int id;
};

double err(double exp_data[], double modl_data[], size_t length);

double fitness_function(double params[]);

void initialize_pop(struct individual parents[], double p_min[], double p_max[], double initial_std_dev, gsl_rng *r);

void mutation(struct individual parents[], struct individual children[], int gen, gsl_rng *r);

void selection(struct individual parents[], struct individual children[], int competition_size, gsl_rng *r);

void independent_assortment(struct individual children[], gsl_rng * r);

int compare_individuals(const void *a, const void *b);

void local_search(struct individual *child, int rank);

double all_spike_fit(double modl_data[], int length);

double variable_winsize_acuracy(double modl_data[], int mdl_length, double exp_data[], int exp_length);

double matrix_fit(double modl_data[], size_t mdl_length, double exp_data[], size_t exp_length);

double sum_nearest_err(
                double modl_data[], size_t mdl_length,
                double exp_data[], size_t exp_length);

double spike_replace_fit(double modl_data[], int mdl_length, double exp_data[], int exp_length);

int load_pop(struct individual parents[]);

int load_bin_pop(int *gen, struct individual parents[], struct individual children[], int done_children[]);

void save_bin_pop(int *gen, struct individual parents[], struct individual children[], int done_children[]);

#endif

---