中科院计算流体力学最新讲义cfd第14讲mpi并行程序设计初步2.pptxVIP

服务器/前端机知识回顾MPI 程序的运行原理： 服务器（前端机）编译 可执行代码复制 N 份，每个节点运行一份 调用MPI库函数 得到每个节点号 my_id 根据my_id 不同，程序执行情况不同 调用MPI 库函数进行通讯a.exea.exea.exeMPI 编程的基本思想： 主从式，对等式计算节点重点：对等式程序设计Copyright by Li Xinliang对等式设计“对等式”程序设计思想a.exea.exe站在其中一个进程的角度思考a.exe如果我是其中一个进程；我应当做……完成我需要完成的任务a.exe计算节点Copyright by Li Xinliang基本的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 XinliangMPI的消息发送机制—— 两步进行MPI_Send( A, … ) 发送MPI_Recv( B, … )接收 配合使用发送 变量A接收 到变量BCopyright by Li Xinliang阻塞接收阻塞发送一、 阻塞式通信与非阻塞式通信开始开始MPI_Send( A, … )MPI_Recv( B , …)消息成功发出消息成功接收结束结束缓冲区数据可使用缓冲区可释放阻塞式发送与接收Copyright by Li Xinliang MPI_Send( ) 返回后缓冲区可释放 sum= …… call MPI_Send(sum,……) sum= ……变量可重复利用 MPI_Recv() 返回后缓冲区数据可使用Call MPI_Recv(sum1,……)Sum=sum0+sum1…… Copyright by Li Xinliang非阻塞发送非阻塞接收启动发送启动接收发送消息接收消息立即返回立即返回计算与通信重叠计算计算通信完成通信完成释放发送缓冲区引用接收数据非阻塞消息发送与接收Copyright by Li Xinliang非阻塞消息发送MPI_ISend(buf,count,datatype,dest,tag,comm,request,ierr)In buf,count,datatype,dest,tag,commOut request,ierr Request (返回的非阻塞通信对象, 整数)非阻塞消息接收MPI_IRecv(buf,count,datatype,source,tag,comm,request,ierr)In buf,count,datatype,source,tag,commOut request,ierr非阻塞通信的完成 MPI_Wait(request,status,ierr)等待消息收发完成 MPI_Test(request, flag,stutus,ierr) MPI_Waitall(const,request_array,status,ierr) 等待多个消息完成 In requestOut status， flag (logical型)Copyright by Li Xinliang非阻塞通信调用后立即返回，缓冲区不能立即使用Sum= ……计算某变量MPI_Isend(sum ….)发送该变量 sum=……不能给变量重新赋值 （发送可能尚未完成）MPI_Irecv(sum1, ……)sum=sum0+sum1数据不能立即使用 （接收可能未完成）MPI_Isend(sum, …, request, …)……Call MPI_Wait(request,status,ierr)Sum=……√MPI_Irecv(sum1, ……, request, …)……Call MPI_Wait(request,status,ierr)Sum=sum0+sum1 √Copyright by Li Xinliang利用通信与计算重叠技术提高效率例： 计算差分串行程序 real A(N,N),B(N,N),h…..Do i=1,NB(I,1)=(A(I,2)-A(I,1))/h B(I,N)=(A(I,N)-A(I,N-1))/henddoDo j=2,N-1Do i=1,NB(i,j)=(A(i,j+1)-A