#include <fstream>

#include <iomanip>

#include <math.h>

#include <iostream>

#include "parabolic_solver.h"

#include "algorithm.h"





const double PI = atan(1.0) * 4.0;

const double a = 1.0;

const double b = 0.0;

const double c = 0.0;



double phi1(double x, double t)

{

	  return exp(-1*a*t);

	//return sin(t);

	//return (  -exp(-a * a * t) * (cos(b * t) + sin (b*t))  );

	//return (  exp((c - a*a)*t) * (cos (b * t) + sin (b * t))  );

}



double phi2(double x, double t)

{

	    return -1*exp(-1*a*t);

	//return -sin(t);

	//return (  exp(-a * a * t) * (cos(b * t) + sin (b*t))  );

	//return (  -exp((c - a*a)*t) * (cos (b * t) + sin (b * t))  );

}



double psi(double x, double t)

{

	    return  sin(x);

	//return 0.0;

	//return cos (x);

	//return sin (x);

}



double f(double x, double t)

{

	//return cos (x) * (cos (t)  + sin(t));

	return 0.0;

}







/*

void main()

{

	CParabolicSolver direct_solve;

	ReverseSolve reverse_solve;



/ *

	direct_solve.setGrid(10, 10, 0.0, PI, 0.3);

	direct_solve.setTask(a, b, c, f, 1.0, 0.0, 1.0, 0.0, phi1, phi2, psi);

	direct_solve.solve_implicit();

	direct_solve.CheckSolve();* /





	int num = direct_solve.GetK();

	reverse_solve.SetTask(a, b, c, 1.0, 0.0, 1.0, 0.0, PI, 0.003, 0.0, 2, 10);

	//reverse_solve.SetLayer(direct_solve.GetLayer(num));

	reverse_solve.GetSolve();

	reverse_solve.CheckSolve("rev_sol.txt");

	//reverse_solve.CheckResult(direct_solve.GetTable(),"test.txt");

}*/