subroutine drfftf (n,r,wsave) c*********************************************************************72 c cc drfftf c double precision r(1), wsave(1) c if (n .eq. 1) return c call drftf1 (n,r,wsave,wsave(n+1),wsave(2*n+1)) c return end