[Chimera-users] Principal axis of inertia

[Tom Goddard]

Oops.  That script I posted for calculating the inertia of a molecule 
had some problems.  The printed matrix of axes was transposed though the 
correct axes were used for the displayed ellipsoid.  Also I was using 
inertia matrix I_ij = sum_over_k(m_k*r_ki*r_kj) omitting the diagonal 
term.  Should have been I_ij = sum_over_k(m_k*(|r_k| - r_ki*r_kj)).  And 
finally I incorrectly computed the ellipsoid size along each axis.  
Everything was a little bit wrong.  Here's the fixed version.  Tested in 
Chimera version 1.2500 but should work with earlier versions.

    Tom

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: inertiapdb.py
URL: <http://plato.cgl.ucsf.edu/pipermail/chimera-users/attachments/20080321/21f01cdf/attachment.ksh>