public void DrawHist(TimeList<Complex[]> list, double t)
{
using (Context context = new Context(Surface)) {
context.Rectangle(0, 0, Width, Height);
context.Color = Background;
context.Fill();
int index = list.Index (t);
int times = 30;
double xl = (double)(Width - 60) / N;
double yl = (double)Height / times;
Complex[] points;
Color c;
int ii;
for (int i = index; i >= Math.Max(0, index - times); i--) {
points = list.Get (i);
N = points.Length;
ii = index - i;
for (int j = 0; j < points.Length; j++) {
DrawHistArea (context, xl, yl, ii, j, _valueExtractor(points [j]));
}
}
}
DrawExpose (null, null);
}
public void CreateAll(string path, TimeList<Complex[]> list, double t)
{
if(!Directory.Exists(path)) {
Directory.CreateDirectory (path);
//throw new ArgumentException("Directory already exists");
}
Complex[] points = list[t];
int index = list.Index (t);
int N = list.Get (index).Length;
string name = Path.GetFileName (path);
using(StreamWriter w = new StreamWriter(Path.Combine(path, name + "_single.dat"))) {
for (int i = 0; i < N; i++) {
w.WriteLine ("{0} {1} {2} {3} {4}", i, points [i].Real.F(), points [i].Imaginary.F(), points [i].Magnitude.F(), points [i].Phase.F());
}
w.Flush ();
w.Close ();
}
double times = 10;
int firstTime;
if(t - times <= 0) {
firstTime = 0;
} else {
firstTime = list.Index (t - times);
}
double[] cmX = new double[N];
double[] cmY = new double[N];
using(StreamWriter w = new StreamWriter(Path.Combine(path, name + "_hist.dat"))) {
for (int ti = firstTime; ti <= index; ti++) {
points = list.Get (ti);
w.Write (points[0].Real.F());
for (int j = 1; j < N; j++) {
w.Write(" ");
w.Write (points[j].Real.F());
cmX [j] += points [j].Real;
}
w.WriteLine ();
w.Write (points[0].Imaginary.F());
for (int j = 1; j < N; j++) {
w.Write(" ");
w.Write (points[j].Imaginary.F());
cmY [j] += points [j].Imaginary;
}
w.WriteLine ();
w.Write (points[0].Magnitude.F());
for (int j = 1; j < N; j++) {
w.Write(" ");
w.Write (points[j].Magnitude.F());
}
w.WriteLine ();
}
w.Flush ();
w.Close ();
}
for (int i = 0; i < N; i++) {
cmX [i] = cmX [i] / (index - firstTime);
cmY [i] = cmY [i] / (index - firstTime);
}
using(StreamWriter w = new StreamWriter(Path.Combine(path, name + "_cm.dat"))) {
for (int j = 0; j < N; j++) {
w.WriteLine ("{0} {1}", j, Math.Sqrt (cmX [j] * cmX [j] + cmY [j] * cmY [j]).F());
}
w.Flush ();
w.Close ();
}
using(StreamWriter w = new StreamWriter(Path.Combine(path, name + ".plt"))) {
w.WriteLine ("set terminal png");
w.WriteLine ("set termopt enhanced");
w.WriteLine ("set autoscale");
w.WriteLine ("set output \"x.png\"");
w.WriteLine ("set title \"Plot 1\"");
w.WriteLine ("set xlabel \"j\"");
w.WriteLine ("set ylabel \"x_j\"");
w.WriteLine ("plot \"" + name + "_single.dat\" using 1:2 notitle");
w.WriteLine ("set output \"y.png\"");
w.WriteLine ("set title \"Plot 2\"");
w.WriteLine ("set xlabel \"j\"");
w.WriteLine ("set ylabel \"y_j\"");
w.WriteLine ("plot \"" + name + "_single.dat\" using 1:3 notitle");
w.WriteLine ("set output \"xy.png\"");
w.WriteLine ("set title \"Plot 3\"");
w.WriteLine ("set xlabel \"j\"");
w.WriteLine ("set ylabel \"y_j, x_j\"");
w.WriteLine ("plot \"" + name + "_single.dat\" using 1:2 title \"x\", \"" + name + "_single.dat\" using 1:3 title \"y\"");
w.WriteLine ("set output \"r.png\"");
w.WriteLine ("set title \"Plot 3\"");
w.WriteLine ("set xlabel \"j\"");
w.WriteLine ("set ylabel \"r_j\"");
w.WriteLine ("plot \"" + name + "_single.dat\" using 1:4 notitle");
w.WriteLine ("set output \"cm.png\"");
w.WriteLine ("set title \"Massenschwerpunkt\"");
w.WriteLine ("set xlabel \"j\"");
w.WriteLine ("set ylabel \"r_{c.m.}\"");
w.WriteLine ("plot \"" + name + "_cm.dat\" using 1:2 notitle");
w.WriteLine ("set output \"Phasenraum.png\"");
w.WriteLine ("set title \"Phasenraum\"");
w.WriteLine ("set xlabel \"Realteil\"");
w.WriteLine ("set ylabel \"Imaginärteil\"");
w.WriteLine ("plot \"" + name + "_single.dat\" using 2:3 notitle");
w.WriteLine ("set terminal png");
w.WriteLine ("set autoscale fix");
w.WriteLine ("set tics out nomirror");
w.WriteLine ("set pm3d map interpolate 0,0");
w.WriteLine ("set xlabel \"j\"");
w.WriteLine ("set ylabel \"t\"");
w.WriteLine ("set output \"hist_x.png\"");
w.WriteLine ("set title \"x\"");
w.WriteLine ("splot '" + name + "_hist.dat' every :3 matrix using 1:($2/" + ((t - list.TimeOf(firstTime)) * 30).F() + "):3 notitle");
w.WriteLine ("set output \"hist_y.png\"");
w.WriteLine ("set title \"y\"");
w.WriteLine ("splot '" + name + "_hist.dat' every :3::1 matrix using 1:(($2 - 1)/" + ((t - list.TimeOf(firstTime)) * 30).F() + "):3 notitle");
w.WriteLine ("set output \"hist_r.png\"");
w.WriteLine ("set title \"r\"");
w.WriteLine ("splot '" + name + "_hist.dat' every :3::2 matrix using 1:(($2 - 2)/" + ((t - list.TimeOf(firstTime)) * 30).F() + "):3 notitle");
w.Flush ();
w.Close ();
}
System.Diagnostics.Process proc = new System.Diagnostics.Process ();
proc.EnableRaisingEvents = false;
proc.StartInfo.CreateNoWindow = true;
proc.StartInfo.WorkingDirectory = path;
proc.StartInfo.FileName = "gnuplot";
proc.StartInfo.Arguments = Path.Combine(path, name + ".plt");
proc.StartInfo.UseShellExecute = false;
proc.StartInfo.RedirectStandardOutput = true;
proc.Start ();
proc.WaitForExit ();
string output = proc.StandardOutput.ReadToEnd ();
proc.Close ();
}
