@@ -175,7 +175,7 @@ double Energy_Amber::Calc_14_Energy(Frame const& fIn, DihedralArray const& Dihed
175175 double rij = sqrt ( rij2 );
176176 // VDW
177177 NonbondType const & LJ = tIn.GetLJparam (d->A1 (), d->A4 ());
178- double e_vdw = Cpptraj::Energy::Ene_LJ_12_6_4 ( rij2, LJ.A (), LJ.B (), LJ. C () );
178+ double e_vdw = Cpptraj::Energy::Ene_LJ_6_12 ( rij2, LJ.A (), LJ.B () );
179179/* double r2 = 1.0 / rij2;
180180 double r6 = r2 * r2 * r2;
181181 double r12 = r6 * r6;
@@ -235,7 +235,7 @@ double Energy_Amber::E_Nonbond(Frame const& fIn, Topology const& tIn, AtomMask c
235235 double rij = sqrt ( rij2 );
236236 // VDW
237237 NonbondType const & LJ = tIn.GetLJparam (atom1, atom2);
238- double e_vdw = Cpptraj::Energy::Ene_LJ_12_6_4 ( rij2, LJ.A (), LJ.B (), LJ. C () );
238+ double e_vdw = Cpptraj::Energy::Ene_LJ_6_12 ( rij2, LJ.A (), LJ.B () );
239239 /* double r2 = 1.0 / rij2;
240240 double r6 = r2 * r2 * r2;
241241 double r12 = r6 * r6;
@@ -295,7 +295,7 @@ double Energy_Amber::E_VDW(Frame const& fIn, Topology const& tIn, AtomMask const
295295 double rij2 = DIST2_NoImage ( crd1, fIn .XYZ ( atom2 ) );
296296 // VDW
297297 NonbondType const & LJ = tIn.GetLJparam (atom1, atom2);
298- double e_vdw = Cpptraj::Energy::Ene_LJ_12_6_4 ( rij2, LJ.A (), LJ.B (), LJ. C () );
298+ double e_vdw = Cpptraj::Energy::Ene_LJ_6_12 ( rij2, LJ.A (), LJ.B () );
299299 /* double r2 = 1.0 / rij2;
300300 double r6 = r2 * r2 * r2;
301301 double r12 = r6 * r6;
0 commit comments