private static void RunSingleTest(ITest test)
{
State state = test.GenerateInitialState();
IMethod method = new Runge_Kutta(state.n);
//IMethod method = new MethodEuler();
const double dt = 3600 * 60;
double time;
for (time = 0; time < test.SuggestedFinalTime; time += dt)
{
bool resulttest = test.Compare(state, time);
if (resulttest == false)
{
//Console.WriteLine("Test {0} failed ", test.GetType().ToString());
// Console.WriteLine(time);
break;
}
Console.WriteLine(time);
state = method.Calculate(state, dt);
}
}
private static void DoCalculations()
{
System.Threading.Thread.CurrentThread.CurrentCulture = new System.Globalization.CultureInfo("en-US");
//State state = State.LoadFromFile("d:/1.txt");
//State state = (new AcceleratedMotionTest()).GenerateInitialState();
//State state = (new StraightMotionTest()).GenerateInitialState();
State state = (new EarthMotionTest()).GenerateInitialState();
// IView view = new TextStreamView(Console.Out);
string fileName = "d:/2.txt";
IView view = new GnuPlotView(@"D:/gnuplot/bin/gnuplot.exe");
using (StreamWriter writer = new StreamWriter(fileName, false))
{
//IView view = new TextStreamView(writer);
double t = 0;
const double tFinal = 2 * 31536000;
const double dt = 60;
view.Show(state, t);
IMethod method = new Runge_Kutta(state.n);
//IMethod method = new MethodEuler();
DateTime before = DateTime.Now;
for (t = 0; t < tFinal;)
{
state = method.Calculate(state, dt);
t += dt;
//вывести раз в секунду
DateTime current = DateTime.Now;
if ((current - before).TotalMilliseconds > 200) //relative to real time (sec)
{
view.Show(state, t);
before = current;
}
}
}
}