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(); |
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> getExcludesForAtom(int atom1); |
89 |
|
bool skipAtomPair(int atom1, int atom2); |
90 |
+ |
bool excludeAtomPair(int atom1, int atom2); |
91 |
|
void addForceToAtomRow(int atom1, Vector3d fg); |
92 |
|
void addForceToAtomColumn(int atom2, Vector3d fg); |
93 |
|
|
94 |
|
// filling interaction blocks with pointers |
95 |
< |
InteractionData fillInteractionData(int atom1, int atom2); |
96 |
< |
void unpackInteractionData(InteractionData idat, int atom1, int atom2); |
96 |
< |
InteractionData fillSkipData(int atom1, int atom2); |
95 |
> |
void fillInteractionData(InteractionData &idat, int atom1, int atom2); |
96 |
> |
void unpackInteractionData(InteractionData &idat, int atom1, int atom2); |
97 |
|
|
98 |
|
private: |
99 |
|
void createGtypeCutoffMap(); |
100 |
|
|
101 |
|
int nLocal_; |
102 |
|
int nGroups_; |
103 |
– |
vector<int> identsLocal; |
103 |
|
vector<int> AtomLocalToGlobal; |
104 |
|
vector<int> cgLocalToGlobal; |
106 |
– |
vector<RealType> massFactorsLocal; |
107 |
– |
vector<vector<int> > skipsForLocalAtom; |
108 |
– |
vector<vector<int> > toposForLocalAtom; |
109 |
– |
vector<vector<int> > topoDistLocal; |
105 |
|
vector<RealType> groupCutoff; |
106 |
|
vector<int> groupToGtype; |
107 |
|
|
136 |
|
// work arrays for assembling potential energy |
137 |
|
vector<potVec> pot_row; |
138 |
|
vector<potVec> pot_col; |
144 |
– |
potVec pot_local; |
139 |
|
|
140 |
|
vector<int> identsRow; |
141 |
|
vector<int> identsCol; |
162 |
|
|
163 |
|
vector<RealType> massFactorsRow; |
164 |
|
vector<RealType> massFactorsCol; |
171 |
– |
|
172 |
– |
vector<vector<int> > skipsForRowAtom; |
173 |
– |
vector<vector<int> > toposForRowAtom; |
174 |
– |
vector<vector<int> > topoDistRow; |
165 |
|
#endif |
166 |
|
|
167 |
|
}; |