高性能计算系列讲座.PDFVIP

高性能计算系列讲座： MPI 并行程序设计初步（三） —— CFD 并行编程实例及优化技术 李新亮 中国科学院力学研究所LNM lixl@ 上一讲回顾 实例（1）用拟谱方法求解不可压N-S方程 实例（2 ） 用流水线方法计算紧致差分 常用的优化方法 回顾 基本的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) • 非阻塞消息发送 • 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型) 发送非连续数据——构建新的数据结构 MPI_TYPE_VECTOR(count,blocklength,stride,oldtype,newtyp e,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) …. 通讯域的分割 MPI_Comm_Split(comm,color, key,New_Comm ) MPI_Comm_World Color 相同的进程在同一组 0 3 6 9 根据key的大小排序 例如： 12个进程，分成3行4列 1 4 7 10 2 5 8 11 Line=mod(myid,3); raw=myid/3 MPI_Comm_Split(MPI_Comm_World, raw, 0,Comm_Raw) MPI_Comm_Split(MPI_Comm_World,line,Co