| 601 |
|
for(cutoffAtom = myCutoffGroup->beginAtom(iterAtom); cutoffAtom != NULL; |
| 602 |
|
cutoffAtom = myCutoffGroup->nextAtom(iterAtom)){ |
| 603 |
|
mfact.push_back(cutoffAtom->getMass()/totalMass); |
| 604 |
+ |
#ifdef IS_MPI |
| 605 |
+ |
groupList.push_back(cutoffAtom->getGlobalIndex() + 1); |
| 606 |
+ |
#else |
| 607 |
|
groupList.push_back(cutoffAtom->getIndex() + 1); |
| 608 |
+ |
#endif |
| 609 |
|
} |
| 610 |
|
|
| 611 |
|
groupStart.push_back(curIndex); |
| 617 |
|
|
| 618 |
|
|
| 619 |
|
//The last cutoff group need more element to indicate the end of the cutoff |
| 620 |
< |
groupStart.push_back(curIndex); |
| 617 |
< |
ngroup = groupStart.size() - 1; |
| 620 |
> |
ngroup = groupStart.size(); |
| 621 |
|
} |