计算流体力学中科院力学所第5讲-mpi图并行程序设计初步1.pptVIP

B 按节拍循环发送数据 （同“全收集”） do step=0,P-1 id_send=mod(myid+step,P) id_recv=mod(myid-step,P) call MPI_Send(B1,N*N/P,MPI_REAL,id_send, ……) call MPI_Recv(B_tmp, N*N/P, MPI_REAL, id_recv, … ) …… 计算出 C1(id_recv) …… (C1(id_recv)=A1*B_tmp) enddo * Copyright by Li Xinliang …… send to ID 0 0 1 2 3 N-1 工作数组，存放接受来的矩阵Bk 服务器/前端机 计算节点 a.exe a.exe a.exe 回顾 MPI 程序的运行原理： 服务器（前端机）编译 可执行代码复制 N 份，每个节点运行一份 调用MPI库函数 得到每个节点号 myid 根据myid 不同，程序不同 调用MPI 库函数进行通讯 MPI 编程的基本思想： 主从式，对等式 推荐采用对等式编程思想 基本的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环境及基本编程 1） 建立MPI运行环境 （有并行机账户或在微机上安装MPI环境）。 2） 编制如下基本的MPI程序 计算S=1+2+3……+1000 要求程序可以实现N个进程的并行运行且负载尽量均衡。 N可变，程序中使用MPI_Comm_Size()函数读入N。由0号进程打印计算结果。 3）在并行环境上运行，输出结果。 要求： 提交源程序及运行情况的屏幕截图 Copyright by Li Xinliang * 2. 实现矩阵相乘的并行计算 矩阵A, B 均为N*N的方阵，试计算矩阵C=AB； 使用P个进程并行计算（N可以被P整除）； 矩阵A，B及C均采用分布式存储； A, C按行分割， B按列分割存储（见本稿 47页）。 要求编写计算C矩阵的MPI程序，并进行计算。 实际计算时，矩阵A, B请采用如下值， N设为100 计算出C矩阵后，请计算 ，并由根节点打印出来。 将S值与串行程序的结果进行对比，校验程序的正确性； 使用1,2,4,10个进程进行计算，并利用MPI_Wtime( )函数计算程序的运行时间；考核加速比及计算效率。 要求： 1）提交计算程序； 2）使用1，2，4，10个进程计算，提交计算结果（S值及计算时间）、计算效率及加速比。 * * * * 3) 得到当前进程标识 call　MPI_Comm_rank(MPI_COMM_WORLD,myid,ierr) (In) Integer:: MPI_COMM_WORLD 为进程所在的通信域 (Out) Integer:: myid, ierr 4) 得到通信域包含的进程数 Call　MPI_Comm_size(MPI_COMM_WORLD,numprocs,ierr) (In) Integer:: MPI_COMM_WORLD (Out) Integer:: numprocs,ierr * Copyright by Li Xinliang 进程的ID号 （从0开始） 最重要的参数！ MPI_COMM_WORLD: MPI预定义的通信域； 可换成创建的通信域Comm