并行计算常用函数库.pptVIP

实一维变换 #include rfftw.h … //专门针对实型数据的rfftw变换 { // 数据类型 fftw_real in[N], out[N]; rfftw_plan p; //计划 int k … //创建计划 ,实到复的变换,复到实则为逆变换 p = rfftw_create_plan(N, FFTW_REAL_TO_COMPLEX, FFTW_ESTIMATE); … rfftw_one(p, in, out); rfftw_destroy_plan(p); } //程序编译时需与-lrfftw –lfftw –lm相连接 FFTW 实多维变换 #include rfftw.h … //专门针对实型数据的rfftw变换 { // 数据类型 fftw_real in[M, N], out[M, N]; rfftw_plan p; //计划 int k … //创建计划 ,实到复的变换,复到实则为逆变换 p = rfftwnd_create_plan(M, N, FFTW_REAL_TO_COMPLEX, FFTW_ESTIMATE); … rfftwnd_one(p, in, out); rfftwnd_destroy_plan(p); } //程序编译时需与-lrfftwnd –lfftw –lm相连接 FFTW FFTW的多线程并行 1.头文件: fftw_threads.h或rfftw_theads.h 2.线程初始化: int fftw_threads_init(void); 3.用到的函数 fftw_threads_one(nthreads, plan, in, out); ?//一维复变换 fftwnd_threads_one(nthreads, plan, in, out); //n维复变换 rfftw_threads_one(nthreads, plan, in, out); //一维实变换 rfftwnd_threads_one(nthreads, plan, in, out); //n维实变换 以一维复变换为例 用 fftw_threads_one (nthreads, plan, in, out) 代替调用单机 fftw_one (plan, in, out) 在Unix上,使用并行复变换的程序应该与 -lfftw_threads - lfftw - lm相连接, 使用并行实变换的程序应与 -lrfftw_threads -lfftw_threads -lrfftw -lfftw –lm 相连接 FFTW FFTW的MPI并行 调用头文件fftw_mpi.h 创建计划: fftw_mpi_plan fftw_mpi_create_plan(MPI_Comm comm, int n FFTW_FORWARD, FFTW_ESTIMATE ); 用完后通过fftw_mpi_destroy_plan(plan)来取消. 3. 变换函数为: void fftw_mpi(fftw_mpi_plan p, int n_fields, fftw_complex *local_data, fftw_complex *work); 返回时, local_data包含局限于当前进程的输出部分, 可调用: void fftw_mpi_local_sizes(fftw_mpi_plan p, int *local_n, int *local_start,int *local_n_after_transform, int *local_start_after_transform,int *total_local_size); 在Unix上, FFTW MPI的程序应与MPI库和-lfftw_mpi -lfftw -lm连接. FFTW 其他软件介绍 Jaguar（rank 1） petsc pspline Scalapack sprng sundials superlu superlu_dist Trilinos umfpack arpack aztec fftpack fftw gsl hypre libsci metis Parmetis mumps 其他软件介绍 JUGENE(rank 4)、JuRoPa(rank 13)、JUMP NAG Parallel Library、ScaLAPACK、ARPACK、PARPACK、 PETSc、MUMPS、SPRNG、ParMETIS、hy