initf.f90

基于linux操作下的fortran快速付立变换的程序

	subroutine initf(idat,m1,m2,m3,md1,md2,md3,iproc,nproc,zf)
	implicit real*8 (a-h,o-z)
	dimension zf(2,0:md1-1,0:md3-1,0:md2/nproc-1)
	
	do 9766,jp2=0,md2/nproc-1
	j2=iproc*(md2/nproc)+jp2
	if (j2.le.m2-1) then

	i2=j2

	do 9763,j3=0,m3-1
	i3=j3

	do 9763,j1=0,m1-1
	i1=j1

	zf(1,j1,j3,jp2) =sin(1.2d0*i1+2.3d0*i2+3.4d0*i3)*sqrt(.5d0*idat)
	zf(2,j1,j3,jp2) =sin(5.6d0*i1+7.8d0*i2+9.0d0*i3)

9763	continue
	endif
9766	continue
	return
	end