[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