240 |
|
} |
241 |
|
} |
242 |
|
|
243 |
< |
void LJ::calcForce(InteractionData idat) { |
243 |
> |
void LJ::calcForce(InteractionData &idat) { |
244 |
|
|
245 |
|
if (!initialized_) initialize(); |
246 |
|
|
280 |
|
} |
281 |
|
|
282 |
|
RealType pot_temp = idat.vdwMult * epsilon * (myPot - myPotC); |
283 |
< |
idat.vpair += pot_temp; |
283 |
> |
idat.vpair[0] += pot_temp; |
284 |
|
|
285 |
|
RealType dudr = idat.sw * idat.vdwMult * epsilon * (myDeriv - |
286 |
|
myDerivC)*sigmai; |
287 |
|
|
288 |
< |
idat.pot += idat.sw * pot_temp; |
289 |
< |
idat.f1 = idat.d * dudr / idat.rij; |
290 |
< |
|
288 |
> |
idat.pot[0] += idat.sw * pot_temp; |
289 |
> |
idat.f1 = idat.d * dudr / idat.rij; |
290 |
|
} |
291 |
|
return; |
292 |
|
} |