45 |
|
#include "brains/SimInfo.hpp" |
46 |
|
#include "brains/SnapshotManager.hpp" |
47 |
|
#include "nonbonded/NonBondedInteraction.hpp" |
48 |
+ |
#include "nonbonded/Cutoffs.hpp" |
49 |
+ |
#include "nonbonded/InteractionManager.hpp" |
50 |
+ |
#include "utils/Tuple.hpp" |
51 |
|
|
52 |
|
using namespace std; |
53 |
|
namespace OpenMD { |
54 |
|
|
55 |
+ |
typedef tuple3<RealType, RealType, RealType> groupCutoffs; |
56 |
+ |
|
57 |
|
/** |
58 |
|
* @class ForceDecomposition |
59 |
|
* |
106 |
|
virtual bool checkNeighborList(); |
107 |
|
virtual vector<pair<int, int> > buildNeighborList() = 0; |
108 |
|
|
109 |
+ |
// how to handle cutoffs: |
110 |
+ |
void setCutoffPolicy(CutoffPolicy cp) {cutoffPolicy_ = cp;} |
111 |
+ |
void setUserCutoff(RealType rcut) {userCutoff_ = rcut; userChoseCutoff_ = true; } |
112 |
+ |
|
113 |
|
// group bookkeeping |
114 |
< |
virtual pair<int, int> getGroupTypes(int cg1, int cg2) = 0; |
114 |
> |
virtual groupCutoffs getGroupCutoffs(int cg1, int cg2) = 0; |
115 |
|
|
116 |
|
// Group->atom bookkeeping |
117 |
|
virtual vector<int> getAtomsInGroupRow(int cg1) = 0; |
144 |
|
SnapshotManager* sman_; |
145 |
|
Snapshot* snap_; |
146 |
|
ForceField* ff_; |
147 |
+ |
InteractionManager* interactionMan_; |
148 |
+ |
|
149 |
|
int storageLayout_; |
150 |
|
RealType skinThickness_; /**< Verlet neighbor list skin thickness */ |
151 |
+ |
RealType largestRcut_; |
152 |
|
|
153 |
|
map<pair<int, int>, int> topoDist; //< topoDist gives the |
154 |
|
//topological distance between |
172 |
|
//method to fill this. |
173 |
|
|
174 |
|
vector<vector<int> > groupList_; |
175 |
+ |
|
176 |
|
vector<Vector3i> cellOffsets_; |
177 |
|
Vector3i nCells_; |
178 |
|
vector<vector<int> > cellList_; |
179 |
|
vector<Vector3d> saved_CG_positions_; |
180 |
|
potVec longRangePot_; |
181 |
|
|
182 |
+ |
bool userChoseCutoff_; |
183 |
+ |
RealType userCutoff_; |
184 |
+ |
CutoffPolicy cutoffPolicy_; |
185 |
+ |
|
186 |
+ |
map<pair<int, int>, tuple3<RealType, RealType, RealType> > gTypeCutoffMap; |
187 |
+ |
|
188 |
|
}; |
189 |
|
} |
190 |
|
#endif |