| 157 |
|
|
| 158 |
|
removeComDrift(); |
| 159 |
|
// Remove angular drift if we are not using periodic boundary conditions. |
| 160 |
< |
if(simParams->getPBC()) removeAngularDrift(); |
| 160 |
> |
if(!simParams->getPBC()) removeAngularDrift(); |
| 161 |
|
|
| 162 |
|
} |
| 163 |
|
|
| 204 |
|
// We now need the inverse of the inertia tensor. |
| 205 |
|
|
| 206 |
|
std::cerr << "Angular Momentum before is " |
| 207 |
< |
<< angularMomentum << std::endl; |
| 208 |
< |
|
| 207 |
> |
<< angularMomentum << std::endl; |
| 208 |
> |
std::cerr << "Inertia Tensor before is " |
| 209 |
> |
<< inertiaTensor << std::endl; |
| 210 |
|
|
| 210 |
– |
inertiaTensor.inverse(); |
| 211 |
|
|
| 212 |
+ |
inertiaTensor =inertiaTensor.inverse(); |
| 213 |
+ |
std::cerr << "Inertia Tensor after inverse is " |
| 214 |
+ |
<< inertiaTensor << std::endl; |
| 215 |
|
|
| 216 |
|
omega = inertiaTensor*angularMomentum; |
| 217 |
|
|
