X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=moldyn.c;h=1b52a5628b3e3fe4df32ec4a76483795a1491c15;hb=d445878769b7fcc15959d2d10f2eeced803c5487;hp=dc2dd8e5a5e89226a9e68bd22ccf22bd2ff25635;hpb=593bb07bf38d6d32e9eca61d9637d9e5b7f428d0;p=physik%2Fposic.git

diff --git a/moldyn.c b/moldyn.c
index dc2dd8e..1b52a56 100644
--- a/moldyn.c
+++ b/moldyn.c
@@ -998,6 +998,14 @@ int tersoff(t_moldyn *moldyn) {
 			v3_sub(&dist_ij,btom,&(atom[i]));
 			d_ij=v3_norm(&dist_ij);
 			if(d_ij<=S) {
+
+				S=;
+				R=;
+				A=;
+				lambda=;
+				B=;
+				mu=;
+
 				if(d_ij<=R) {
 					f_c=1.0;
 					df_c=0.0;
@@ -1007,6 +1015,8 @@ int tersoff(t_moldyn *moldyn) {
 					arg1=PI*(d_ij-R)/s_r;
 					f_c=0.5+0.5*cos(arg1);
 					df_c=-0.5*sin(arg1)*(PI/(s_r*d_ij));
+					f_r=A*exp(-lambda*d_ij);
+					f_a=-B*exp(-mu*d_ij);
 				}
 			}
 			else