中科院计算流体力学最新讲义CFD2011-第14讲-MPI并行程序设计初步2.pptVIP

回顾 基本的MPI函数（6个） MPI初始化 MPI_Init(ierr) ； MPI结束 MPI_Finalize(ierr) 得到当前进程标识 MPI_Comm_rank(MPI_COMM_WORLD,myid,ierr) 得到通信域包含的进程数 MPI_Comm_size(MPI_COMM_WORLD,numprocs,ierr) 消息发送 MPI_Send(buf,count,datatype,dest,tag,comm, ierr) 消息接收 MPI_Recv(buf,count,datatype,source,tag,comm,status,ierr) * Copyright by Li Xinliang 非阻塞消息发送 MPI_ISend(buf,count,datatype,dest,tag,comm,request,ierr) In buf,count,datatype,dest,tag,comm Out request,ierr Request (返回的非阻塞通信对象, 整数) 非阻塞消息接收 MPI_IRecv(buf,count,datatype,source,tag,comm,request,ierr) In buf,count,datatype,source,tag,comm Out request,ierr 非阻塞通信的完成 MPI_Wait(request,status,ierr) 等待消息收发完成 MPI_Test(request, flag,stutus,ierr) MPI_Waitall(const,request_array,status,ierr) 等待多个消息完成 In request Out status， flag (logical型) * Copyright by Li Xinliang 发送非连续数据——构建新的数据结构 MPI_TYPE_VECTOR(count,blocklength,stride,oldtype,newtype,ierr) Count: 块的数量； blocklength: 每块的元素个数 Stride: 跨度 （各块起始元素之间的距离） Oldtype: 旧数据类型， Newtype: 新数据类型 （整数） 例：integer MY_TYPE Call MPI_TYPE_VECTOR(50,1,100,MPI_REAL,MY_TYPE,ierr) Call MPI_TYPE_Commit(MY_TYPE,ierr) A(1,1), A(2,1), ……… A(1,2), A(2,2) ……. A(1,3)…. * Copyright by Li Xinliang 通讯域的分割 MPI_Comm_Split(comm,color, key,New_Comm ) 0 2 1 4 3 5 7 6 8 9 10 11 Color 相同的进程在同一组 根据key的大小排序 例如： 12个进程， 分成 3行4列 Line=mod(myid,3); raw=myid/3 MPI_Comm_Split(MPI_Comm_World, raw, 0,Comm_Raw) MPI_Comm_Split(MPI_Comm_World,line,Comm_Line) Call MPI_Comm_rank(Comm_Raw,myid_raw,ierr) Call MPI_Comm_rank(Comm_line, myid_line,ierr) MPI_Comm_World * Copyright by Li Xinliang 实例 1. 用（拟）谱方法求解二维不可压N-S方程 2p 物理模型 周期性边界条件 按照给定能谱布置初始流动 研究流动的演化规律 * Copyright by Li Xinliang Fourier 变换 (1D) Fourier 变换 的特点： 求导数 - 乘积 困难： 非线性项 卷积 计算量巨大 在物理空间计算 Fourier 变换的快速算法FFT * Copyright by Li Xinliang 二维 Fourier 变换 两次一维 Fourier 变换 * Copyright by Li Xinliang 求解步骤： 1） 读入初值 2） 调用FFT 得到 3） 计算