--- trunk/OOPSE-3.0/src/primitives/RigidBody.cpp 2005/10/05 19:34:01 2346 +++ trunk/OOPSE-3.0/src/primitives/RigidBody.cpp 2005/10/05 19:35:28 2347 @@ -173,18 +173,15 @@ 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 <