#include <stdio.h>
#include <math.h>
#include "stats.h"
#define N 372            /* number of data points */     
#define ORDER 15         /* size of smoothing sample */

/* General declarations */

double dummy[N];         /* dummy weights series */
int i, j, offset;
double bits_allocated[N], marg_error[N], me_average[N];
double fitted_err[N], bitvec[ORDER], mevec[ORDER];
double slope_0, intercept_0;
FILE *datafile;

int main(void)
{
  datafile = fopen("smooth0.dat", "r");
  if (datafile == NULL) 
    { printf("\nCannot open data file.\n"); exit(1);
    }
  for (i = 0; i < N; i++)
    { fscanf(datafile, "%lf", &bits_allocated[i]);
      fscanf(datafile, "%lf", &marg_error[i]);
      fscanf(datafile, "%lf", &me_average[i]);
    }

  fclose(datafile);
  
  for (j=0;j<ORDER;j++) dummy[j]=1.0/ORDER;
  offset = ORDER - 1;

/* Just checking...
  for (i = 0; i < N; i++)
    { printf("%lf %lf %lf\n", bits_allocated[i], marg_error[i], me_average[i]);
    }  
*/

  for (i = offset; i < N; i++)
    { for (j = 0; j < ORDER; j++)
        { bitvec[j] = bits_allocated[i-j];
	  mevec[j] = marg_error[i-j];
	}
       slope_0 = slope(bitvec, mevec, dummy, ORDER);
       intercept_0 = intercept(bitvec, mevec, dummy, ORDER);
       /* printf("%f\t%f\n", slope_0, intercept_0); */
       fitted_err[i] = intercept_0 + slope_0 * bits_allocated[i];
       printf("%f %f %f\n", 
          bits_allocated[i], marg_error[i], fitted_err[i]);
    }
  
  return 0;
  
 }