| 425 |
|
//find the storage layout |
| 426 |
|
|
| 427 |
|
int storageLayout = computeStorageLayout(info); |
| 428 |
– |
|
| 429 |
– |
cerr << "computed Storage Layout = " << storageLayout << "\n"; |
| 428 |
|
|
| 429 |
|
//allocate memory for DataStorage(circular reference, need to |
| 430 |
|
//break it) |
| 744 |
|
storageLayout |= DataStorage::dslFlucQForce; |
| 745 |
|
} |
| 746 |
|
} |
| 747 |
+ |
|
| 748 |
+ |
// if the user has asked for them, make sure we've got the memory for the |
| 749 |
+ |
// objects defined. |
| 750 |
|
|
| 751 |
|
if (simParams->getOutputParticlePotential()) { |
| 752 |
|
storageLayout |= DataStorage::dslParticlePot; |
| 753 |
|
} |
| 754 |
+ |
if (simParams->getOutputElectricField()) { |
| 755 |
+ |
storageLayout |= DataStorage::dslElectricField; |
| 756 |
+ |
} |
| 757 |
+ |
if (simParams->getOutputFluctuatingCharges()) { |
| 758 |
+ |
storageLayout |= DataStorage::dslFlucQPosition; |
| 759 |
+ |
storageLayout |= DataStorage::dslFlucQVelocity; |
| 760 |
+ |
storageLayout |= DataStorage::dslFlucQForce; |
| 761 |
+ |
} |
| 762 |
|
|
| 763 |
|
return storageLayout; |
| 764 |
|
} |