srDFT_SC/Data/H10/plot_data.sh

54 lines
2.3 KiB
Bash
Raw Normal View History

2019-06-03 12:44:05 +02:00
#!/bin/bash
2019-06-04 13:59:23 +02:00
basis=tz
2019-06-03 12:44:05 +02:00
atom=H
echo data_${atom}10_v${basis}_cas1010
# DZ
FCI=-4.992784034175441
FCI_LDA=-4.992784034175441
FCI_PBE=-4.992784034175441
FCI_PBEont=-4.992784034175441
FCI_PBEontden=-4.992784034175441
2019-10-01 19:13:31 +02:00
## TZ
2019-06-03 12:44:05 +02:00
FCI=-4.998098113018433
FCI_LDA=-4.998098113018433
FCI_PBE=-4.998098113018433
FCI_PBEont=-4.998098113018433
FCI_PBEontden=-4.998098113018433
exact=-5.00000
cat <<EOF > plot.plt
2019-06-04 13:59:23 +02:00
plot "data_${atom}10_v${basis}_cas1010" u 1:2 w l lw 3 title "FCI ${basis}", $FCI w l lw 3 title "2*${atom} FCI", "" u 1:(\$2+\$6) w l lw 3 title "FCI ${basis}+LDA", $FCI_LDA w l lw 3 title "2*(FCI ${atom} + LDA)", 'exact-${atom}10' u (\$1 ):2 w l lw 3 smooth csplines title "Exact spline"
2019-06-03 12:44:05 +02:00
set term eps
set output "${atom}2_${basis}_LDA.eps"
replot
2019-06-04 13:59:23 +02:00
plot "data_${atom}10_v${basis}_cas1010" u 1:2 w l lw 3 title "FCI ${basis}", $FCI w l lw 3 title "2*${atom} FCI", "" u 1:(\$2+\$7) w l lw 3 title "FCI ${basis}+PBE", $FCI_PBE w l lw 3 title "2*(FCI ${atom} + PBE)", 'exact-${atom}10' u (\$1 ):2 w l lw 3 smooth csplines title "Exact spline"
2019-06-03 12:44:05 +02:00
set term eps
set output "${atom}2_${basis}_PBE.eps"
replot
2019-06-04 13:59:23 +02:00
plot "data_${atom}10_v${basis}_cas1010" u 1:2 w l lw 3 title "FCI ${basis}", $FCI w l lw 3 title "2*${atom} FCI", "" u 1:(\$2+\$8) w l lw 3 title "FCI ${basis}+PBE on top", $FCI_PBEont w l lw 3 title "2*(FCI ${atom} + PBE on top)", 'exact-${atom}10' u (\$1 ):2 w l lw 3 smooth csplines title "Exact spline"
2019-06-03 12:44:05 +02:00
set term eps
set output "${atom}2_${basis}_PBE_on_top.eps"
replot
2019-06-04 13:59:23 +02:00
plot "data_${atom}10_v${basis}_cas1010" u 1:2 w l lw 3 title "FCI ${basis}", $FCI w l lw 3 title "2*${atom} FCI", "" u 1:(\$2+\$9) w l lw 3 title "FCI ${basis}+PBE on top / density", $FCI_PBEontden w l lw 3 title "2*(FCI ${atom} + PBE on top / density)", 'exact-${atom}10' u (\$1 ):2 w l lw 3 smooth csplines title "Exact spline"
2019-06-03 12:44:05 +02:00
set term eps
set output "${atom}2_${basis}_PBE_on_top_den.eps"
replot
set yrange [:0.05]
2019-06-04 13:59:23 +02:00
plot "data_${atom}10_v${basis}_cas1010" u 1:(\$2 - $FCI) w l lw 3 title "FCI ${basis}", "" u 1:(\$2+\$6 -$FCI_PBEontden ) w l lw 3 title "FCI ${basis}+PBE on top / density" , 'exact-${atom}10' u (\$1 ):(\$2 - $exact) w l lw 3 smooth csplines title "Exact spline", 'exact-${atom}10' u (\$1 ):(\$2 - $exact) w p ps 0.5 notitle
2019-06-03 12:44:05 +02:00
set term eps
set output "${atom}2_${basis}_PBE_on_top_den_vs_exact.eps"
replot
EOF
gnuplot plot.plt