dffti.f90 Source File

      subroutine dffti(n, wsave)
         use fftpack_kind, only: dp => rk
         implicit none
         integer, intent(in) :: n
         real(dp), intent(out) :: wsave(*)
         if (n == 1) return
         call rffti1(n, wsave(n + 1), wsave(2*n + 1))
      end subroutine dffti