# include "grid/cartesian1D.h"
#include "utils.h"
    
double uadv = 1.; 

void flux_upwind (scalar f, scalar df, double dt) 
{
    foreach()
        df[] = f[]*uadv*dt;
}

int main() {
    int n = 100;
    double CFL = 0.1;
    double t = 0., tend = 1.; 
    double dt = CFL/n/uadv;

    init_grid (n);

    periodic (right); 

    scalar f[], df[];
    foreach() {
        double k = 2.*pi;
        f[] = cos(k*x);
    }   


    FILE * fp    = fopen ("function.dat", "w");

    for (int i = 0; t <= tend; i++) {
    
      // we obtain the flux
      flux_upwind (f, df, dt);

      // numerical solution at t+dt
      foreach()
        f[] -= (df[]-df[-1])/Delta;
 
      if (i % 100 == 0) {
        foreach ()  
          fprintf(fp, "%g %g %g \n", t, x, f[]);
        fprintf(fp, " \n");
      }   

      t += dt; 
    }   


    fclose(fp);

    free_grid();
}