unswitchreal.f90

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

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

!       Decompose symmetric and antisymmetric parts
        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,n2eff
         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

        end