| # | Line 177 | Line 177 | namespace oopse { | |
|---|---|---|
| 177 | ||
| 178 | //project the inertial moment of directional atoms into this rigid body | |
| 179 | if (atoms_[i]->isDirectional()) { | |
| 180 | – | Itmp += IAtom; |
| 180 | Itmp += refOrients_[i].transpose() * atoms_[i]->getI() * refOrients_[i]; | |
| 181 | } | |
| 182 | } | |
| 183 | ||
| 184 | + | std::cout << Itmp << std::endl; |
| 185 | + | |
| 186 | //diagonalize | |
| 187 | Vector3d evals; | |
| 188 | Mat3x3d::diagonalize(Itmp, evals, sU_); | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |