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