set.f90

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

        subroutine set(n,zmpi)
        implicit real*8 (a-h,o-z)
        dimension zmpi(n)
        do 10,i=1,n-1,2
        if (rand().gt..5d0) then
        sign=1.d0
        else
        sign=-1.d0
        endif
        zmpi(i)=1.797693134862084d+308*sign
10      zmpi(i+1)=1.797693134862084d+308
        return
        end