| 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 | 
  | 
} |