X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=potentials%2Falbe.c;fp=potentials%2Falbe.c;h=9d8336e7ca99e7b1fedcaba8803aa8d735196f9a;hb=520d9dd22e61fb0b303733421d760da609dd1b3c;hp=e8aee74d4be12cf9abc3821acb5abb064b12fbc2;hpb=2076649dc4067af11bfeabc4ef3f353420b3df14;p=physik%2Fposic.git

diff --git a/potentials/albe.c b/potentials/albe.c
index e8aee74..9d8336e 100644
--- a/potentials/albe.c
+++ b/potentials/albe.c
@@ -448,7 +448,7 @@ printf("    t: %.15f %.15f %.15f\n",aj->f.x,aj->f.y,aj->f.z);
 #endif
 
 	/* virial */
-	virial_calc(aj,&force,dist);
+	virial_calc(ai,&force,dist);
 
 	/* dzeta prefactor = - f_c f_a db, (* -0.5 due to force calc) */
 	exchange->pre_dzeta=0.5*f_a*f_c*db;
@@ -464,6 +464,9 @@ printf("    t: %.15f %.15f %.15f\n",aj->f.x,aj->f.y,aj->f.z);
 }
 #endif
 
+	/* virial */
+	virial_calc(ai,&force,dist);
+
 	v3_scale(&force,&force,-1.0);
 	v3_add(&(ai->f),&(ai->f),&force);
 
@@ -474,9 +477,6 @@ printf("    t: %.15f %.15f %.15f\n",ai->f.x,ai->f.y,ai->f.z);
 }
 #endif
 
-	/* virial */
-	virial_calc(ai,&force,dist);
-
 	/* reset k counter for second k loop */
 	exchange->kcnt=0;
 		
@@ -520,6 +520,9 @@ printf("    t: %.15f %.15f %.15f\n",ak->f.x,ak->f.y,ak->f.z);
 }
 #endif
 
+	/* virial */
+	virial_calc(ai,&force,&(exchange->dist[k]));
+
 	v3_scale(&force,&force,-1.0);
 	v3_add(&(ai->f),&(ai->f),&force);
 
@@ -531,9 +534,6 @@ printf("  ## %f\n",exchange->d[k]);
 }
 #endif
 
-	/* virial */
-	virial_calc(ai,&force,&(exchange->dist[k]));
-
 	/* increase k counter */
 	exchange->kcnt+=1;