#include "grid/cartesian1D.h"
#include "run.h"

scalar T[], dT[];
double dt;
double alpha = 1;

int main()
{
    // periodic(left);
    L0 = 4 * pi;
    X0 = -L0 / 2.;
    N = 1 << 9;
    // stable condition of Lax - Friedrichs method dt < dx/c
    DT = L0 / N / alpha;
    run();
}

T[left] = dirichlet(0.);
T[right] = dirichlet(0.);

event init(t = 0)
{
    foreach ()
    {
        T[] = 0.5 * (fabs(x) < 1);
    }
    boundary({T});
}

// event printdata (t = 0; t <= 1000 * DT; t += 100 * DT) {
event printdata(t = 0; t <= 1; t += 0.1)
{
    // event printdata (t = 0.5) {
    static FILE *fp = fopen("AD1D", "w");
    foreach ()
        fprintf(fp, "%g %g %g %g\n", x, T[], dT[], t);
    fprintf(fp, "\n\n");
    fflush(fp);
}

event integration(i++)
{
    double dt = DT;
    dt = dtnext(dt);
    foreach ()
        // Lax - Friedrich method
        dT[] = -alpha * (T[1] - T[-1]) / (2 * Delta);
    foreach ()
        T[] = 0.5 * (T[1] + T[-1]) + dt * dT[];
    boundary({T});
}