--- trunk/src/primitives/RigidBody.cpp 2005/10/05 19:35:28 648 +++ trunk/src/primitives/RigidBody.cpp 2005/10/12 21:00:59 663 @@ -177,11 +177,12 @@ namespace oopse { //project the inertial moment of directional atoms into this rigid body if (atoms_[i]->isDirectional()) { - Itmp += IAtom; Itmp += refOrients_[i].transpose() * atoms_[i]->getI() * refOrients_[i]; } } + std::cout << Itmp << std::endl; + //diagonalize Vector3d evals; Mat3x3d::diagonalize(Itmp, evals, sU_);