X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?p=unres.git;a=blobdiff_plain;f=ctest%2Fcsa_energy_check.sh;fp=ctest%2Fcsa_energy_check.sh;h=af2ba4affb2f19092a6eb5ba22be0f21f164e230;hp=0000000000000000000000000000000000000000;hb=9453fc761eb545fcb727824c94d012dbf3931951;hpb=6f521277aa2a382d409f5189957283b0998b0d07

diff --git a/ctest/csa_energy_check.sh b/ctest/csa_energy_check.sh
new file mode 100755
index 0000000..af2ba4a
--- /dev/null
+++ b/ctest/csa_energy_check.sh
@@ -0,0 +1,51 @@
+#!/bin/bash
+
+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
+    echo "CRITICAL: out $file do not exist"
+    exit 2
+fi
+
+emin=`grep jlee $file |awk '{print $10}'`
+echo 'Emin= ' $emin
+
+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
+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