附錄 3 – 驅(qū)動器程序

　　下面是一個用于計時 DFT 代碼的 Fortran 驅(qū)動器程序。

　　program　main

　　interface

　　subroutine　dft(x,a,phi,n)

　　real*8　x(n),a(n),phi(n)

　　integer　n

　　end　subroutine

　　subroutine　init(a,n)

　　real*8　a(n)

　　integer　n

　　end　subroutine

　　subroutine　consume(a,b,c)

　　real*8　a(*),b(*),c(*)

　　end　subroutine

　　end　interface

　　!　Parameters:

　　!　　nmax　is　the　problem　size.

　　!　　nrep　is　the　number　of　repetitions　of　the

　　!　　　problem.　This　should　be　chosen　so　that

　　!　　　the　elapsed　time　is　long　enough　to　give

　　!　　　sufficient　timing　resolution.

　　!　　cyc　is　the　clock　frequency　in　Hz　for　the

　　!　　　processor　that　the　program　is　to　be　run　on.

　　!　　　(Can　be　found　from　AIX　command　pmcycles.)

　　integer,　parameter　::　nmax=1000

　　integer,　parameter　::　nrep=100

　　real*8,　parameter　::　cyc=4704000000.d0

　　real*8　x(nmax),　a(nmax),　phi(nmax)

　　real*8　tx,　ty,　accum,　del(4)

　　intrinsic　sin,　sqrt

　　real*8　rtc

　　del(4)=0.d0

　　acc　=　0.d0