/// <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(); } }