unswitchreal_cent.f90

来自「基于linux操作下的fortran快速付立变换的程序」· F90 代码 · 共 42 行

        subroutine unswitchreal_cent(n1dfft,max2,n2,lot,n1zt,lzt,zw,zt)
        implicit real*8 (a-h,o-z)
        dimension zw(2,lot,n2),zt(2,lzt,n1zt)

        do j=1,n1dfft
         zt(1,1,2*j-1)=zw(1,j,1)
         zt(2,1,2*j-1)=0.0d0
         zt(1,1,2*j)  =zw(2,j,1)
         zt(2,1,2*j)  =0.0d0
        enddo

        do i=2,max2+1
         do j=1,n1dfft
          zt(1,i,2*j-1)= (zw(1,j,i)+zw(1,j,n2+2-i))*0.5d0
          zt(2,i,2*j-1)= (zw(2,j,i)-zw(2,j,n2+2-i))*0.5d0
          zt(1,i,2*j)  = (zw(2,j,i)+zw(2,j,n2+2-i))*0.5d0
          zt(2,i,2*j)  =-(zw(1,j,i)-zw(1,j,n2+2-i))*0.5d0
         enddo
        enddo



!       Here, zero and positive frequencies
!        do 90,j=1,n1dfft
!        do 90,i=1,max2+1
!        zt(1,i,j)=zw(1,j,i)
!        zt(2,i,j)=zw(2,j,i)
!90      continue

!       Here, negative frequencies
!        if(m2>=max2+2)then
!         do 110,j=1,n1dfft
!         do 110,i=max2+2,m2
!         zt(1,i,j)=zw(1,j,i+n2-m2)
!         zt(2,i,j)=zw(2,j,i+n2-m2)
!110      continue
!        endif

        return
        end