2018-12-29 15:24:35 +01:00
|
|
|
#!/bin/bash
|
|
|
|
file=$1
|
|
|
|
|
2019-01-14 15:20:51 +01:00
|
|
|
if [[ -z ${QP_ROOT} ]] ; then
|
2019-01-14 19:52:41 +01:00
|
|
|
>&2 echo "please source quantum_package.rc"
|
2019-01-14 15:20:51 +01:00
|
|
|
exit 1
|
|
|
|
fi
|
|
|
|
source ${QP_ROOT}/quantum_package.rc
|
2018-12-29 15:24:35 +01:00
|
|
|
|
2019-01-14 15:20:51 +01:00
|
|
|
qp_run print_e_conv $1
|
2019-01-04 02:17:55 +01:00
|
|
|
nstates=`cat ${1}/determinants/n_states`
|
|
|
|
echo $nstates
|
2018-12-29 15:24:35 +01:00
|
|
|
|
|
|
|
|
2019-01-14 15:20:51 +01:00
|
|
|
for i in $(seq 1 $nstates) ; do
|
|
|
|
out=${1}.${i}.conv
|
2019-01-04 02:17:55 +01:00
|
|
|
cat << EOF > ${out}.plt
|
2019-01-04 02:18:54 +01:00
|
|
|
set term pdf
|
|
|
|
set output "$out.pdf"
|
2019-01-04 02:17:55 +01:00
|
|
|
set log x
|
|
|
|
set xlabel "Number of determinants"
|
|
|
|
set ylabel "Total Energy (a.u.)"
|
2018-12-29 15:24:35 +01:00
|
|
|
|
2019-01-04 02:17:55 +01:00
|
|
|
plot "$out" w lp title "E_{var} state $i", "$out" u 1:3 w lp title "E_{var} + PT2 state $i"
|
2018-12-29 15:24:35 +01:00
|
|
|
|
2019-01-04 02:17:55 +01:00
|
|
|
EOF
|
2018-12-29 15:24:35 +01:00
|
|
|
|
2019-01-04 02:17:55 +01:00
|
|
|
gnuplot ${out}.plt
|
|
|
|
rm ${out}.plt
|
2018-12-29 15:24:35 +01:00
|
|
|
|
|
|
|
done
|
|
|
|
|
2019-01-14 15:20:51 +01:00
|
|
|
for i in $(seq 2 $nstates) ; do
|
|
|
|
out=${1}.${i}.delta_e.conv
|
2018-12-29 15:24:35 +01:00
|
|
|
cat << EOF > ${out}.plt
|
2019-01-04 02:18:54 +01:00
|
|
|
set term pdf
|
|
|
|
set output "$out.pdf"
|
2018-12-29 15:24:35 +01:00
|
|
|
set log x
|
|
|
|
set xlabel "Number of determinants"
|
2019-01-04 02:17:55 +01:00
|
|
|
set ylabel "Energy difference (a.u.)"
|
2018-12-29 15:24:35 +01:00
|
|
|
|
2019-01-04 02:17:55 +01:00
|
|
|
plot "$out" w lp title "Delta E_{var} state $i", "$out" u 1:3 w lp title "Delta E_{var} + PT2 state $i"
|
2018-12-29 15:24:35 +01:00
|
|
|
|
|
|
|
EOF
|
2019-01-14 15:20:51 +01:00
|
|
|
gnuplot ${out}.plt
|
|
|
|
rm ${out}.plt
|
2019-01-04 02:17:55 +01:00
|
|
|
done
|