X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;ds=inline;f=ctest%2Fprota_unres_energy_check.sh;h=8497f9da2c8ddba71a8852a21394e477af63c10d;hb=6f0e7ef07e02a56fff97df63bd7f07d4b6fb50ce;hp=88ded4bd13fe4fd24d971550aec095fa9445f053;hpb=f51580b07599cf085ce9867770100304ff4f371b;p=unres.git diff --git a/ctest/prota_unres_energy_check.sh b/ctest/prota_unres_energy_check.sh index 88ded4b..8497f9d 100755 --- a/ctest/prota_unres_energy_check.sh +++ b/ctest/prota_unres_energy_check.sh @@ -65,13 +65,17 @@ elif [ "$1" == "1l2y_micro" ]; then 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.1"|bc -l` != "0" ]; then - echo 'standard deviation greater than 0.1' + elif [ `echo "a=${array[1]};a>0.15"|bc -l` != "0" ]; then + echo 'standard deviation greater than 0.15' exit 1 else exit 0 fi elif [ "$1" == "1L2Y_L" ] || [ "$1" == "1L2Y_NH" ]; then + if [ ! -s $file_stat ]; then + echo 'FATAL error - stat file empty' + exit 2 + fi chi2=`./matplotlib_fit_hist.py $file_stat 111` echo 'Chi2 for fitting theoretical temperature distribution ' ${chi2} echo "`pwd`/${file_stat}.png" @@ -83,6 +87,10 @@ elif [ "$1" == "1L2Y_L" ] || [ "$1" == "1L2Y_NH" ]; then exit 0 fi elif [ "$1" == "1L2Y_B" ]; then + if [ ! -s $file_stat ]; then + echo 'FATAL error - stat file empty' + exit 2 + fi chi2=`./matplotlib_fit_hist.py $file_stat 111` echo 'Chi2 for fitting theoretical temperature distribution ' ${chi2} echo "`pwd`/${file_stat}.png" @@ -98,6 +106,10 @@ elif [ "$1" == "1L2Y_B" ]; then elif [ "$1" == "1L2Y_remd" ]; then rm -rf remd_all.stat tail -q -n +1000 1L2Y_remd*.stat >remd_all.stat + if [ ! -s remd_all.stat ]; then + echo 'FATAL error - stat files empty' + exit 2 + fi ./matplotlib_hist.py echo "`pwd`/1L2Y_remd_ene_hist.png" @@ -173,7 +185,7 @@ elif [ "$1" == "1DKZcut-micro" ]; then elif [ "$1" == "1ei0_min" ]; then extremediff="10.0" # extreme energy difference, comething went terribly wrong expectenergy="151.3218" # - expected total energy - cutoffdiff="5.0" # energy cutoff variation - more then this rises warning + cutoffdiff="6.0" # energy cutoff variation - more then this rises warning refe="134.8382" startene=`grep ETOT $file|head -1| awk '{print $2*1.0}'`