Fast C routines to compute DFTs (long double precision)
