vgl.f90

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

        subroutine vgl(idat,scal,m1,m2,m3,md1,md2,md3,iproc,nproc,zf,sum)
        implicit real*8 (a-h,o-z)
        dimension zf(2,0:md1-1,0:md3-1,0:md2/nproc-1)

        sum=0.d0
        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

        t1=sin(1.2d0*i1+2.3d0*i2+3.4d0*i3)*sqrt(.5d0*idat)
        t2=sin(5.6d0*i1+7.8d0*i2+9.0d0*i3)
        diff=(scal*zf(1,j1,j3,jp2)-t1)**2+(scal*zf(2,j1,j3,jp2)-t2)**2
        sum=sum+diff
        if (diff.gt.1.d-10) write(6,'(3(i3),6(x,e12.5))')  &
           j1,j3,jp2,scal*zf(1,j1,j3,jp2),t1,scal*zf(2,j1,j3,jp2),t2
9763    continue
        endif
9766    continue
        return
        end