Fast C routines to compute DFTs (quad precision)
