[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

More information about the mmtk mailing list