35 |
|
* |
36 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
37 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
38 |
< |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
38 |
> |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
39 |
|
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
40 |
|
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
41 |
|
*/ |
60 |
|
} |
61 |
|
|
62 |
|
void InteractionManager::initialize() { |
63 |
< |
|
63 |
> |
|
64 |
> |
if (initialized_) return; |
65 |
> |
|
66 |
|
ForceField* forceField_ = info_->getForceField(); |
67 |
|
|
68 |
|
lj_->setForceField(forceField_); |
111 |
|
|
112 |
|
bool vdwExplicit = false; |
113 |
|
bool metExplicit = false; |
114 |
< |
bool hbExplicit = false; |
114 |
> |
// bool hbExplicit = false; |
115 |
|
|
116 |
|
key = make_pair(atype1, atype2); |
117 |
|
|
304 |
|
eam_->setCutoffRadius(rcut); |
305 |
|
} |
306 |
|
|
305 |
– |
void InteractionManager::setSwitchingRadius(RealType rswitch) { |
306 |
– |
electrostatic_->setSwitchingRadius(rswitch); |
307 |
– |
} |
308 |
– |
|
307 |
|
void InteractionManager::doPrePair(InteractionData idat){ |
308 |
|
|
309 |
|
if (!initialized_) initialize(); |
377 |
|
|
378 |
|
RealType InteractionManager::getSuggestedCutoffRadius(int *atid) { |
379 |
|
if (!initialized_) initialize(); |
380 |
< |
|
380 |
> |
|
381 |
|
AtomType* atype = typeMap_[*atid]; |
382 |
|
|
383 |
|
pair<AtomType*, AtomType*> key = make_pair(atype, atype); |