@@ -96,7 +96,7 @@ void PotentialTerm_LJ_Coulomb::CalcForce(Frame& frameIn, CharMask const& maskIn)
9696 {
9797 NonbondType LJ = GetLJparam (*atoms_, *nonbond_, *idx, *jdx);
9898 // NonBondKernel(frameIn, *idx, *jdx, LJ.A(), LJ.B(), maskIn, *E_vdw_, *E_elec_);
99- nonbond.Calc_F_E (frameIn, *idx, *jdx, LJ.A (), LJ.B (), LJ. C (), QFAC_, (*atoms_)[*idx].Charge (), (*atoms_)[*jdx].Charge (), 1 , 1 , maskIn, *E_vdw_, *E_elec_);
99+ nonbond.Calc_F_E (frameIn, *idx, *jdx, LJ.A (), LJ.B (), QFAC_, (*atoms_)[*idx].Charge (), (*atoms_)[*jdx].Charge (), 1 , 1 , maskIn, *E_vdw_, *E_elec_);
100100 } // END idx not bonded to jdx
101101 } // END inner loop over jdx
102102 } // END outer loop over idx
@@ -112,7 +112,7 @@ void PotentialTerm_LJ_Coulomb::CalcForce(Frame& frameIn, CharMask const& maskIn)
112112 if (excluded.find ( *idx1 ) == excluded.end ())
113113 {
114114 NonbondType LJ = GetLJparam (*atoms_, *nonbond_, *idx0, *idx1);
115- nonbond.Calc_F_E (frameIn, *idx0, *idx1, LJ.A (), LJ.B (), LJ. C (), QFAC_, (*atoms_)[*idx0].Charge (), (*atoms_)[*idx1].Charge (), 1 , 1 , maskIn, *E_vdw_, *E_elec_);
115+ nonbond.Calc_F_E (frameIn, *idx0, *idx1, LJ.A (), LJ.B (), QFAC_, (*atoms_)[*idx0].Charge (), (*atoms_)[*idx1].Charge (), 1 , 1 , maskIn, *E_vdw_, *E_elec_);
116116 // NonBondKernel(frameIn, *idx0, *idx1, LJ.A(), LJ.B(), maskIn, *E_vdw_, *E_elec_);
117117 }
118118 } // END inner loop over idx1
0 commit comments