unswitch_corn.f90

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

        subroutine unswitch_corn(n1dfft,n2,lot,n1,lzt,zw,zt)
        implicit real*8 (a-h,o-z)
        dimension zw(2,lot,n2),zt(2,lzt,n1)
! WARNING: Assuming that high frequencies are in the corners
!          and that n2 is multiple of 4

! Low frequencies
        do 100,j=1,n1dfft
        do 100,i=n2/4+1,3*n2/4
        zt(1,i-n2/4,j)=zw(1,j,i)
        zt(2,i-n2/4,j)=zw(2,j,i)
100     continue
        return
        end