X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=ctest%2Fcsa_energy_check.sh;h=af2ba4affb2f19092a6eb5ba22be0f21f164e230;hb=HEAD;hp=56507fb0899a14fcf8369f4afed6422e5e19e1ef;hpb=aac701bc87a599b2d57e90641bd288fbd5bb368f;p=unres.git diff --git a/ctest/csa_energy_check.sh b/ctest/csa_energy_check.sh index 56507fb..af2ba4a 100755 --- a/ctest/csa_energy_check.sh +++ b/ctest/csa_energy_check.sh @@ -2,6 +2,15 @@ if [ "$1" == "1l2y_csa_4P" ]; then file=1l2y_csa_4P.CSA.history +elif [ "$1" == "1l2y_csa_local_4P" ]; then + file=1l2y_csa_local_4P.CSA.history +elif [ "$1" == "1l2y_csa" ]; then + file=1l2y_csa.CSA.history +elif [ "$1" == "1l2y_csa_local" ]; then + file=1l2y_csa_local.CSA.history +else + exit 1 +fi # Check if file exist if [ ! -f $file ]; then @@ -11,10 +20,32 @@ fi emin=`grep jlee $file |awk '{print $10}'` echo 'Emin= ' $emin -echo 'difference ' `echo "a=$emin+215.41208;if(0>a)a*=-1;a"|bc -l` ' from reference emin -215.41208' -if [ `echo "a=$emin+215.41208;if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then + +if [ "$1" == "1l2y_csa_4P" ]; then + echo 'difference ' `echo "a=$emin+215.41208;if(0>a)a*=-1;a"|bc -l` ' from reference emin -215.41208' + if [ `echo "a=$emin+215.41208;if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then echo 'ERROR difference from reference emin -215.41208 greater than 5.0' exit 1 -fi + fi +elif [ "$1" == "1l2y_csa_local_4P" ]; then + echo 'difference ' `echo "a=$emin+200.13486;if(0>a)a*=-1;a"|bc -l` ' from reference emin -200.13486' + grep rms 1l2y_csa_local_4P*.pdb + if [ `echo "a=$emin+200.13486;if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then + echo 'ERROR difference from reference emin 200.13486 greater than 5.0' + exit 1 + fi +elif [ "$1" == "1l2y_csa" ]; then + echo 'difference ' `echo "a=$emin-29.92824;if(0>a)a*=-1;a"|bc -l` ' from reference emin 29.92824' + if [ `echo "a=$emin-29.92824;if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then + echo 'ERROR difference from reference emin 29.92824 greater than 5.0' + exit 1 + fi +elif [ "$1" == "1l2y_csa_local" ]; then + echo 'difference ' `echo "a=$emin-31.52677;if(0>a)a*=-1;a"|bc -l` ' from reference emin 31.52677' + grep rms 1l2y_csa_local_4P*.pdb + if [ `echo "a=$emin-31.52677;if(0>a)a*=-1;a>5.0"|bc -l` != "0" ]; then + echo 'ERROR difference from reference emin 31.52677 greater than 5.0' + exit 1 + fi fi