| 91 |
|
/** |
| 92 |
|
* Constructor of SimInfo |
| 93 |
|
* |
| 94 |
< |
* @param molStampPairs MoleculeStamp Array. The first element of |
| 95 |
< |
* the pair is molecule stamp, the second element is the total |
| 96 |
< |
* number of molecules with the same molecule stamp in the system |
| 94 |
> |
* @param ff pointer to a concrete ForceField instance |
| 95 |
|
* |
| 96 |
< |
* @param ff pointer of a concrete ForceField instance |
| 99 |
< |
* |
| 100 |
< |
* @param simParams |
| 96 |
> |
* @param simParams pointer to the simulation parameters in a Globals object |
| 97 |
|
*/ |
| 98 |
|
SimInfo(ForceField* ff, Globals* simParams); |
| 99 |
|
virtual ~SimInfo(); |
| 104 |
|
* @return return true if adding successfully, return false if the |
| 105 |
|
* molecule is already in SimInfo |
| 106 |
|
* |
| 107 |
< |
* @param mol molecule to be added |
| 107 |
> |
* @param mol Molecule to be added |
| 108 |
|
*/ |
| 109 |
|
bool addMolecule(Molecule* mol); |
| 110 |
|
|
| 282 |
|
Globals* getSimParams() { |
| 283 |
|
return simParams_; |
| 284 |
|
} |
| 289 |
– |
|
| 290 |
– |
/** Returns the velocity of center of mass of the whole system.*/ |
| 291 |
– |
Vector3d getComVel(); |
| 292 |
– |
|
| 293 |
– |
/** Returns the center of the mass of the whole system.*/ |
| 294 |
– |
Vector3d getCom(); |
| 295 |
– |
/** Returns the center of the mass and Center of Mass velocity of |
| 296 |
– |
the whole system.*/ |
| 297 |
– |
void getComAll(Vector3d& com,Vector3d& comVel); |
| 285 |
|
|
| 299 |
– |
/** Returns intertia tensor for the entire system and system |
| 300 |
– |
Angular Momentum.*/ |
| 301 |
– |
void getInertiaTensor(Mat3x3d &intertiaTensor,Vector3d &angularMomentum); |
| 302 |
– |
|
| 303 |
– |
/** Returns system angular momentum */ |
| 304 |
– |
Vector3d getAngularMomentum(); |
| 305 |
– |
|
| 306 |
– |
/** Returns volume of system as estimated by an ellipsoid defined |
| 307 |
– |
by the radii of gyration*/ |
| 308 |
– |
void getGyrationalVolume(RealType &vol); |
| 309 |
– |
/** Overloaded version of gyrational volume that also returns |
| 310 |
– |
det(I) so dV/dr can be calculated*/ |
| 311 |
– |
void getGyrationalVolume(RealType &vol, RealType &detI); |
| 312 |
– |
|
| 286 |
|
void update(); |
| 287 |
|
/** |
| 288 |
|
* Do final bookkeeping before Force managers need their data. |
| 381 |
|
|
| 382 |
|
/** |
| 383 |
|
* Sets GlobalGroupMembership |
| 411 |
– |
* @see #SimCreator::setGlobalIndex |
| 384 |
|
*/ |
| 385 |
|
void setGlobalGroupMembership(const vector<int>& globalGroupMembership) { |
| 386 |
|
assert(globalGroupMembership.size() == static_cast<size_t>(nGlobalAtoms_)); |
| 389 |
|
|
| 390 |
|
/** |
| 391 |
|
* Sets GlobalMolMembership |
| 420 |
– |
* @see #SimCreator::setGlobalIndex |
| 392 |
|
*/ |
| 393 |
|
void setGlobalMolMembership(const vector<int>& globalMolMembership) { |
| 394 |
|
assert(globalMolMembership.size() == static_cast<size_t>(nGlobalAtoms_)); |
| 661 |
|
|
| 662 |
|
/** |
| 663 |
|
* Set MolToProcMap array |
| 693 |
– |
* @see #SimCreator::divideMolecules |
| 664 |
|
*/ |
| 665 |
|
void setMolToProcMap(const vector<int>& molToProcMap) { |
| 666 |
|
molToProcMap_ = molToProcMap; |
