Fast C routines to compute the Discrete Fourier Transform