Move calculation of spectrum to utils
Use fftw