| 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 |
|
} |