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