/// <summary>
/// Utility/Debug function, plots numbers at specific coordinates.
/// </summary>
/// <param name="X"></param>
/// <param name="filename"></param>
public static void PlotCoordinateLabels(this MultidimensionalArray X, string filename)
{
using (var gp = new Gnuplot()) {
int L = X.NoOfRows;
if (X.NoOfCols != 2)
{
throw new NotSupportedException("works only for 2D grid");
}
gp.Cmd("set terminal png");
gp.Cmd("set output \"{0}\"", filename);
gp.SetXRange(-4, 4);
gp.SetYRange(-4, 4);
for (int l = 0; l < L; l++)
{
gp.Cmd("set label \"{0}\" at {1},{2}", l, X[l, 0], X[l, 1]);
}
gp.PlotSlope(0, 0);
gp.Execute();
}
}
/// <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();
}
}
