| 139 |
|
void Morse::calcForce(InteractionData &idat) { |
| 140 |
|
|
| 141 |
|
if (!initialized_) initialize(); |
| 142 |
< |
|
| 142 |
> |
|
| 143 |
> |
|
| 144 |
|
MorseInteractionData &mixer = MixingMap[Mtids[idat.atid1]][Mtids[idat.atid2]]; |
| 145 |
|
|
| 146 |
|
RealType myPot = 0.0; |
| 172 |
|
|
| 173 |
|
switch(variant) { |
| 174 |
|
case mtShifted : { |
| 174 |
– |
|
| 175 |
|
myPot = De * (expfnc2 - 2.0 * expfnc); |
| 176 |
|
myDeriv = 2.0 * De * beta * (expfnc - expfnc2); |
| 177 |
|
|
| 190 |
|
break; |
| 191 |
|
} |
| 192 |
|
case mtRepulsive : { |
| 193 |
– |
|
| 193 |
|
myPot = De * expfnc2; |
| 194 |
|
myDeriv = -2.0 * De * beta * expfnc2; |
| 195 |
|
|
| 222 |
|
|
| 223 |
|
(*(idat.pot))[VANDERWAALS_FAMILY] += *(idat.sw) * pot_temp; |
| 224 |
|
*(idat.f1) = *(idat.d) * dudr / *(idat.rij); |
| 225 |
< |
|
| 225 |
> |
|
| 226 |
|
return; |
| 227 |
|
} |
| 228 |
|
|