| 201 |
|
// electrostaticSummationMethod keyword. |
| 202 |
|
|
| 203 |
|
if (simParams_->haveElectrostaticSummationMethod()) { |
| 204 |
< |
std::string myMethod = simParams_->getElectrostaticSummationMethod(); |
| 204 |
> |
string myMethod = simParams_->getElectrostaticSummationMethod(); |
| 205 |
|
toUpper(myMethod); |
| 206 |
|
|
| 207 |
|
if (myMethod == "SHIFTED_POTENTIAL") { |
| 256 |
|
stringToCutoffPolicy["MAX"] = MAX; |
| 257 |
|
stringToCutoffPolicy["TRADITIONAL"] = TRADITIONAL; |
| 258 |
|
|
| 259 |
< |
std::string cutPolicy; |
| 259 |
> |
string cutPolicy; |
| 260 |
|
if (forceFieldOptions_.haveCutoffPolicy()){ |
| 261 |
|
cutPolicy = forceFieldOptions_.getCutoffPolicy(); |
| 262 |
|
}else if (simParams_->haveCutoffPolicy()) { |
| 388 |
|
setupCutoffs(); |
| 389 |
|
|
| 390 |
|
info_->prepareTopology(); |
| 391 |
+ |
|
| 392 |
+ |
doParticlePot_ = info_->getSimParams()->getOutputParticlePotential(); |
| 393 |
+ |
cerr << "dPP = " << doParticlePot_ << "\n"; |
| 394 |
+ |
|
| 395 |
|
} |
| 396 |
|
|
| 397 |
|
ForceFieldOptions& fopts = forceField_->getForceFieldOptions(); |
| 654 |
|
idat.sw = &sw; |
| 655 |
|
idat.shiftedPot = (cutoffMethod_ == SHIFTED_POTENTIAL) ? true : false; |
| 656 |
|
idat.shiftedForce = (cutoffMethod_ == SHIFTED_FORCE) ? true : false; |
| 657 |
+ |
idat.doParticlePot = doParticlePot_; |
| 658 |
+ |
sdat.doParticlePot = doParticlePot_; |
| 659 |
|
|
| 660 |
|
loopEnd = PAIR_LOOP; |
| 661 |
|
if (info_->requiresPrepair() ) { |
| 860 |
|
MPI_Allreduce(tmpTau.getArrayPointer(), tau.getArrayPointer(), |
| 861 |
|
9, MPI_REALTYPE, MPI_SUM, MPI_COMM_WORLD); |
| 862 |
|
#endif |
| 863 |
< |
curSnapshot->statData.setTau(tau); |
| 863 |
> |
curSnapshot->setTau(tau); |
| 864 |
|
} |
| 865 |
|
|
| 866 |
|
} //end namespace OpenMD |