Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | n | |||
real(kind=dp), | intent(inout) | :: | c(*) | |||
real(kind=dp), | intent(inout) | :: | wsave(*) |
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
integer, | public | :: | iw1 | ||||
integer, | public | :: | iw2 |
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