| 55 |  |  | 
| 56 |  | class ForceMatrixDecomposition : public ForceDecomposition { | 
| 57 |  | public: | 
| 58 | < | ForceMatrixDecomposition(SimInfo* info) : ForceDecomposition(info) {}; | 
| 58 | > | ForceMatrixDecomposition(SimInfo* info, InteractionManager* iMan) : ForceDecomposition(info, iMan) {}; | 
| 59 |  |  | 
| 60 |  | void distributeInitialData(); | 
| 61 |  | void zeroWorkArrays(); | 
| 68 |  | vector<pair<int, int> >  buildNeighborList(); | 
| 69 |  |  | 
| 70 |  | // group bookkeeping | 
| 71 | < | pair<int, int> getGroupTypes(int cg1, int cg2); | 
| 71 | > | groupCutoffs getGroupCutoffs(int cg1, int cg2); | 
| 72 |  |  | 
| 73 |  | // Group->atom bookkeeping | 
| 74 |  | vector<int> getAtomsInGroupRow(int cg1); | 
| 84 |  |  | 
| 85 |  | // atom bookkeeping | 
| 86 |  | int getNAtomsInRow(); | 
| 87 | < | vector<int> getSkipsForRowAtom(int atom1); | 
| 88 | < | int getTopoDistance(int atom1, int atom2); | 
| 87 | > | int getTopologicalDistance(int atom1, int atom2); | 
| 88 | > | vector<int> getSkipsForAtom(int atom1); | 
| 89 |  | bool skipAtomPair(int atom1, int atom2); | 
| 90 |  | void addForceToAtomRow(int atom1, Vector3d fg); | 
| 91 |  | void addForceToAtomColumn(int atom2, Vector3d fg); | 
| 96 |  | InteractionData fillSkipData(int atom1, int atom2); | 
| 97 |  |  | 
| 98 |  | private: | 
| 99 | + | void createGtypeCutoffMap(); | 
| 100 | + |  | 
| 101 |  | int nLocal_; | 
| 102 |  | int nGroups_; | 
| 103 |  | vector<int> identsLocal; | 
| 104 |  | vector<int> AtomLocalToGlobal; | 
| 105 |  | vector<int> cgLocalToGlobal; | 
| 106 |  | vector<RealType> massFactorsLocal; | 
| 107 | < | vector<vector<int> > skipsForLocalAtom; | 
| 108 | < | vector<vector<int> > toposForLocalAtom; | 
| 107 | < | vector<vector<int> > topoDistLocal; | 
| 107 | > | vector<RealType> groupCutoff; | 
| 108 | > | vector<int> groupToGtype; | 
| 109 |  |  | 
| 110 |  | #ifdef IS_MPI | 
| 111 |  | DataStorage atomRowData; | 
| 152 |  | vector<int> cgRowToGlobal; | 
| 153 |  | vector<int> cgColToGlobal; | 
| 154 |  |  | 
| 155 | + | vector<RealType> groupCutoffRow; | 
| 156 | + | vector<RealType> groupCutoffCol; | 
| 157 | + | vector<int> groupColToGtype; | 
| 158 | + | vector<int> groupRowToGtype; | 
| 159 | + |  | 
| 160 |  | vector<vector<int> > cellListRow_; | 
| 161 |  | vector<vector<int> > cellListCol_; | 
| 162 |  |  | 
| 165 |  |  | 
| 166 |  | vector<RealType> massFactorsRow; | 
| 167 |  | vector<RealType> massFactorsCol; | 
| 162 | – |  | 
| 163 | – | vector<vector<int> > skipsForRowAtom; | 
| 164 | – | vector<vector<int> > toposForRowAtom; | 
| 165 | – | vector<vector<int> > topoDistRow; | 
| 168 |  | #endif | 
| 169 |  |  | 
| 170 |  | }; |