/home/kevin/Documents/PYLikeC/fitness_function.c

Go to the documentation of this file.

#include <stdio.h>
#include "currents.h"
#include "pga.h"

#ifdef SINGLE_BUILD
int main(void)
{
   double params[NUM_PARAMS];
   double read_dummy;
#endif
#ifndef SINGLE_BUILD
double fitness_function(double params[])
{
#endif
   struct input_data model_data;
   int exp_delay_count;
   double exp_time[MAX_DAT];
   double exp_delay[MAX_DAT];
   int model_delay_count;
   double model_delay[MAX_DAT];
   double spike_times[MAX_SPIKES];
   int model_spike_count;
   int i;
   int min_delay_count;
   
   model_data.position = 0;
   model_data.filled = 0;
#ifdef SINGLE_BUILD
   scanf("%lf", &read_dummy);
   for (i = 0; i < NUM_PARAMS; i++)
     scanf("%lf %lf", &model_data.params[i], &read_dummy);
#endif
#ifndef SINGLE_BUILD
   for (i = 0; i < NUM_PARAMS; i++)
     model_data.params[i] = params[i];
#endif
   
   fill_struct(&model_data); 
    
   model_spike_count = model(&model_data, spike_times);
   /* removed for all spike runs 
   exp_delay_count = get_exp_delay(exp_time, exp_delay);
   model_delay_count = get_delay(spike_times, model_spike_count, exp_time, exp_delay_count, model_delay);
   if (model_delay_count < exp_delay_count)
     min_delay_count = model_delay_count;
   else
     min_delay_count = exp_delay_count;
   if (min_delay_count > 500)
     return 1 / err(exp_delay, model_delay,  min_delay_count);
   else 
     {
        fprintf(stderr, "too few spikes %i\n", model_delay_count);
        return 0.0;
     }
    */
   
#ifdef SINGLE_BUILD
   fprintf(stderr, "Fit = %f\n", all_spike_fit(spike_times, model_spike_count));
   return 0;
#endif
#ifndef SINGLE_BUILD
   return all_spike_fit(spike_times, model_spike_count);
#endif
}

---