blob: 30b6de9f2a9ab2db9000c9d402fd4f3aaffe782e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
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.06:1.1]
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
|
