Fast C routines to compute the Non-equispaced Discrete Fourier Transform
