diff options
Diffstat (limited to 'plot.gnu')
| -rw-r--r-- | plot.gnu | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/plot.gnu b/plot.gnu new file mode 100644 index 0000000..2c7765a --- /dev/null +++ b/plot.gnu @@ -0,0 +1,27 @@ +set term qt persist + +exp_keff=real(system("cat results/explicit.out | grep 'Combined k-effective' | cut -d '=' -f 2 | cut -d '+' -f 1 | sed 's/ //g'")) +exp_pm=real(system("cat results/explicit.out | grep 'Combined k-effective' | cut -d '=' -f 2 | cut -d '-' -f 2 | sed 's/ //g'")) +exp_ppm=exp_keff+exp_pm +exp_mpm=exp_keff-exp_pm + +vwh_keff=real(system("cat results/vwh.out | grep 'Combined k-effective' | cut -d '=' -f 2 | cut -d '+' -f 1 | sed 's/ //g'")) +vwh_pm=real(system("cat results/vwh.out | grep 'Combined k-effective' | cut -d '=' -f 2 | cut -d '-' -f 2 | sed 's/ //g'")) +vwh_ppm=vwh_keff+vwh_pm +vwh_mpm=vwh_keff-vwh_pm + +set arrow from graph 0, first exp_keff to graph 1, first exp_keff nohead front lc rgb "red" dt 2 +set arrow from graph 0, first vwh_keff to graph 1, first vwh_keff nohead front lc rgb "blue" dt 2 + +set object 1 rectangle from graph 0, first exp_mpm to graph 1, first exp_ppm back fc rgb "red" fillstyle pattern 4 noborder +set object 2 rectangle from graph 0, first vwh_mpm to graph 1, first vwh_ppm back fc rgb "blue" fillstyle pattern 4 noborder + +set yrange [1.45:1.7] +set format y "%.5f" +set style fill pattern 4 noborder + +plot "results/rpt_table" using 1:2:3 with yerrorlines dt 2 title "RPT K_{eff}", \ + "results/rrpt_table" using 1:2:3 with yerrorlines dt 2 title "RRPT K_{eff}", \ + keyentry with lines title "Explicit K_{eff}" lc rgb "red" dt 2, \ + keyentry with lines title "VWH K_{eff}" lc rgb "blue" dt 2 + |