cutoffdiff="0.01" # energy cutoff variation - more then this rises warning
elif [ "$1" == "prota_MIN_CART" ]; then
extremediff="10.0" # extreme energy difference, comething went terribly wrong
cutoffdiff="0.01" # energy cutoff variation - more then this rises warning
elif [ "$1" == "prota_MIN_CART" ]; then
extremediff="10.0" # extreme energy difference, comething went terribly wrong
cutoffdiff="0.1" # energy cutoff variation - more then this rises warning
sumsl_return=`grep SUMSL $file|awk '{print $4}'`
echo 'SUMSL return code' $sumsl_return
if [ "$sumsl_return" != "4" ]; then
echo 'ERROR = SUMSL return code' $sumsl_return 'is not 4'
cutoffdiff="0.1" # energy cutoff variation - more then this rises warning
sumsl_return=`grep SUMSL $file|awk '{print $4}'`
echo 'SUMSL return code' $sumsl_return
if [ "$sumsl_return" != "4" ]; then
echo 'ERROR = SUMSL return code' $sumsl_return 'is not 4'
fi
elif [ "$1" == "1l2y_MIN_INT" ]; then
extremediff="10.0" # extreme energy difference, comething went terribly wrong
fi
elif [ "$1" == "1l2y_MIN_INT" ]; then
extremediff="10.0" # extreme energy difference, comething went terribly wrong
cutoffdiff="0.1" # energy cutoff variation - more then this rises warning
sumsl_return=`grep SUMSL $file|awk '{print $4}'`
echo 'SUMSL return code' $sumsl_return
if [ "$sumsl_return" != "4" ]; then
echo 'ERROR = SUMSL return code' $sumsl_return 'is not 4'
cutoffdiff="0.1" # energy cutoff variation - more then this rises warning
sumsl_return=`grep SUMSL $file|awk '{print $4}'`
echo 'SUMSL return code' $sumsl_return
if [ "$sumsl_return" != "4" ]; then
echo 'ERROR = SUMSL return code' $sumsl_return 'is not 4'
fi
elif [ "$1" == "1l2y_MIN_REGULAR_INT" ]; then
extremediff="10.0" # extreme energy difference, comething went terribly wrong
fi
elif [ "$1" == "1l2y_MIN_REGULAR_INT" ]; then
extremediff="10.0" # extreme energy difference, comething went terribly wrong
#
# something wrong with REGULAR and sometimes gives code 8 and different energy
# for now 8 is only warning and cutoffdiff is large
#
# something wrong with REGULAR and sometimes gives code 8 and different energy
# for now 8 is only warning and cutoffdiff is large
if [ `echo "a=${array[0]}-${refe};if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then
echo 'difference ' `echo "a=${array[0]}-${refe};if(0>a)a*=-1;a"|bc -l` "from reference ave etot ${refe} greater than 5.0"
exit 1
if [ `echo "a=${array[0]}-${refe};if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then
echo 'difference ' `echo "a=${array[0]}-${refe};if(0>a)a*=-1;a"|bc -l` "from reference ave etot ${refe} greater than 5.0"
exit 1
- elif [ `echo "a=${array[1]};a>0.15"|bc -l` != "0" ]; then
- echo 'standard deviation greater than 0.15'
+ elif [ `echo "a=${array[1]};a>0.30"|bc -l` != "0" ]; then
+ echo 'standard deviation greater than 0.30'
elif [ "$1" == "1DKZcut-min" ]; then
extremediff="10.0" # extreme energy difference, comething went terribly wrong
elif [ "$1" == "1DKZcut-min" ]; then
extremediff="10.0" # extreme energy difference, comething went terribly wrong
- expectenergy="332.1213" # - expected total energy
- cutoffdiff="3.0" # energy cutoff variation - more then this rises warning
+ expectenergy="323.6931" # - expected total energy
+ cutoffdiff="6.0" # energy cutoff variation - more then this rises warning
elif [ "$1" == "1DKZcut-ene" ]; then
extremediff="10000.0" # extreme energy difference, comething went terribly wrong
elif [ "$1" == "1DKZcut-ene" ]; then
extremediff="10000.0" # extreme energy difference, comething went terribly wrong
cutoffdiff="0.01" # energy cutoff variation - more then this rises warning
elif [ "$1" == "1DKZcut-micro" ]; then
cutoffdiff="0.01" # energy cutoff variation - more then this rises warning
elif [ "$1" == "1DKZcut-micro" ]; then
- extremediff="10.0" # extreme energy difference, comething went terribly wrong
- expectenergy="151.3218" # - expected total energy
- cutoffdiff="6.0" # energy cutoff variation - more then this rises warning
+ extremediff="20.0" # extreme energy difference, comething went terribly wrong
+ expectenergy="141.9976" # - expected total energy
+ cutoffdiff="10.0" # energy cutoff variation - more then this rises warning
startene=`grep ETOT $file|head -1| awk '{print $2*1.0}'`
echo "initial energy=${startene} reference=${refe}"
if [ `echo "a=${startene}-(${refe});if(0>a)a*=-1;a>0.01"|bc -l` != "0" ]; then
startene=`grep ETOT $file|head -1| awk '{print $2*1.0}'`
echo "initial energy=${startene} reference=${refe}"
if [ `echo "a=${startene}-(${refe});if(0>a)a*=-1;a>0.01"|bc -l` != "0" ]; then
diff=`gawk -f checkgrad.awk $file |grep 'Largest abs(1-numerical/analytical)='|awk '{printf "%15.10f",$3}'`
gawk -f checkgrad.awk $file
diff=`gawk -f checkgrad.awk $file |grep 'Largest abs(1-numerical/analytical)='|awk '{printf "%15.10f",$3}'`
gawk -f checkgrad.awk $file
ene=`grep ETOT $file|head -1| awk '{print $2*1.0}'`
hene=`grep H_CONS $file|head -1| awk '{print $2*1.0}'`
ene=`grep ETOT $file|head -1| awk '{print $2*1.0}'`
hene=`grep H_CONS $file|head -1| awk '{print $2*1.0}'`
echo "ETOT= " ${ene} " reference= " ${ene_ref}
if [ `echo "a=${ene}-(${ene_ref});if(0>a)a*=-1;a>0.01"|bc -l` != "0" ]; then
echo "ETOT= " ${ene} " reference= " ${ene_ref}
if [ `echo "a=${ene}-(${ene_ref});if(0>a)a*=-1;a>0.01"|bc -l` != "0" ]; then