do　k=1,nrep

　　tx=rtc()

　　call　init(x,nmax)

　　call　dft(x,a,phi,nmax)

　　ty=rtc()

　　call　consume(x,a,phi)

　　do　j=1,nmax

　　acc　=　acc　+　a(j)　+　phi(j)

　　end　do

　　del(4)　=　del(4)　+　(ty-tx)

　　end　do

　　del(1)　=　del(4)/real(nmax,8)

　　del(2)　=　del(1)/real(nrep,8)

　　del(3)　=　cyc*del(2)

　　print　*,'acc=',acc/real(nrep,8),'　n=',nmax,

　　&　　'　r=',nrep,　'　a=',del(1),'　b=',del(2),

　　&　　'　c=',del(3),'　w=',del(4)

　　end　program

　　聲明

　　獲得的精確性可能會隨著使用處理器的模型及其配置而發(fā)生變化，同時依賴的因素還有所用匯編器的版本。因此，在不同的操作條件下您將獲得不同的使用經驗。