subroutine cfft1d(n,c,dir) integer*4 n, dir, ier complex*8 c(*) call fftww(n,c,dir) return end