dcosti Subroutine

subroutine dcosti(n, Wsave)

Uses

Arguments

Type IntentOptional Attributes Name
integer :: n
real(kind=rk) :: Wsave

Contents

Variables

Source Code

dcosti

Variables

Type Visibility Attributes Name Initial
real(kind=rk), public :: dt
real(kind=rk), public :: fk
integer, public :: k
integer, public :: kc
integer, public :: nm1
integer, public :: np1
integer, public :: ns2
real(kind=rk), public, parameter :: pi = acos(-1.0_rk)

Source Code

      subroutine dcosti(n,Wsave)
      use fftpack_kind
      implicit none
      real(rk) :: dt , fk , Wsave
      integer :: k , kc , n , nm1 , np1 , ns2
      dimension Wsave(*)
      real(rk),parameter :: pi = acos(-1.0_rk)
      if ( n<=3 ) return
      nm1 = n - 1
      np1 = n + 1
      ns2 = n/2
      dt = pi/real(nm1, rk)
      fk = 0.0_rk
      do k = 2 , ns2
         kc = np1 - k
         fk = fk + 1.0_rk
         Wsave(k) = 2.0_rk*sin(fk*dt)
         Wsave(kc) = 2.0_rk*cos(fk*dt)
      enddo
      call dffti(nm1,Wsave(n+1))
      end subroutine dcosti