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). |
39 |
< |
* [4] Vardeman & Gezelter, in progress (2009). |
39 |
> |
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
40 |
> |
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
41 |
|
*/ |
42 |
|
|
43 |
|
/** |
167 |
|
return nAtoms_; |
168 |
|
} |
169 |
|
|
170 |
+ |
/** Returns the number of effective cutoff groups on local processor */ |
171 |
+ |
unsigned int getNLocalCutoffGroups(); |
172 |
+ |
|
173 |
|
/** Returns the number of local bonds */ |
174 |
|
unsigned int getNBonds(){ |
175 |
|
return nBonds_; |
223 |
|
*/ |
224 |
|
Molecule* nextMolecule(MoleculeIterator& i); |
225 |
|
|
226 |
+ |
/** Returns the total number of fluctuating charges that are present */ |
227 |
+ |
int getNFluctuatingCharges() { |
228 |
+ |
return nGlobalFluctuatingCharges_; |
229 |
+ |
} |
230 |
+ |
|
231 |
|
/** Returns the number of degrees of freedom */ |
232 |
|
int getNdf() { |
233 |
|
return ndf_ - getFdf(); |
522 |
|
int nIntegrableObjects_; /**< number of integrable objects in local processor */ |
523 |
|
int nCutoffGroups_; /**< number of cutoff groups in local processor */ |
524 |
|
int nConstraints_; /**< number of constraints in local processors */ |
525 |
+ |
int nFluctuatingCharges_; /**< number of fluctuating charges in local processor */ |
526 |
|
|
527 |
|
/// Counts of global objects |
528 |
|
int nGlobalMols_; /**< number of molecules in the system (GLOBAL) */ |
530 |
|
int nGlobalCutoffGroups_; /**< number of cutoff groups in this system (GLOBAL) */ |
531 |
|
int nGlobalIntegrableObjects_; /**< number of integrable objects in this system */ |
532 |
|
int nGlobalRigidBodies_; /**< number of rigid bodies in this system (GLOBAL) */ |
533 |
+ |
int nGlobalFluctuatingCharges_;/**< number of fluctuating charges in this system (GLOBAL) */ |
534 |
+ |
|
535 |
|
|
536 |
|
/// Degress of freedom |
537 |
|
int ndf_; /**< number of degress of freedom (excludes constraints) (LOCAL) */ |
546 |
|
bool usesDirectionalAtoms_; /**< are there atoms with position AND orientation? */ |
547 |
|
bool usesMetallicAtoms_; /**< are there transition metal atoms? */ |
548 |
|
bool usesElectrostaticAtoms_; /**< are there electrostatic atoms? */ |
549 |
+ |
bool usesFluctuatingCharges_; /**< are there fluctuating charges? */ |
550 |
|
bool usesAtomicVirial_; /**< are we computing atomic virials? */ |
551 |
|
bool requiresPrepair_; /**< does this simulation require a pre-pair loop? */ |
552 |
|
bool requiresSkipCorrection_; /**< does this simulation require a skip-correction? */ |
555 |
|
public: |
556 |
|
bool usesElectrostaticAtoms() { return usesElectrostaticAtoms_; } |
557 |
|
bool usesDirectionalAtoms() { return usesDirectionalAtoms_; } |
558 |
< |
bool usesMetallicAtoms() { return usesMetallicAtoms_; } |
558 |
> |
bool usesFluctuatingCharges() { return usesFluctuatingCharges_; } |
559 |
|
bool usesAtomicVirial() { return usesAtomicVirial_; } |
560 |
|
bool requiresPrepair() { return requiresPrepair_; } |
561 |
|
bool requiresSkipCorrection() { return requiresSkipCorrection_;} |
612 |
|
public: |
613 |
|
vector<RealType> getMassFactors() { return massFactors_; } |
614 |
|
|
615 |
< |
PairList getExcludedInteractions() { return excludedInteractions_; } |
616 |
< |
PairList getOneTwoInteractions() { return oneTwoInteractions_; } |
617 |
< |
PairList getOneThreeInteractions() { return oneThreeInteractions_; } |
618 |
< |
PairList getOneFourInteractions() { return oneFourInteractions_; } |
615 |
> |
PairList* getExcludedInteractions() { return &excludedInteractions_; } |
616 |
> |
PairList* getOneTwoInteractions() { return &oneTwoInteractions_; } |
617 |
> |
PairList* getOneThreeInteractions() { return &oneThreeInteractions_; } |
618 |
> |
PairList* getOneFourInteractions() { return &oneFourInteractions_; } |
619 |
|
|
620 |
|
private: |
608 |
– |
|
621 |
|
|
622 |
|
/// lists to handle atoms needing special treatment in the non-bonded interactions |
623 |
|
PairList excludedInteractions_; /**< atoms excluded from interacting with each other */ |