Python bindings to a subset of the NUFFT algorithm. Only the 1D/type 3 case is implemented but it's way faster than Lomb-Scargle!
To install, run python setup.py install. Then, there is only one function:
nufft. Assuming that you have a time series in t and y and you
want to evaluate it at (angular) frequencies f:
import nufft
nufft.nufft(t, y, f)You can specify your required precision using eps=1e-10. The default is
1e-10.
Python bindings by Dan Foreman-Mackey but the code that actually does the work is from the Greengard lab at NYU (see the website). The Fortran code is BSD licensed and the Python bindings are MIT licensed.