| 182 |
|
for (mol = info_->beginMolecule(i); mol != NULL; mol = info_->nextMolecule(i)) { |
| 183 |
|
for (integrableObject = mol->beginIntegrableObject(j); integrableObject != NULL; |
| 184 |
|
integrableObject = mol->nextIntegrableObject(j)) { |
| 185 |
– |
oldVel_[index] = integrableObject->getVel(); |
| 186 |
– |
oldJi_[index] = integrableObject->getJ(); |
| 185 |
|
|
| 186 |
+ |
oldVel_[index] = integrableObject->getVel(); |
| 187 |
+ |
|
| 188 |
+ |
if (integrableObject->isDirectional()) |
| 189 |
+ |
oldJi_[index] = integrableObject->getJ(); |
| 190 |
+ |
|
| 191 |
|
++index; |
| 192 |
< |
} |
| 190 |
< |
|
| 192 |
> |
} |
| 193 |
|
} |
| 194 |
|
|
| 195 |
|
// do the iteration: |