unfill_corn.f90

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

        subroutine unfill_corn(md1,md3,lot,n1dfft,n3,zw,zf)
        implicit real*8 (a-h,o-z)
        dimension zw(2,lot,n3),zf(2,md1,md3)
! WARNING: Assuming that high frequencies are in the corners
!          and that n3 is multiple of 4

        do 100,i3=n3/4+1,3*n3/4
        do i1=1,n1dfft
        zf(1,i1,i3-n3/4)=zw(1,i1,i3)
        zf(2,i1,i3-n3/4)=zw(2,i1,i3)
        enddo
100     continue

        return
        end