/// <summary>
/// Simple plotting interface
/// </summary>
/// <returns>Output of <see cref="GnuplotExtensions.PlotNow(Gnuplot)"/></returns>
static public object Plot(IEnumerable <double> X1, IEnumerable <double> Y1, string Name1 = null, string Format1 = null,
IEnumerable <double> X2 = null, IEnumerable <double> Y2 = null, string Name2 = null, string Format2 = null,
IEnumerable <double> X3 = null, IEnumerable <double> Y3 = null, string Name3 = null, string Format3 = null,
IEnumerable <double> X4 = null, IEnumerable <double> Y4 = null, string Name4 = null, string Format4 = null,
IEnumerable <double> X5 = null, IEnumerable <double> Y5 = null, string Name5 = null, string Format5 = null,
IEnumerable <double> X6 = null, IEnumerable <double> Y6 = null, string Name6 = null, string Format6 = null,
IEnumerable <double> X7 = null, IEnumerable <double> Y7 = null, string Name7 = null, string Format7 = null,
bool logX = false, bool logY = false)
{
using (var gp = new Gnuplot()) {
IEnumerable <double>[] Xs = new[] { X1, X2, X3, X4, X5, X6, X7 };
IEnumerable <double>[] Ys = new[] { Y1, Y2, Y3, Y4, Y5, Y6, Y7 };
string[] Ns = new string[] { Name1, Name2, Name3, Name4, Name5, Name6, Name7 };
string[] Fs = new string[] { Format1, Format2, Format3, Format4, Format5, Format6, Format7 };
for (int i = 0; i < 7; i++)
{
if (Ys[i] != null)
{
var f1 = new PlotFormat();
if (Fs[i] != null)
{
f1.FromString(Fs[i]);
}
gp.PlotXY(Xs[i], Ys[i], title: Ns[i], format: f1, logX: logX, logY: logY);
}
}
return(gp.PlotNow());
}
}
/// <summary>
/// Visualizes the stored data using gnuplot.
/// </summary>
public void Plot()
{
Gnuplot gp = new Gnuplot();
gp.SetXLabel("x");
gp.SetYLabel("y");
gp.Cmd("set terminal wxt noraise");
if (LogX)
{
gp.Cmd("set logscale x");
gp.Cmd("set format x \"10^{%L}\"");
}
if (LogY)
{
gp.Cmd("set logscale y");
gp.Cmd("set format y \"10^{%L}\"");
}
gp.Cmd("set grid xtics ytics");
int lineColor = 0;
foreach (var group in dataGroups)
{
gp.PlotXY(group.Abscissas, group.Values, group.Name,
new PlotFormat(lineColor: ((LineColors)(++lineColor)), pointType: ((PointTypes)4), pointSize: 1.5, Style: Styles.LinesPoints));
}
gp.Execute();
}
/// <summary>
/// Plots a line plot from a given residual log.
/// The y-axis will be scaled logarithmically.
/// </summary>
private Image drawGnuplotGraph()
{
Image graph;
using (Gnuplot gp = new Gnuplot()) {
gp.SetXLabel("Iteration");
gp.SetYLabel(log.Norm + " Norm");
//gp.Cmd("set terminal wxt noraise title 'SessionGuid: " + log.SessionGuid.ToString() + "'");
gp.Cmd("set title " + titleOfPlot);
gp.Cmd("set key outside under horizontal box");
gp.Cmd("set logscale y");
gp.Cmd("set format y \"10^{%L}\"");
gp.Cmd("set grid xtics ytics");
List <string> KeysToPlot = log.Values.Keys.ToList();
KeysToPlot.Remove("#Line");
KeysToPlot.Remove("#Time");
KeysToPlot.Remove("#Iter");
for (int i = 0; i < KeysToPlot.Count; i++)
{
gp.PlotXY(log.Values["#Time"], log.Values[KeysToPlot[i]], KeysToPlot[i],
new PlotFormat(lineColor: ((LineColors)i), Style: Styles.Lines));
}
graph = gp.PlotGIF();
}
return(graph);
}
/// <summary>
/// Draws a graph with Gnuplot and saves it as a gif. Subsequently sets Picture of Form to the plot.
/// </summary>
private void drawGnuplotGraph()
{
using (Gnuplot gp = new Gnuplot()) {
gp.SetXLabel("Iteration");
gp.SetYLabel(log.Norm + " Norm");
gp.Cmd("set title \"Residual Plot\"");
gp.Cmd("set key outside under horizontal box");
gp.Cmd("set logscale y");
gp.Cmd("set format y \"10^{%L}\"");
gp.Cmd("set grid xtics ytics");
List <string> KeysToPlot = log.Values.Keys.ToList();
KeysToPlot.Remove("#Line");
KeysToPlot.Remove("#Time");
KeysToPlot.Remove("#Iter");
for (int i = 0; i < KeysToPlot.Count; i++)
{
gp.PlotXY(log.Values["#Time"], log.Values[KeysToPlot[i]], KeysToPlot[i],
new PlotFormat(lineColor: ((LineColors)i), Style: Styles.Lines));
}
Image graph = gp.PlotGIF();
pictureBox1.Image = graph;
pictureBox1.Refresh();
}
}
/// <summary>
/// plot a 2D grid
/// </summary>
public static void Plot2DGrid(this GridCommons grd, params int[] Cells)
{
using (var gp = new Gnuplot()) {
Console.WriteLine("Plotting 2D grid with gnuplot...");
if (grd.SpatialDimension != 2)
{
throw new NotSupportedException("works only for 2D grid");
}
int J = grd.Cells.Length;
if (Cells == null)
{
Cells = J.ForLoop(j => j);
}
foreach (int j in Cells)
{
var Cell_j = grd.Cells[j];
//var Kref = grd.RefElements.Single(KK => KK.SupportedCellTypes.Contains(Cell_j.Type));
//var Vtx = Kref.Vertices;
//var _Vtx = Kref.GetSubDivTree(3).GlobalVertice;
//var Vtx = MultidimensionalArray.Create(_Vtx.GetLength(0), _Vtx.GetLength(1));
//Vtx.SetA2d(_Vtx);
//Vtx.Scale(1);
var Vtx_glob = Cell_j.TransformationParams;
double[] xNodes = Vtx_glob.ExtractSubArrayShallow(-1, 0).To1DArray();
double[] yNodes = Vtx_glob.ExtractSubArrayShallow(-1, 1).To1DArray();
double xC = xNodes.Sum() / xNodes.Length;
double yC = yNodes.Sum() / yNodes.Length;
for (int k = 0; k < xNodes.Length; k++)
{
double dx = xNodes[k] - xC;
dx *= 0.95;
xNodes[k] = xC + dx;
double dy = yNodes[k] - yC;
dy *= 0.95;
yNodes[k] = yC + dy;
}
double hy = yNodes.Max() - yNodes.Min();
gp.PlotXY(xNodes, yNodes, title: null, format: (new PlotFormat(Style: Styles.LinesPoints, lineColor: ((LineColors)j))));
gp.Cmd("set label \"{0}\" at {2},{3} font \"Arial,5\"", j.ToString(), Cell_j.GlobalID, xC.ToStringDot(), (yC + hy * 0.21).ToStringDot());
gp.Cmd("set label \"[{1}]\" at {2},{3} font \"Arial,5\"", j.ToString(), Cell_j.GlobalID, xC.ToStringDot(), (yC - hy * 0.21).ToStringDot());
}
gp.Execute();
Console.WriteLine("Press any key to continue...");
Console.ReadKey();
}
}
/// <summary>
/// Writes plot commands to a gnuplot object.
/// </summary>
public void ToGnuplot(Gnuplot gp)
{
//gp.Cmd($"set terminal enhanced font \"Helvetica, 30\" ");
// =======
// margins
// =======
{
if (lmargin != null)
{
gp.Cmd("set lmargin {0:0.####e-00}", this.lmargin.Value);
}
else
{
gp.Cmd("unset lmargin ");
}
if (rmargin != null)
{
gp.Cmd("set rmargin {0:0.####e-00}", this.rmargin.Value);
}
else
{
gp.Cmd("unset rmargin ");
}
if (tmargin != null)
{
gp.Cmd("set tmargin {0:0.####e-00}", this.tmargin.Value);
}
else
{
gp.Cmd("unset tmargin ");
}
if (bmargin != null)
{
gp.Cmd("set bmargin {0:0.####e-00}", this.bmargin.Value);
}
else
{
gp.Cmd("unset bmargin ");
}
}
// ============
// log settings
// ============
{
if (this.LogX)
{
gp.Cmd("set logscale x " + this.LogBaseX);
}
else
{
gp.Cmd("unset logscale x");
}
if (this.LogY)
{
gp.Cmd("set logscale y " + this.LogBaseY);
}
else
{
gp.Cmd("unset logscale y");
}
if (this.LogX2)
{
gp.Cmd("set logscale x2 " + this.LogBaseX);
}
else
{
gp.Cmd("unset logscale x2");
}
if (this.LogY2)
{
gp.Cmd("set logscale y2 " + this.LogBaseY);
}
else
{
gp.Cmd("unset logscale y2");
}
}
// ==========
// axis range
// ==========
{
if ((this.XrangeMax != null) != (this.XrangeMin != null))
{
throw new ArgumentException("X range minimum and maximum must be set either both or none.");
}
if ((this.YrangeMax != null) != (this.YrangeMin != null))
{
throw new ArgumentException("Y range minimum and maximum must be set either both or none.");
}
if (this.XrangeMin != null)
{
if (this.XrangeMin.Value >= this.XrangeMax.Value)
{
throw new ArgumentException("X range maximum must be grater than minimum.");
}
gp.SetXRange(this.XrangeMin.Value, this.XrangeMax.Value);
}
else
{
gp.SetXAutorange();
}
if (this.YrangeMin != null)
{
if (this.YrangeMin.Value >= this.YrangeMax.Value)
{
throw new ArgumentException("Y range maximum must be grater than minimum.");
}
gp.SetYRange(this.YrangeMin.Value, this.YrangeMax.Value);
}
else
{
gp.SetYAutorange();
}
if ((this.X2rangeMax != null) != (this.X2rangeMin != null))
{
throw new ArgumentException("X2 range minimum and maximum must be set either both or none.");
}
if ((this.Y2rangeMax != null) != (this.Y2rangeMin != null))
{
throw new ArgumentException("Y2 range minimum and maximum must be set either both or none.");
}
if (this.X2rangeMin != null)
{
if (this.X2rangeMin.Value >= this.X2rangeMax.Value)
{
throw new ArgumentException("X range maximum must be grater than minimum.");
}
gp.SetX2Range(this.X2rangeMin.Value, this.X2rangeMax.Value);
}
else
{
gp.SetX2Autorange();
}
if (this.Y2rangeMin != null)
{
if (this.Y2rangeMin.Value >= this.Y2rangeMax.Value)
{
throw new ArgumentException("Y2 range maximum must be grater than minimum.");
}
gp.SetY2Range(this.Y2rangeMin.Value, this.Y2rangeMax.Value);
}
else
{
gp.SetY2Autorange();
}
}
// ========================
// labels, title, legend...
// ========================
{
gp.SetXLabel(this.Xlabel);
gp.SetYLabel(this.Ylabel);
gp.SetX2Label(this.X2label);
gp.SetY2Label(this.Y2label);
gp.SetTitle(this.Title);
if (this.ShowLegend)
{
gp.Cmd("unset key");
string command = $"set key font \",{this.LegendFont}\"";
if ((this.LegendPosition != null) & (this.LegendAlignment != null))
{
System.Console.WriteLine("legend position and legend alignment is set. Choose only one of them! Ignoring alignment ...");
}
if (this.LegendPosition != null)
{
command += String.Format("at {1:0.####e-00},{2:0.####e-00} vertical maxrows {0} ", this.dataGroups.Length, this.LegendPosition[0], this.LegendPosition[1]);
}
else if (this.LegendAlignment != null)
{
Dictionary <string, string> alignments = new Dictionary <string, string>();
alignments.Add("r", "right");
alignments.Add("c", "center");
alignments.Add("l", "left");
alignments.Add("b", "bottom");
alignments.Add("t", "top");
switch (LegendAlignment[0])
{
case "o":
//alignment within plotboundary
command += "outside ";
break;
case "i":
//alignement within graphyboundary
command += "inside ";
break;
default:
throw new ArgumentException("this style is not specified: use [i/o]");
}
for (int i = 1; i <= 2; i++)
{
foreach (KeyValuePair <string, string> kvp in alignments)
{
if (LegendAlignment[i] == kvp.Key)
{
command += kvp.Value + " ";
}
}
}
}
else
{
//throw new ArgumentNullException("no alignment or position chosen");
}
if (this.LegendHorizontal == true)
{
command += "horizontal ";
}
if (this.LegendSwap == true)
{
command += "Left reverse ";
}
System.Console.WriteLine(command);
gp.Cmd(command);
}
else
{
gp.Cmd("set key off");
}
}
// ====
// tics
// ====
{
if (this.ShowXtics)
{
if (this.LogX)
{
gp.Cmd("set xtics format \"$" + this.LogBaseX + "^{%L}$\" ");
gp.Cmd("set xtics offset 0, 0-0.4 font \"sans, 18\" ");
gp.Cmd($"set xtics font \"sans, {this.LabelFont}\" ");
}
else
{
gp.Cmd("set xtics ");
}
}
else
{
gp.Cmd("set xtics format \" \" "); // shows ticks/marks, but hides numbers
}
if (this.ShowX2tics)
{
if (this.LogX2)
{
gp.Cmd("set x2tics format \"$" + this.LogBaseX + "^{%L}$\" ");
}
else
{
gp.Cmd("set x2tics ");
}
}
else
{
gp.Cmd("set x2tics format \" \" "); // shows ticks/marks, but hides numbers
}
if (this.ShowYtics)
{
if (this.LogY)
{
gp.Cmd("set ytics format \"$" + this.LogBaseY + "^{%L}$\" ");
gp.Cmd($"set ytics font \"sans, {this.LabelFont}\" ");
}
else
{
gp.Cmd("set ytics ");
}
}
else
{
gp.Cmd("set ytics format \" \" "); // shows ticks/marks, but hides numbers
}
if (this.ShowY2tics)
{
if (this.LogY2)
{
gp.Cmd("set y2tics format \"$" + this.LogBaseY + "^{%L}$\" ");
}
else
{
gp.Cmd("set y2tics ");
}
}
else
{
//gp.Cmd("unset y2tics");
gp.Cmd("set y2tics format \" \" "); // shows ticks/marks, but hides numbers
}
}
// =================
// finally, plotting
// =================
{
foreach (var xyData in this.dataGroups)
{
gp.PlotXY(xyData.Abscissas, xyData.Values, xyData.Name, xyData.Format, useX2: xyData.UseX2, useY2: xyData.UseY2, deferred: true);
}
gp.WriteDeferredPlotCommands();
}
}
static void Plot2dGridGnuplot(GridCommons grd)
{
using (var gp = new Gnuplot()) {
int J = grd.Cells.Length;
double xmin = double.MaxValue, xmax = double.MinValue, ymin = double.MaxValue, ymax = double.MinValue;
for (int j = 0; j < J; j++)
{
var Cell = grd.Cells[j];
var Kref = grd.GetRefElement(Cell.Type);
int[] R;
if (Kref.GetType() == typeof(Triangle))
{
R = new int[] { 0, 1, 2, 0 };
}
else if (Kref.GetType() == typeof(Square))
{
R = new int[] { 0, 1, 3, 2, 0 };
}
else
{
throw new NotSupportedException();
}
int I = 4;
int K = 20;
NodeSet LocNodes = new NodeSet(Kref, I * K * (R.Length - 1), 2);
var vtx = Kref.Vertices;
double alpha = 1.0 / (K - 1);
for (int iFace = 0; iFace < R.Length - 1; iFace++)
{
for (int k = 0; k < K; k++)
{
double a = alpha * k;
for (int d = 0; d < 2; d++)
{
LocNodes[iFace * K + k, d] = vtx[R[iFace], d] * (1 - a) + vtx[R[iFace + 1], d] * a;
}
}
}
for (int i = 0; i < I; i++)
{
int ind0 = K * (R.Length - 1) * i;
int indE = K * (R.Length - 1) * (i + 1) - 1;
int indp0 = K * (R.Length - 1) * (i - 1);
int indpE = K * (R.Length - 1) * i - 1;
var LocNodes_i = LocNodes.ExtractSubArrayShallow(new int[] { ind0, 0 }, new int[] { indE, 1 });
if (i > 0)
{
var LocNodes_iP = LocNodes.ExtractSubArrayShallow(new int[] { indp0, 0 }, new int[] { indpE, 1 });
LocNodes_i.Set(LocNodes_iP);
LocNodes_i.Scale(0.65);
}
else
{
LocNodes_i.Scale(0.9);
}
}
LocNodes.LockForever();
MultidimensionalArray GlobalNodes = Transform(Kref, Cell, LocNodes);
xmin = Math.Min(xmin, GlobalNodes.ExtractSubArrayShallow(-1, 0).Min());
xmax = Math.Max(xmax, GlobalNodes.ExtractSubArrayShallow(-1, 0).Max());
ymin = Math.Min(ymin, GlobalNodes.ExtractSubArrayShallow(-1, 1).Min());
ymax = Math.Max(ymax, GlobalNodes.ExtractSubArrayShallow(-1, 1).Max());
gp.PlotXY(GlobalNodes.GetColumn(0), GlobalNodes.GetColumn(1), title: ("tri" + j),
format: new PlotFormat(lineColor: ((LineColors)j)));
}
gp.SetXRange(xmin - 0.1, xmax + 0.1);
gp.SetYRange(ymin - 0.1, ymax + 0.1);
gp.Execute();
Console.WriteLine("press any key to continue...");
Console.ReadKey();
}
}
/// <summary>
/// Writes plot commands to a gnuplot object.
/// </summary>
public void ToGnuplot(Gnuplot gp)
{
// =======
// margins
// =======
{
if (lmargin != null)
{
gp.Cmd("set lmargin {0:0.####e-00}", this.lmargin.Value);
}
else
{
gp.Cmd("unset lmargin ");
}
if (rmargin != null)
{
gp.Cmd("set rmargin {0:0.####e-00}", this.rmargin.Value);
}
else
{
gp.Cmd("unset rmargin ");
}
if (tmargin != null)
{
gp.Cmd("set tmargin {0:0.####e-00}", this.tmargin.Value);
}
else
{
gp.Cmd("unset tmargin ");
}
if (bmargin != null)
{
gp.Cmd("set bmargin {0:0.####e-00}", this.bmargin.Value);
}
else
{
gp.Cmd("unset bmargin ");
}
}
// ============
// log settings
// ============
{
if (this.LogX)
{
gp.Cmd("set logscale x");
}
else
{
gp.Cmd("unset logscale x");
}
if (this.LogY)
{
gp.Cmd("set logscale y");
}
else
{
gp.Cmd("unset logscale y");
}
if (this.LogX2)
{
gp.Cmd("set logscale x2");
}
else
{
gp.Cmd("unset logscale x2");
}
if (this.LogY2)
{
gp.Cmd("set logscale y2");
}
else
{
gp.Cmd("unset logscale y2");
}
}
// ==========
// axis range
// ==========
{
if ((this.XrangeMax != null) != (this.XrangeMin != null))
{
throw new ArgumentException("X range minimum and maximum must be set either both or none.");
}
if ((this.YrangeMax != null) != (this.YrangeMin != null))
{
throw new ArgumentException("Y range minimum and maximum must be set either both or none.");
}
if (this.XrangeMin != null)
{
if (this.XrangeMin.Value >= this.XrangeMax.Value)
{
throw new ArgumentException("X range maximum must be grater than minimum.");
}
gp.SetXRange(this.XrangeMin.Value, this.XrangeMax.Value);
}
else
{
gp.SetXAutorange();
}
if (this.YrangeMin != null)
{
if (this.YrangeMin.Value >= this.YrangeMax.Value)
{
throw new ArgumentException("Y range maximum must be grater than minimum.");
}
gp.SetYRange(this.YrangeMin.Value, this.YrangeMax.Value);
}
else
{
gp.SetYAutorange();
}
if ((this.X2rangeMax != null) != (this.X2rangeMin != null))
{
throw new ArgumentException("X2 range minimum and maximum must be set either both or none.");
}
if ((this.Y2rangeMax != null) != (this.Y2rangeMin != null))
{
throw new ArgumentException("Y2 range minimum and maximum must be set either both or none.");
}
if (this.X2rangeMin != null)
{
if (this.X2rangeMin.Value >= this.X2rangeMax.Value)
{
throw new ArgumentException("X range maximum must be grater than minimum.");
}
gp.SetX2Range(this.X2rangeMin.Value, this.X2rangeMax.Value);
}
else
{
gp.SetX2Autorange();
}
if (this.Y2rangeMin != null)
{
if (this.Y2rangeMin.Value >= this.Y2rangeMax.Value)
{
throw new ArgumentException("Y2 range maximum must be grater than minimum.");
}
gp.SetY2Range(this.Y2rangeMin.Value, this.Y2rangeMax.Value);
}
else
{
gp.SetY2Autorange();
}
}
// ========================
// labels, title, legend...
// ========================
{
gp.SetXLabel(this.Xlabel);
gp.SetYLabel(this.Ylabel);
gp.SetX2Label(this.X2label);
gp.SetY2Label(this.Y2label);
gp.SetTitle(this.Title);
if (this.ShowLegend)
{
gp.Cmd("unset key");
if (this.LegendPosition != null)
{
gp.Cmd("set key at {1:0.####e-00},{2:0.####e-00} vertical maxrows {0} ", this.dataGroups.Length, this.LegendPosition[0], this.LegendPosition[1]);
}
else
{
gp.Cmd("set key outside right vertical maxrows {0} ", this.dataGroups.Length);
}
}
else
{
gp.Cmd("set key off");
}
}
// ====
// tics
// ====
{
if (this.ShowXtics)
{
if (this.LogX)
{
gp.Cmd("set xtics format \"$10^{%T}$\" ");
}
else
{
gp.Cmd("set xtics ");
}
}
else
{
gp.Cmd("set xtics format \" \" "); // shows ticks/marks, but hides numbers
}
if (this.ShowX2tics)
{
if (this.LogX2)
{
gp.Cmd("set x2tics format \"$10^{%T}$\" ");
}
else
{
gp.Cmd("set x2tics ");
}
}
else
{
gp.Cmd("set x2tics format \" \" "); // shows ticks/marks, but hides numbers
}
if (this.ShowYtics)
{
if (this.LogY)
{
gp.Cmd("set ytics format \"$10^{%T}$\" ");
}
else
{
gp.Cmd("set ytics ");
}
}
else
{
gp.Cmd("set ytics format \" \" "); // shows ticks/marks, but hides numbers
}
if (this.ShowY2tics)
{
if (this.LogY2)
{
gp.Cmd("set y2tics format \"$10^{%T}$\" ");
}
else
{
gp.Cmd("set y2tics ");
}
}
else
{
//gp.Cmd("unset y2tics");
gp.Cmd("set y2tics format \" \" "); // shows ticks/marks, but hides numbers
}
}
// =================
// finally, plotting
// =================
{
foreach (var xyData in this.dataGroups)
{
gp.PlotXY(xyData.Abscissas, xyData.Values, xyData.Name, xyData.Format, useX2: xyData.UseX2, useY2: xyData.UseY2, deferred: true);
}
gp.WriteDeferredPlotCommands();
}
}