| 142 |
|
} |
| 143 |
|
} |
| 144 |
|
|
| 145 |
< |
void Morse::calcForce(InteractionData idat) { |
| 145 |
> |
void Morse::calcForce(InteractionData &idat) { |
| 146 |
|
|
| 147 |
|
if (!initialized_) initialize(); |
| 148 |
|
|
| 221 |
|
} |
| 222 |
|
|
| 223 |
|
RealType pot_temp = idat.vdwMult * (myPot - myPotC); |
| 224 |
< |
idat.vpair += pot_temp; |
| 224 |
> |
idat.vpair[0] += pot_temp; |
| 225 |
|
|
| 226 |
|
RealType dudr = idat.sw * idat.vdwMult * (myDeriv - myDerivC); |
| 227 |
|
|
| 228 |
< |
idat.pot += idat.sw * pot_temp; |
| 228 |
> |
idat.pot[0] += idat.sw * pot_temp; |
| 229 |
|
idat.f1 = idat.d * dudr / idat.rij; |
| 230 |
|
} |
| 231 |
|
return; |