#include <stdio.h>
#include <stdlib.h>
#include "timing.h"
#include <math.h>

int main(int argc, char** argv) {

  double wct_start,wct_end;
  int id,nt,k,i;
  int NITER; 
  int N;
  double t = 0;
  double flops;
  

  if(argc != 2) {
	printf("Usage: %s <size>\n",argv[0]);
	exit(1);
  }
  N = atoi(argv[1]);
  int size = (int)(pow(1.5, N));

  double *a = (double*)malloc(size*sizeof(double));
  double *b = (double*)malloc(size*sizeof(double));
  double *d = (double*)malloc(size*sizeof(double));
  for(k=0; k<N; ++k) {
	a[k]=0.0;
        b[k] = 1.0;
        d[k] = 1.0;  
  }
 const double s = 1.00000000001;

  NITER=1;
  do {
 wct_start = getTimeStamp();

	// repeat measurement often enough
	for(k=0; k<NITER; ++k) {
  	// This is the benchmark loop
  	for(i=0; i<(int)(pow(1.5,N)); ++i) 
       {
    	  //a[i] = b[i] + s * d[i];
         a[i] = s * b[i] + a[i];
 	}
  	// end of benchmark loop
  	if(a[(int)(pow(1.5,N))/2]<0.) printf("%lf",a[(int)(pow(1.5,N))/2]);
	}
	wct_end = getTimeStamp();
	NITER = NITER*2;
  } while (wct_end-wct_start<0.1); // at least 100 ms

  NITER = NITER/2;
  //Stream flops
  flops = NITER * 2 / (wct_end-wct_start);

  // Adapt this to print any other interesting metrics, e.g., performance
  printf("N: %d, MFLOP/s: %f\n",(int)(pow(1.5,N)),flops);
  

  free(a);
  free(b);
  free(d);
  return 0;
}