dcosqf Subroutine

subroutine dcosqf(n, x, wsave)

Uses

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: n
real(kind=dp), intent(inout) :: x(*)
real(kind=dp), intent(in) :: wsave(*)

Contents

Variables

Source Code

dcosqf

Variables

Type Visibility Attributes Name Initial
real(kind=dp), public, parameter :: sqrt2 = sqrt(2.0_dp)
real(kind=dp), public :: tsqx

Source Code

      subroutine dcosqf(n, x, wsave)
         use fftpack_kind, only: dp => rk
         implicit none
         integer, intent(in) :: n
         real(dp), intent(in) :: wsave(*)
         real(dp), intent(inout) :: x(*)
         real(dp) :: tsqx
         real(dp), parameter :: sqrt2 = sqrt(2.0_dp)
         if (n < 2) then
            return
         elseif (n == 2) then
            tsqx = sqrt2*x(2)
            x(2) = x(1) - tsqx
            x(1) = x(1) + tsqx
         else
            call cosqf1(n, x, wsave, wsave(n + 1))
         end if
      end subroutine dcosqf