[MMTK] Bug in MMTK.Collections.rotationalConstants()?
Andreas Kring
akring1729 at gmail.com
Tue Feb 12 14:53:28 UTC 2008
> My preferred fix is this:
>
> def rotationalConstants(self, conf=None):
> """Returns a sorted array of rotational constants A, B, C
> in internal units."""
> com, i = self.centerAndMomentOfInertia(conf)
> pmi = i.eigenvalues()
> return Numeric.sort(Units.h / (8.*Numeric.pi*Numeric.pi*pmi))[::-1]
>
> pi comes from Numeric or NumPy (which is imported globally in the module
> Collections), so there is no reason to have another copy in MMTK.Units.
>
Ok, I'll fix it like you show above. Thanks for the quick answer.
Kind regards
Andreas
More information about the mmtk
mailing list