| 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 |
| 97 |
< |
* |
| 98 |
< |
* @param ff pointer of a concrete ForceField instance |
| 94 |
> |
* @param ff pointer to a concrete ForceField instance |
| 95 |
|
* |
| 96 |
< |
* @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 |
|
|
| 270 |
|
SnapshotManager* getSnapshotManager() { |
| 271 |
|
return sman_; |
| 272 |
|
} |
| 273 |
< |
|
| 273 |
> |
/** Returns the storage layout (computed by SimCreator) */ |
| 274 |
> |
int getStorageLayout() { |
| 275 |
> |
return storageLayout_; |
| 276 |
> |
} |
| 277 |
> |
/** Sets the storage layout (computed by SimCreator) */ |
| 278 |
> |
void setStorageLayout(int sl) { |
| 279 |
> |
storageLayout_ = sl; |
| 280 |
> |
} |
| 281 |
> |
|
| 282 |
|
/** Sets the snapshot manager. */ |
| 283 |
|
void setSnapshotManager(SnapshotManager* sman); |
| 284 |
|
|
| 290 |
|
Globals* getSimParams() { |
| 291 |
|
return simParams_; |
| 292 |
|
} |
| 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); |
| 293 |
|
|
| 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 |
– |
|
| 294 |
|
void update(); |
| 295 |
|
/** |
| 296 |
|
* Do final bookkeeping before Force managers need their data. |
| 389 |
|
|
| 390 |
|
/** |
| 391 |
|
* Sets GlobalGroupMembership |
| 411 |
– |
* @see #SimCreator::setGlobalIndex |
| 392 |
|
*/ |
| 393 |
|
void setGlobalGroupMembership(const vector<int>& globalGroupMembership) { |
| 394 |
|
assert(globalGroupMembership.size() == static_cast<size_t>(nGlobalAtoms_)); |
| 397 |
|
|
| 398 |
|
/** |
| 399 |
|
* Sets GlobalMolMembership |
| 420 |
– |
* @see #SimCreator::setGlobalIndex |
| 400 |
|
*/ |
| 401 |
|
void setGlobalMolMembership(const vector<int>& globalMolMembership) { |
| 402 |
|
assert(globalMolMembership.size() == static_cast<size_t>(nGlobalAtoms_)); |
| 612 |
|
|
| 613 |
|
PropertyMap properties_; /**< Generic Properties can be added */ |
| 614 |
|
SnapshotManager* sman_; /**< SnapshotManager (handles particle positions, etc.) */ |
| 615 |
+ |
int storageLayout_; /**< Bits to tell how much data to store on each object */ |
| 616 |
|
|
| 617 |
|
/** |
| 618 |
|
* The reason to have a local index manager is that when molecule |
| 632 |
|
string dumpFileName_; |
| 633 |
|
string statFileName_; |
| 634 |
|
string restFileName_; |
| 655 |
– |
|
| 635 |
|
|
| 636 |
|
bool topologyDone_; /** flag to indicate whether the topology has |
| 637 |
|
been scanned and all the relevant |
| 669 |
|
|
| 670 |
|
/** |
| 671 |
|
* Set MolToProcMap array |
| 693 |
– |
* @see #SimCreator::divideMolecules |
| 672 |
|
*/ |
| 673 |
|
void setMolToProcMap(const vector<int>& molToProcMap) { |
| 674 |
|
molToProcMap_ = molToProcMap; |