| 52 |
|
|
| 53 |
|
namespace OpenMD { |
| 54 |
|
|
| 55 |
< |
RepulsivePower::RepulsivePower() : name_("RepulsivePower"), |
| 56 |
< |
initialized_(false), forceField_(NULL) {} |
| 55 |
> |
RepulsivePower::RepulsivePower() : initialized_(false), forceField_(NULL), |
| 56 |
> |
name_("RepulsivePower") {} |
| 57 |
|
|
| 58 |
|
void RepulsivePower::initialize() { |
| 59 |
|
|
| 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); |