| 68 |
|
using namespace std; |
| 69 |
|
namespace OpenMD { |
| 70 |
|
|
| 71 |
< |
ForceManager::ForceManager(SimInfo * info) : info_(info), switcher_(NULL), |
| 72 |
< |
initialized_(false) { |
| 71 |
> |
ForceManager::ForceManager(SimInfo * info) : info_(info), |
| 72 |
> |
initialized_(false), |
| 73 |
> |
switcher_(NULL) { |
| 74 |
|
forceField_ = info_->getForceField(); |
| 75 |
|
interactionMan_ = new InteractionManager(); |
| 76 |
|
fDecomp_ = new ForceMatrixDecomposition(info_, interactionMan_); |
| 627 |
|
Snapshot* curSnapshot = info_->getSnapshotManager()->getCurrentSnapshot(); |
| 628 |
|
DataStorage* config = &(curSnapshot->atomData); |
| 629 |
|
DataStorage* cgConfig = &(curSnapshot->cgData); |
| 629 |
– |
int jstart, jend; |
| 630 |
|
|
| 631 |
|
//calculate the center of mass of cutoff group |
| 632 |
|
|
| 659 |
|
RealType electroMult, vdwMult; |
| 660 |
|
RealType vij; |
| 661 |
|
Vector3d fij, fg, f1; |
| 662 |
– |
tuple3<RealType, RealType, RealType> cuts; |
| 663 |
– |
RealType rCut, rCutSq, rListSq; |
| 662 |
|
bool in_switching_region; |
| 663 |
|
RealType sw, dswdr, swderiv; |
| 664 |
|
vector<int> atomListColumn, atomListRow; |
| 723 |
|
} |
| 724 |
|
} |
| 725 |
|
|
| 726 |
< |
for (unsigned int cg1 = 0; cg1 < point_.size() - 1; cg1++) { |
| 726 |
> |
for (cg1 = 0; cg1 < point_.size() - 1; cg1++) { |
| 727 |
|
|
| 728 |
|
atomListRow = fDecomp_->getAtomsInGroupRow(cg1); |
| 729 |
|
newAtom1 = true; |