| 268 |
|
ndf_local -= nConstraints_; |
| 269 |
|
|
| 270 |
|
#ifdef IS_MPI |
| 271 |
< |
MPI_Allreduce(&ndf_local,&ndf_,1,MPI_INT,MPI_SUM, MPI_COMM_WORLD); |
| 272 |
< |
MPI_Allreduce(&nfq_local,&nGlobalFluctuatingCharges_,1, MPI_INT, MPI_SUM, MPI_COMM_WORLD); |
| 271 |
> |
MPI::COMM_WORLD.Allreduce(&ndf_local, &ndf_, 1, MPI::INT,MPI::SUM); |
| 272 |
> |
MPI::COMM_WORLD.Allreduce(&nfq_local, &nGlobalFluctuatingCharges_, 1, |
| 273 |
> |
MPI::INT, MPI::SUM); |
| 274 |
|
#else |
| 275 |
|
ndf_ = ndf_local; |
| 276 |
|
nGlobalFluctuatingCharges_ = nfq_local; |
| 284 |
|
|
| 285 |
|
int SimInfo::getFdf() { |
| 286 |
|
#ifdef IS_MPI |
| 287 |
< |
MPI_Allreduce(&fdf_local,&fdf_,1,MPI_INT,MPI_SUM, MPI_COMM_WORLD); |
| 287 |
> |
MPI::COMM_WORLD.Allreduce(&fdf_local, &fdf_, 1, MPI::INT, MPI::SUM); |
| 288 |
|
#else |
| 289 |
|
fdf_ = fdf_local; |
| 290 |
|
#endif |
| 340 |
|
} |
| 341 |
|
|
| 342 |
|
#ifdef IS_MPI |
| 343 |
< |
MPI_Allreduce(&ndfRaw_local,&ndfRaw_,1,MPI_INT,MPI_SUM, MPI_COMM_WORLD); |
| 343 |
> |
MPI::COMM_WORLD.Allreduce(&ndfRaw_local, &ndfRaw_, 1, MPI::INT, MPI::SUM); |
| 344 |
|
#else |
| 345 |
|
ndfRaw_ = ndfRaw_local; |
| 346 |
|
#endif |
| 353 |
|
|
| 354 |
|
|
| 355 |
|
#ifdef IS_MPI |
| 356 |
< |
MPI_Allreduce(&ndfTrans_local,&ndfTrans_,1,MPI_INT,MPI_SUM, MPI_COMM_WORLD); |
| 356 |
> |
MPI::COMM_WORLD.Allreduce(&ndfTrans_local, &ndfTrans_, 1, |
| 357 |
> |
MPI::INT, MPI::SUM); |
| 358 |
|
#else |
| 359 |
|
ndfTrans_ = ndfTrans_local; |
| 360 |
|
#endif |
| 780 |
|
#endif |
| 781 |
|
|
| 782 |
|
return atomTypes; |
| 783 |
+ |
} |
| 784 |
+ |
|
| 785 |
+ |
|
| 786 |
+ |
int getGlobalCountOfType(AtomType* atype) { |
| 787 |
+ |
/* |
| 788 |
+ |
set<AtomType*> atypes = getSimulatedAtomTypes(); |
| 789 |
+ |
map<AtomType*, int> counts_; |
| 790 |
+ |
|
| 791 |
+ |
for(mol = beginMolecule(mi); mol != NULL; mol = nextMolecule(mi)) { |
| 792 |
+ |
for(atom = mol->beginAtom(ai); atom != NULL; |
| 793 |
+ |
atom = mol->nextAtom(ai)) { |
| 794 |
+ |
atom->getAtomType(); |
| 795 |
+ |
} |
| 796 |
+ |
} |
| 797 |
+ |
*/ |
| 798 |
+ |
return 0; |
| 799 |
|
} |
| 800 |
|
|
| 801 |
|
void SimInfo::setupSimVariables() { |
| 895 |
|
|
| 896 |
|
|
| 897 |
|
void SimInfo::prepareTopology() { |
| 880 |
– |
int nExclude, nOneTwo, nOneThree, nOneFour; |
| 898 |
|
|
| 899 |
|
//calculate mass ratio of cutoff group |
| 900 |
|
SimInfo::MoleculeIterator mi; |
| 943 |
|
|
| 944 |
|
//scan topology |
| 945 |
|
|
| 929 |
– |
nExclude = excludedInteractions_.getSize(); |
| 930 |
– |
nOneTwo = oneTwoInteractions_.getSize(); |
| 931 |
– |
nOneThree = oneThreeInteractions_.getSize(); |
| 932 |
– |
nOneFour = oneFourInteractions_.getSize(); |
| 933 |
– |
|
| 946 |
|
int* excludeList = excludedInteractions_.getPairList(); |
| 947 |
|
int* oneTwoList = oneTwoInteractions_.getPairList(); |
| 948 |
|
int* oneThreeList = oneThreeInteractions_.getPairList(); |
| 1019 |
|
|
| 1020 |
|
|
| 1021 |
|
StuntDouble* SimInfo::getIOIndexToIntegrableObject(int index) { |
| 1022 |
< |
if (index >= IOIndexToIntegrableObject.size()) { |
| 1022 |
> |
if (index >= int(IOIndexToIntegrableObject.size())) { |
| 1023 |
|
sprintf(painCave.errMsg, |
| 1024 |
|
"SimInfo::getIOIndexToIntegrableObject Error: Integrable Object\n" |
| 1025 |
|
"\tindex exceeds number of known objects!\n"); |
| 1037 |
|
int SimInfo::getNGlobalConstraints() { |
| 1038 |
|
int nGlobalConstraints; |
| 1039 |
|
#ifdef IS_MPI |
| 1040 |
< |
MPI_Allreduce(&nConstraints_, &nGlobalConstraints, 1, MPI_INT, MPI_SUM, |
| 1041 |
< |
MPI_COMM_WORLD); |
| 1040 |
> |
MPI::COMM_WORLD.Allreduce(&nConstraints_, &nGlobalConstraints, 1, |
| 1041 |
> |
MPI::INT, MPI::SUM); |
| 1042 |
|
#else |
| 1043 |
|
nGlobalConstraints = nConstraints_; |
| 1044 |
|
#endif |