srDFT_SC/new/H10_vqz/plots/plot.sh

63 lines
1.8 KiB
Bash
Raw Normal View History

2019-10-02 18:43:27 +02:00
WF=FCI
BASIS=vqz
METHOD=DFT
TYPE=relat
FILE=data_${METHOD}_${BASIS}E_${TYPE}
OUT=${METHOD}_${BASIS}E_${TYPE}
#lt -1
cat << EOF > pouet.gp
2019-10-11 18:15:50 +02:00
set yrange [:0]
set xrange [:3.5]
set key bottom right
plot '${FILE}' using 1:2 smooth cspline notitle lt 2 , "" using 1:2 w p lt 2 ps 0.5 title "(MRCI+Q)/$BASIS"
replot '${FILE}' using 1:3 smooth cspline notitle lt 3 , "" using 1:3 w p lt 3 ps 0.5 title "(MRCI+Q)+PBE/$BASIS"
2019-10-02 18:43:27 +02:00
replot '${FILE}' using 1:5 smooth cspline notitle lt 7 , "" using 1:5 w p lt 7 ps 0.5 title "Exact"
set terminal eps linewidth 3
set output "${OUT}.eps"
replot
EOF
gnuplot pouet.gp
if [[ $METHOD == "DFT" ]]; then
OUT=${METHOD}_${BASIS}E_${TYPE}_zoom
cat << EOF > pouet.gp
2019-10-11 18:15:50 +02:00
set xrange [1.5:2.4]
set yrange [:-0.5]
set key bottom right
#set yrange [-0.38:-0.2]
plot '${FILE}' using 1:2 smooth cspline notitle lt 2 , "" using 1:2 w p lt 2 ps 0.5 title "(MRCI+Q)/$BASIS"
replot '${FILE}' using 1:3 smooth cspline notitle lt 3 , "" using 1:3 w p lt 3 ps 0.5 title "(MRCI+Q)+PBE/$BASIS"
replot '${FILE}' using 1:4 smooth cspline notitle lt 4 , "" using 1:4 w p lt 4 ps 0.5 title "(MRCI+Q)+PBEot/$BASIS"
2019-10-02 18:43:27 +02:00
replot '${FILE}' using 1:5 smooth cspline notitle lt 7 , "" using 1:5 w p lt 7 ps 0.5 title "Exact"
set terminal eps linewidth 3
set output "${OUT}.eps"
replot
EOF
fi
gnuplot pouet.gp
TYPE=error
FILE=data_${METHOD}_${BASIS}E_${TYPE}
OUT=${METHOD}_${BASIS}E_${TYPE}
#lt -1
cat << EOF > pouet.gp
plot '${FILE}' using 1:2 smooth cspline notitle lt 2 , "" using 1:2 w p lt 2 title "${WF}/$BASIS"
replot '${FILE}' using 1:3 smooth cspline notitle lt 3 , "" using 1:3 w p lt 3 title "${WF}+PBE/$BASIS"
replot '${FILE}' using 1:4 smooth cspline notitle lt 4 , "" using 1:4 w p lt 4 title "${WF}+PBEot/$BASIS"
set terminal eps linewidth 3
set output "${OUT}.eps"
replot
EOF
gnuplot pouet.gp