| # | Line 157 | Line 157 | namespace oopse { | |
|---|---|---|
| 157 | ||
| 158 | removeComDrift(); | |
| 159 | // Remove angular drift if we are not using periodic boundary conditions. | |
| 160 | < | if(!simParams->getPBC()) removeAngularDrift(); |
| 160 | > | if(!simParams->getUsePeriodicBoundaryConditions()) removeAngularDrift(); |
| 161 | ||
| 162 | } | |
| 163 | ||
| # | Line 210 | Line 210 | namespace oopse { | |
| 210 | */ | |
| 211 | ||
| 212 | inertiaTensor =inertiaTensor.inverse(); | |
| 213 | + | /* |
| 214 | std::cerr << "Inertia Tensor after inverse is " | |
| 215 | << inertiaTensor << std::endl; | |
| 216 | < | |
| 216 | > | */ |
| 217 | omega = inertiaTensor*angularMomentum; | |
| 218 | ||
| 219 | SimInfo::MoleculeIterator i; | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |