CFD2011-第14讲-MPI并行程序设计初步2.pptxVIP

计算流体力学讲义 第六讲 MPI并行程序设计 （2） 李新亮 lixl@imech.ac.cn ；力学所主楼219； ;服务器/前端机;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) ;MPI的消息发送机制—— 两步进行 MPI_Send( A, … ) 发送 MPI_Recv( B, … ) 接收 ;; MPI_Send( ) 返回后缓冲区可释放 sum= …… call MPI_Send(sum,……) sum= …… 变量可重复利用 MPI_Recv() 返回后缓冲区数据可使用 Call MPI_Recv(sum1,……) Sum=sum0+sum1 …… ;;非阻塞消息发送 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型) ;非阻塞通信调用后立即返回，缓冲区不能立即使用 Sum= …… 计算某变量 MPI_Isend(sum ….) 发送该变量 sum=…… 不能给变量重新赋值 （发送可能尚未完成） MPI_Irecv(sum1, ……) sum=sum0+sum1 数据不能立即使用 （接收可能未完成） MPI_Isend(sum, …, request, …) …… Call MPI_Wait(request,status,ierr) Sum=…… √ ;;并行程序—— 以两个进程并行为例 real A(N,N/2),B(N,N/2),A1(N),h …… …… If(myid .eq. 0) then call MPI_send(A(1,N/2),N,MPI_real,1,99,MPI_Comm_world,ierr) call MPI_recv(A1,N,MPI_real,1,99,MPI_Comm_World,status,ierr) Else call MPI_recv(A1,N,MPI_real,0,99,MPI_Comm_World,status,ierr) call MPI_send(A(1,1),N,MPI_real,0,99,MPI_Comm_world,ierr) endif ; If(myid .eq. 0) then Do i=1,N B(i,1)=(A(i,2)-A(i,1))/h B(i,N)=(A1(i)-A(i,N-1))/(2.*h) Enddo Else Do i=1,N B(i,1)=(A(i,2)-A1(i))/(2.*h) B(i,N)=(A(i,N)-A(i,N-1))/h Enddo endif Do j=2,N-1 Do i=1,N B(i,j)=(A(i,j+