Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | n | |||
real(kind=dp), | intent(inout) | :: | x(*) | |||
real(kind=dp), | intent(in) | :: | wsave(*) |
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
integer, | public | :: | iw1 | ||||
integer, | public | :: | iw2 | ||||
integer, | public | :: | iw3 | ||||
integer, | public | :: | np1 |
subroutine dsint(n, x, wsave)
use fftpack_kind, only: dp => rk
implicit none
integer, intent(in) :: n
real(dp), intent(inout) :: x(*)
real(dp), intent(in) :: wsave(*)
integer :: iw1, iw2, iw3, np1
np1 = n + 1
iw1 = n/2 + 1
iw2 = iw1 + np1
iw3 = iw2 + np1
call sint1(n, x, wsave, wsave(iw1), wsave(iw2), wsave(iw3))
end subroutine dsint