ene=`grep ETOT $file|head -1| awk '{print $2*1.0}'`
hene=`grep H_CONS $file|head -1| awk '{print $2*1.0}'`
- ene_ref=3124.347
+ ene_ref=3124.502
hene_ref=1872.981
echo "ETOT= " ${ene} " reference= " ${ene_ref}
fi
- if [ `echo "a=${diff};a>0.0020"|bc -l` != "0" ]; then
+ if [ `echo "a=${diff};a>0.0080"|bc -l` != "0" ]; then
echo 'ERROR largest abs(1-numerical/analytical)' ${diff}
- echo ' greater than 0.0020'
+ echo ' greater than 0.0080'
exit_error=1
fi
fi
# Check if energy value is extremely different from expected
-blaene=`floating $grepene`
-myene=`echo $blaene |bc`
-diff=`echo $myene"+("$expectenergy")"|bc`
-absdiff=`echo $diff| sed 's/-//'`
+#blaene=`floating $grepene`
+#myene=`echo $blaene |bc`
+#diff=`echo $myene"+("$expectenergy")"|bc`
+#absdiff=`echo $diff| sed 's/-//'`
+absdiff=`awk "function abs(v) {return v < 0 ? -v : v}BEGIN{print $grepene+$expectenergy}"`
true1=`echo "$absdiff >= $extremediff"|bc`
if [ $true1 -eq 1 ]; then