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); |