X-Git-Url: https://hackdaworld.org/gitweb/?p=physik%2Fposic.git;a=blobdiff_plain;f=calc_delta_e;fp=calc_delta_e;h=eff07c93e2eb870a97abe4ab838646ec0774049e;hp=2b4130e2ff8b8a033ac6d3c497d07579080a8e7b;hb=d24930a9135ffd5927fbf65a71494f5cf772f264;hpb=8524173a28f2c22a539ef1b0910a1136d9cb254b

diff --git a/calc_delta_e b/calc_delta_e
index 2b4130e..eff07c9 100755
--- a/calc_delta_e
+++ b/calc_delta_e
@@ -1,8 +1,21 @@
 #!/bin/bash
 
+# guess potential
+pot=`grep ^potential $1/config | awk '{ print $2 }'`
+
+echo
+if [ "$pot" = "tersoff" ]; then
+echo "potential: $pot"
+musi=-4.629595
+muc=-7.421824
+music=-12.37363
+else
+echo "potential: $pot"
 musi=-4.628414
 muc=-7.373091
 music=-12.679618
+fi
+echo
 
 file=`ls $1/atomic_conf_* | tail -1`
 atom_cnt=`grep '# \[P\]' $file | awk '{ print $3 }'`
@@ -21,5 +34,7 @@ echo "$si_cnt $c_cnt $e1 $e0" | \
 	awk '{ print "   Gao: "($3-$4)*($1+$2) }'
 echo "$si_cnt $c_cnt $musi $muc $e1 $music" | \
 	#awk '{ print "   Posselt: "$5*($1+$2)-0.5*($1+$2)*$6-0.5*($1-$2)*($3-$4) }'
-	awk '{ print "   Posselt: "$5*($1+$2)-$1*$3-$2*$4) }'
+	awk '{ print "   Posselt: "$5*($1+$2)-$1*$3-$2*$4 }'
+echo "$si_cnt $c_cnt $musi $music $e1 $music" | \
+	awk '{ print "   Tersoff: "$5*($1+$2)-$1*$3-$2*($4-$3) }'