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