[Numpy-discussion] Multi-taper spectral analysis package: pymutt. Critical update.

[martin smith]

The current release, version 0.82.0, contains fixes for two major bugs.
 The first bug is a show-stopping segmentation fault under some versions
of Linux and arises from a variable type mismatch in calls to the numpy
api.  The second bug causes bad spectral values at the Nyquist frequency
for series with even lengths and arises from an error in interfacing to
a local version of fftpack.

I strongly recommend upgrading to the current release.



Pymutt is a module that implements Thomson's
multi-taper spectral analysis algorithms.  The code is based on a
subroutine from Lees and Park and has, of course, a python interface.
References are provided in the readme file.

The code has seen substantial usage and should be fairly reliable.
Examples are included.

It's available at http://code.google.com/p/pymutt/.

Martin L. Smith
Blindgoat Geophysics