fnorm.f90

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

        subroutine fnorm(ndat,m1,m2,m3,md1,md2,md3,iproc,nproc,zf,sum)
        implicit real*8 (a-h,o-z)
        dimension zf(2,md1,md3,md2/nproc,ndat)

        sum=0.d0
        do 9766,idat=1,ndat
        do 9766,jp2=1,md2/nproc
        j2=iproc*(md2/nproc)+jp2
        if (j2.le.m2) then
        do 9763,j3=1,m3
        do 9763,j1=1,m1
        sum=sum+zf(1,j1,j3,jp2,idat)**2+zf(2,j1,j3,jp2,idat)**2
9763    continue
        endif
9766    continue
        return
        end