--- branches/development/src/nonbonded/LJ.cpp 2011/06/09 20:26:29 1579 +++ branches/development/src/nonbonded/LJ.cpp 2011/06/14 20:41:44 1582 @@ -261,9 +261,7 @@ namespace OpenMD { RealType myDeriv = 0.0; RealType myDerivC = 0.0; - ros = *(idat.rij) * sigmai; - - cerr << "ros = " << ros << "\n"; + ros = *(idat.rij) * sigmai; getLJfunc(ros, myPot, myDeriv); @@ -280,17 +278,12 @@ namespace OpenMD { myDerivC = 0.0; } - cerr << "myPot = " << myPot << "\n"; - cerr << "myPotC = " << myPotC << "\n"; - cerr << "epsilon = " << epsilon << "\n"; - cerr << "vdwm = " << *(idat.vdwMult) << "\n"; RealType pot_temp = *(idat.vdwMult) * epsilon * (myPot - myPotC); *(idat.vpair) += pot_temp; RealType dudr = *(idat.sw) * *(idat.vdwMult) * epsilon * (myDeriv - - myDerivC)*sigmai; - - (idat.pot)[VANDERWAALS_FAMILY] += *(idat.sw) * pot_temp; + myDerivC)*sigmai; + (*(idat.pot))[VANDERWAALS_FAMILY] += *(idat.sw) * pot_temp; *(idat.f1) = *(idat.d) * dudr / *(idat.rij); } return;