fftw使用說明王武中科院超級計算中心emailwangwu@-資料下載頁
2025-07-15 18:38本頁面
【正文】 w_one(plan, in, out)。 fftw_destroy_plan(plan)。 // Fortran編譯環(huán)境 double plex in, out dimension in(N), out(N) integer plan call fftw_f77_create_plan(plan, N, FFTW_FORWARD, FFTW_ESTIMATE) call fftw_f77_one(plan, in, out) call fftw_f77_destroy_plan(plan) 四 FFTW 的調(diào)用算例 SCCAS 21 . 并行 一維復(fù)變換的例子 include int main(int argc, char **argv) { const int N = 128。 fftwnd_mpi_plan plan。 fftw_plex *data。 MPI_Comm m=MPI_COMM_WORLD。 MPI_Init(amp。argc, amp。argv)。 MPI_Comm_rank(m,amp。rank)。MPI_Comm_size(m,amp。np)。 plan = fftwnd_mpi_creat_plan(m, N, 1, 8)。 fftw_mpi_local_sizes(plan,amp。local_n,amp。nx_out, amp。start_x,amp。start_x_out, amp。size)。 …allocate and initialize data . . . //分配并初始化數(shù)據(jù) data= (fftw_plex*) fftw_malloc(sizeof(fftw_plex)* size)。 preal=(fftw_real *)in。srand(rank*100)。 for(i=0。i2*local_n。i++) preal[i]=(fftw_real)rand()/RAND_MAX。 fftwnd_mpi(p, 1, data, NULL, FFTW_NORMAL_ORDER)。 fftw_mpi_destroy_plan(plan)。 MPI_Finalize( )。 } 四 FFTW 的調(diào)用算例 SCCAS 22 . FFTW的 數(shù)值測試 [scwang@LB]$ locate fftw /soft/fftw215 fftw_test fftw_mpi_test [scydf@LB270108 ~/fftw/fftwf90/2ds_test]$ %Demo 表 1：串行性能（深騰 7000, x86_64, 時間 : 秒 , 浮點性能 : Mflops） 問題規(guī)模 1024x 1024 2048 x 2048 5120 x 5120 10240 x 10240 計算時間 浮點性能 1 310 695 788 670 646 024 647 917 表 2：并行效率（深騰 7000, x86_64, 4096 x 4096 ） 處理機數(shù) 1 2 4 8 16 32 時 間 效 率 1 0..2670 四 FFTW 的調(diào)用算例 計算時間： T(Forward_plan)+T(Forward_FFT)+T(Inverse_plan)+T(Inverse_FFT) SCCAS 23 1) FFTW帶有一個 GNU形式的 configure程序 ,安裝很簡單 : tar –xvfz ./configure make make install 這將建立單機形式的復(fù)和實變換庫以及測試程序 2) 對于某些系統(tǒng) , configure腳本知道好的 CFLAG. 如果你的系統(tǒng)不知道 , 可用如下命令編譯 FFTW: ./ configure CC=“Compiler”, CFLAGS=”Compiler CFLAGS ” configure 也可接受一些 FFTW說明的標(biāo)志 : enablefloat 生成一個單精度版的 FFTW enablethreads 使能夠編譯和安裝 FFTW線程庫 enablempi 使能夠編譯和安裝 FFTW的 MPI庫 disablefortran 禁止 FFTW庫中包含 Fortran可調(diào)用程序 3) 比如調(diào)用 dfftw_one或 sfftw_one 之前的安裝選項為 : ./configure prefix=/home_soft/soft/ia64/lib/Mathlib/ CC=icc MPICC=icc CFLAGS=O3 F77=ifort FFLAGS=O3 enabledouble enabletypeprefix enablempi enablefortran enableshared (double/float對應(yīng)于雙精度或單精度 ) make make install 五 . FFTW在 Unix上的安裝 SCCAS 24
公司管理相關(guān)推薦
鄂ICP備17016276號-1