#define MAXT 100
#define PI 3.1416
int hist = 2; //filter type --see MakeConvlFunc
int conv_start, conv_tot; /*conv_tot is the total time of the filter and
conv_start is the
beginning point of the filter
with respect to the time of the spike (i.e., it is
zero for the PSP filter, but -3*stddev for a
Gaussian filter */
int binw = 10; //bin width for regular histogram
int stddev = 10; //standard deviation for Gaussian filter
int growth = 1; //growth rate of PSP filter
int decay = 20; //decay rate of PSP filter
float convf[MAXT]; /* array that contains the values of the convolution filter
just make sure big enough for the time length of filter
which is conv_tot and depends on filter variables such as
decay rate */
/***** this function makes the convolution filter **************************/
void MakeConvlFunc(void)
{
register int i;
float temp;
for(i=0; i (TOTTIME-conv_start)) break;
for(j=0; j=0 && (t+j+conv_start)