| 53 |  | nObjects_.push_back(info_->getNGlobalBends()); | 
| 54 |  | nObjects_.push_back(info_->getNGlobalTorsions()); | 
| 55 |  | nObjects_.push_back(info_->getNGlobalInversions()); | 
| 56 | + | nObjects_.push_back(info_->getNGlobalMolecules()); | 
| 57 |  |  | 
| 58 |  | stuntdoubles_.resize(nObjects_[STUNTDOUBLE]); | 
| 59 |  | bonds_.resize(nObjects_[BOND]); | 
| 60 |  | bends_.resize(nObjects_[BEND]); | 
| 61 |  | torsions_.resize(nObjects_[TORSION]); | 
| 62 |  | inversions_.resize(nObjects_[INVERSION]); | 
| 63 | + | molecules_.resize(nObjects_[MOLECULE]); | 
| 64 |  |  | 
| 65 |  | SimInfo::MoleculeIterator mi; | 
| 66 |  | Molecule::IntegrableObjectIterator ioi; | 
| 82 |  |  | 
| 83 |  | for (mol = info_->beginMolecule(mi); mol != NULL; | 
| 84 |  | mol = info_->nextMolecule(mi)) { | 
| 85 | < |  | 
| 85 | > |  | 
| 86 | > | molecules_[mol->getGlobalIndex()] = mol; | 
| 87 | > |  | 
| 88 |  | // Hull is constructed from all known integrable objects. | 
| 89 |  | for (sd = mol->beginIntegrableObject(ioi); | 
| 90 |  | sd != NULL; |