--- trunk/src/primitives/RigidBody.cpp 2005/10/05 19:12:02 646 +++ trunk/src/primitives/RigidBody.cpp 2005/10/12 21:00:59 663 @@ -173,18 +173,16 @@ namespace oopse { IAtom(0, 0) += mtmp * r2; IAtom(1, 1) += mtmp * r2; IAtom(2, 2) += mtmp * r2; + Itmp += IAtom; //project the inertial moment of directional atoms into this rigid body if (atoms_[i]->isDirectional()) { - //IAtom += atoms_[i]->getI(); - Itmp += IAtom; Itmp += refOrients_[i].transpose() * atoms_[i]->getI() * refOrients_[i]; - } else { - Itmp += IAtom; - } + } } - std::cout << Itmp <