X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=potentials%2Flennard_jones.c;h=1277c466c919d98f3137c63dbab487bd2e70bc1e;hb=a70de3dccbf0a01c39c2643818ec86c0b465caba;hp=7d45a177134c13080cb5002ad474117f2450eab7;hpb=caa3bc828974c35df2462fde737c31c0a618ee4e;p=physik%2Fposic.git

diff --git a/potentials/lennard_jones.c b/potentials/lennard_jones.c
index 7d45a17..1277c46 100644
--- a/potentials/lennard_jones.c
+++ b/potentials/lennard_jones.c
@@ -16,8 +16,8 @@
 #include <math.h>
 
 #include "../moldyn.h"
-#inlcude "../math/math.h"
-//#include "lennard_jones.h"
+#include "../math/math.h"
+#include "lennard_jones.h"
 
 int lennard_jones(t_moldyn *moldyn,t_atom *ai,t_atom *aj,u8 bc) {
 
@@ -26,7 +26,7 @@ int lennard_jones(t_moldyn *moldyn,t_atom *ai,t_atom *aj,u8 bc) {
 	double d,h1,h2;
 	double eps,sig6,sig12;
 
-	params=moldyn->pot2b_params;
+	params=moldyn->pot_params;
 	eps=params->epsilon4;
 	sig6=params->sigma6;
 	sig12=params->sigma12;
@@ -35,10 +35,11 @@ int lennard_jones(t_moldyn *moldyn,t_atom *ai,t_atom *aj,u8 bc) {
 
 	v3_sub(&distance,&(aj->r),&(ai->r));
 	if(bc) check_per_bound(moldyn,&distance);
-	d=v3_absolute_square(&distance);	/* 1/r^2 */
+	d=v3_absolute_square(&distance);	/* r^2 */
 	if(d<=moldyn->cutoff_square) {
 		d=1.0/d; 			/* 1/r^2 */
 		h2=d*d;				/* 1/r^4 */
+		h2*=d;				/* 1/r^6 */
 		h1=h2*h2;			/* 1/r^12 */
 		moldyn->energy+=(eps*(sig12*h1-sig6*h2)-params->uc);
 		h2*=d;				/* 1/r^8 */