Fast C routines to compute DFTs (single and double precision)
