/home/kevin/Documents/PYLikeC/all_spike_fit.c

Go to the documentation of this file.

#define _GNU_SOURCE
#include <math.h>
#include <stdio.h>
#include "pga.h"
#include "currents.h"

double all_spike_fit(double modl_data[], int length)
{
   FILE *spike_file_ptr;
   double exp_data[MAX_SPIKES];
   int count;
   double difference;
   int i, j, done;
   
   if (length < 2 )
     {
        fprintf(stderr, "too few spikes %i\n", length);
        return -INFINITY;
     }
   else
     {
        
//      spike_file_ptr = fopen("250_upt_all_spikes.txt", "r");
        
//      spike_file_ptr = fopen("model_spike_times.txt", "r");

        spike_file_ptr = fopen("model_spike_times_random.txt", "r");
        
        count = 0;
        while((fscanf(spike_file_ptr, "%lf", &exp_data[count]) == 1) && (count < MAX_SPIKES))
          count++;
        fclose(spike_file_ptr);
/*      i = 0;
        j = 0;
        done = 0;
        difference = 0;
        while (done == 0)
          {
             difference += fabs(modl_data[i] - exp_data[j]);
             if (i + 1 < length)
               i++;
             if (j + 1 < count)
               j++;
             if ((i + 1 >= length) && (j + 1 >= count))
               done = 1;
          }
        return -difference * abs(length - count);
 */
        /* return matrix_fit(modl_data, length, exp_data, count); */
        return spike_replace_fit(modl_data, length, exp_data, count);
     }
}

---