| 132 |
|
//equal to the total number of atoms minus number of atoms belong to |
| 133 |
|
//cutoff group defined in meta-data file plus the number of cutoff |
| 134 |
|
//groups defined in meta-data file |
| 135 |
+ |
std::cerr << "nGA = " << nGlobalAtoms_ << "\n"; |
| 136 |
+ |
std::cerr << "nCA = " << nCutoffAtoms << "\n"; |
| 137 |
+ |
std::cerr << "nG = " << nGroups << "\n"; |
| 138 |
+ |
|
| 139 |
|
nGlobalCutoffGroups_ = nGlobalAtoms_ - nCutoffAtoms + nGroups; |
| 140 |
+ |
|
| 141 |
+ |
std::cerr << "nGCG = " << nGlobalCutoffGroups_ << "\n"; |
| 142 |
|
|
| 143 |
|
//every free atom (atom does not belong to rigid bodies) is an |
| 144 |
|
//integrable object therefore the total number of integrable objects |
| 976 |
|
Molecule* mol; |
| 977 |
|
RigidBody* rb; |
| 978 |
|
Atom* atom; |
| 979 |
+ |
CutoffGroup* cg; |
| 980 |
|
SimInfo::MoleculeIterator mi; |
| 981 |
|
Molecule::RigidBodyIterator rbIter; |
| 982 |
< |
Molecule::AtomIterator atomIter;; |
| 982 |
> |
Molecule::AtomIterator atomIter; |
| 983 |
> |
Molecule::CutoffGroupIterator cgIter; |
| 984 |
|
|
| 985 |
|
for (mol = beginMolecule(mi); mol != NULL; mol = nextMolecule(mi)) { |
| 986 |
|
|
| 991 |
|
for (rb = mol->beginRigidBody(rbIter); rb != NULL; rb = mol->nextRigidBody(rbIter)) { |
| 992 |
|
rb->setSnapshotManager(sman_); |
| 993 |
|
} |
| 994 |
+ |
|
| 995 |
+ |
for (cg = mol->beginCutoffGroup(cgIter); cg != NULL; cg = mol->nextCutoffGroup(cgIter)) { |
| 996 |
+ |
cg->setSnapshotManager(sman_); |
| 997 |
+ |
} |
| 998 |
|
} |
| 999 |
|
|
| 1000 |
|
} |
