rnorm.f90

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

        subroutine rnorm(n1,n2,n3,nd1,nd2,nd3,iproc,nproc,rho,sum)
        implicit real*8 (a-h,o-z)
        dimension rho(nd1,nd2,nd3/nproc)

        sum=0.d0
        do 9766,jp3=1,nd3/nproc
        j3=iproc*(nd3/nproc)+jp3
        if (j3.le.n3) then
        do 9763,j2=1,n2
        do 9763,j1=1,n1
        sum=sum+rho(j1,j2,jp3)
9763    continue
        endif
9766    continue
        return
        end