65 |
|
|
66 |
|
public: |
67 |
|
static InteractionManager* Instance(); |
68 |
< |
static void setForceField(ForceField *ff) {forceField_ = ff;}; |
68 |
> |
static void setForceField(ForceField *ff) {forceField_ = ff;} |
69 |
|
|
70 |
+ |
static void setCutoffRadius(RealType rcut) {rCut_ = rcut;} |
71 |
+ |
static void setSwitchingRadius(RealType rsw) {rSwitch_ = rsw;} |
72 |
+ |
static void useShiftedForce() ; |
73 |
+ |
static void useShiftedPot(); |
74 |
|
|
71 |
– |
|
75 |
|
// Fortran support routines |
76 |
|
|
77 |
|
static void doPrePair(int *atid1, int *atid2, RealType *rij, RealType *rho_i_at_j, RealType *rho_j_at_i); |
80 |
|
static void doSkipCorrection(int *atid1, int *atid2, RealType *d, RealType *r, RealType *skippedCharge1, RealType *skippedCharge2, RealType *sw, RealType *electroMult, RealType *pot, RealType *vpair, RealType *f1, RealType *eFrame1, RealType *eFrame2, RealType *t1, RealType *t2); |
81 |
|
static void doSelfCorrection(int *atid, RealType *eFrame, RealType *skippedCharge, RealType *pot, RealType *t); |
82 |
|
static RealType getSuggestedCutoffRadius(int *atid1); |
83 |
+ |
static RealType getSuggestedCutoffRadius(AtomType *atype); |
84 |
|
|
85 |
|
private: |
86 |
|
virtual ~InteractionManager() { } |
102 |
|
static Morse* morse_; |
103 |
|
static Electrostatic* electrostatic_; |
104 |
|
|
105 |
+ |
static RealType rCut_; |
106 |
+ |
static RealType rSwitch_; |
107 |
+ |
|
108 |
|
static map<int, AtomType*> typeMap_; |
109 |
|
/** |
110 |
|
* Each pair of atom types can have multiple interactions, so the |
111 |
|
* natural data structures are a map between the pair, and a set |
112 |
|
* of non-bonded interactions. |
113 |
|
*/ |
114 |
< |
static map<pair<AtomType*, AtomType*>, set<NonBondedInteraction*> > interactions_; |
108 |
< |
|
114 |
> |
static map<pair<AtomType*, AtomType*>, set<NonBondedInteraction*> > interactions_; |
115 |
|
}; |
116 |
|
} |
117 |
|
#endif |