| 174 |
|
*(idat.vpair) += pot_temp; |
| 175 |
|
|
| 176 |
|
RealType dudr = *(idat.sw) * *(idat.vdwMult) * epsilon * (myDeriv - |
| 177 |
< |
myDerivC)*sigmai; |
| 177 |
> |
myDerivC)*sigmai; |
| 178 |
|
|
| 179 |
|
(*(idat.pot))[VANDERWAALS_FAMILY] += *(idat.sw) * pot_temp; |
| 180 |
|
*(idat.f1) = *(idat.d) * dudr / *(idat.rij); |
| 182 |
|
return; |
| 183 |
|
} |
| 184 |
|
|
| 185 |
< |
void RepulsivePower::getNRepulsionFunc(const RealType &r, int &n, RealType &pot, RealType &deriv) { |
| 185 |
> |
void RepulsivePower::getNRepulsionFunc(const RealType &r, int &n, |
| 186 |
> |
RealType &pot, RealType &deriv) { |
| 187 |
|
|
| 188 |
|
RealType ri = 1.0 / r; |
| 189 |
|
RealType rin = pow(ri, n); |