fill_corn.f90

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

        subroutine fill_corn(md1,md3,lot,n1dfft,n3,zf,zw)
	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 90,i3=1,n3/4
	do i1=1,n1dfft
	zw(1,i1,i3)=0.d0
	zw(2,i1,i3)=0.d0
	enddo
90      continue
        do 100,i3=n3/4+1,3*n3/4
	do i1=1,n1dfft
	zw(1,i1,i3)=zf(1,i1,i3-n3/4)
	zw(2,i1,i3)=zf(2,i1,i3-n3/4)
	enddo
100     continue
        do 110,i3=3*n3/4+1,n3
	do i1=1,n1dfft
	zw(1,i1,i3)=0.d0
	zw(2,i1,i3)=0.d0
	enddo
110     continue

	return
	end