ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/freemol/trunk/src/mengine/src/ehal.c
(Generate patch)
# Line 194 | Line 194
194                      if (skip[j] == -i)
195                          eps /= units.v14scale;
196                      rik = sqrt(rik2);
197 <                                         rv7 = rv7*rv7*rv7*rv7*rv7*rv7*rv7;  //   pow(rv,7.0);
197 >                    rv7 = rv*rv*rv*rv*rv*rv*rv;  //   pow(rv,7.0);
198                      rv14 = rv7*rv7;
199                      rik6 = rik2*rik2*rik2;
200                      rik7 = rik6*rik;
201                      rho = rik7 + (sigma-1.00)*rv7;
202                      tau = rik + (kappa-1.00)*rv;
203 <                                        tau7 = tau*tau*tau*tau*tau*tau*tau;
204 <                                        tau8 = tau7*tau;
203 >                    tau7 = tau*tau*tau*tau*tau*tau*tau;
204 >                    tau8 = tau7*tau;
205                      e = eps*kappa7*(rv7/tau7)*(sigma*rv7/rho-2.00);
206          
207                      de = -7.00*eps*kappa7*(rv7/tau8)*(sigma*rv7/rho-2.00)

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines