--- trunk/src/primitives/RigidBody.cpp 2005/10/05 19:35:28 648 +++ trunk/src/primitives/RigidBody.cpp 2005/12/02 15:38:03 770 @@ -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_); @@ -482,7 +483,7 @@ namespace oopse { "RigidBody error.\n" "\tAtom %s does not have a position specified.\n" "\tThis means RigidBody cannot set up reference coordinates.\n", - ats->getType() ); + ats->getType().c_str() ); painCave.isFatal = 1; simError(); } @@ -502,7 +503,7 @@ namespace oopse { "RigidBody error.\n" "\tAtom %s does not have an orientation specified.\n" "\tThis means RigidBody cannot set up reference orientations.\n", - ats->getType() ); + ats->getType().c_str() ); painCave.isFatal = 1; simError(); }