zfftf.f90 Source File

      subroutine zfftf(n, c, wsave)
         use fftpack_kind, only: dp => rk
         implicit none
         integer, intent(in) :: n
         real(dp), intent(inout) :: c(*)
         real(dp), intent(inout) :: wsave(*)
         integer :: iw1, iw2
         if (n == 1) return
         iw1 = n + n + 1
         iw2 = iw1 + n + n
         call cfftf1(n, c, wsave, wsave(iw1), wsave(iw2))
      end subroutine zfftf