计算物理-TPG2.ppt

计算物理 并行计算方法 简介(1/5) 简介(2/5) 简介(3/5) 简介(4/5) 简介(5/5) Hello World(1/5) Hello World(2/5) Hello World(3/5) Hello World(4/5) Hello World(5/5) 六个接口的MPI子集(1/12) 六个接口的MPI子集(2/12) 六个接口的MPI子集(3/12) 六个接口的MPI子集(4/12) 六个接口的MPI子集(5/12) 六个接口的MPI子集(6/12) 六个接口的MPI子集(7/12) 六个接口的MPI子集(8/12) 六个接口的MPI子集(9/12) 六个接口的MPI子集(10/12) 六个接口的MPI子集(11/12) 六个接口的MPI子集(12/12) 安全MPI编程(1/3) 安全MPI编程(2/3) 安全MPI编程(3/3) 对等模式的程序设计(1/5) 对等模式的程序设计(2/5) 对等模式的程序设计(3/5) 对等模式的程序设计(4/5) 对等模式的程序设计(5/5) 主从模式的程序设计(1/4) 主从模式的程序设计(2/4) 主从模式的程序设计(3/4) 主从模式的程序设计(4/4) 纯并行的对等模式(1/4) 纯并行的对等模式(2/4) 纯并行的对等模式(3/4) 纯并行的对等模式(4/4) MPI通信模式(1/5) MPI通信模式(2/5) MPI通信模式(3/5) MPI通信模式(4/5) MPI通信模式(5/5) 作业 √ √ 主从模式 主进程：分派任务，收集结果，输出 次进程：执行任务，递交结果 离散动力系统的 Feigenbaum 状态图 平方迭代方程：xn+1 = r xn (1-xn) 选择一个参数 r 初始点 x0 = 0.5，进行迭代2?10 3次 利用后面的10 3次结果，得到点序列 (r, xn) 以 r 为横轴，xn 为纵轴，画点序列 选择另一个 r，重复以上过程 串行程序 并行程序 结束 开始 迭代计算 xn+1 初始化点序列数 组PP，距离 dR IR++够大？ PP(IR,N)=xn+1 Y Y N N N++够大？ xn+1越界？ 放弃？ 输出PP Y N √ 迭代计算 xn+1 并记录在 PC 迭代计算 xn+1 并记录在 PC 开始 主进程0 次进程1 Y PP(R)=PC MPISEND PC MPISEND R 到 D MPIRECV PC MPIRECV R R R1 ？ N N Y 输出 PP MPI结束 MPI结束 结束 提取源 D 和 R 次进程2 MPISEND PC MPIRECV R R R1 ？ N Y MPI结束 MPI和系数初始化 MPI和系数初始化 MPI和系数初始化 IR++在 [NUM-1,MR+NUM-2] 范围？ idX Recv R/PC from idY : R=IR idX Send R/PC to idY : R=IR √ MR=8: IR=0?7 Mathematica 作图代码 √ Feigenbaum 状态图 √ 纯并行的对等模式 在程序设计阶段分派进程 不需要任何同步或等待 结果又其它程序汇集 离散动力系统的 Feigenbaum 状态图 在程序设计阶段，通过MYID选择不同的参数 r 每个进程输出各自的数据文件 由其它程序/软件汇集成总的结果 并行程序 结束 开始 迭代计算 xn+1 并记录在 PP MPI和系数初始化 IR够大？ Y N N++够大？ 输出PP Y N IR=MYID IR=IR+NUMPROCS √ √ √ Mathematica 作图代码 √ 优化通信的更多考虑——重叠计算与I/O 缓存发送的数据 接收调用执行后才执行发送操作 正确返回发送调用 重用发送缓冲区 两种通信机制：阻塞(消息送出-执行后续)和非阻塞 MPI的四种通信模式 MPI_RSEND 就绪通信模式 MPI_SSEND 同步通信模式 MPI_BSEND 缓存通信模式 MPI_RECV MPI_SEND 标准通信模式 接收 发送 通信模式 非标准的通信模式，只有发送操作而没有相应的接收操作 √ 标准通信模式 由MPI自身决定：是否缓存发送的数据 MPI缓冲区充足 缓存将要发出的数据 发送操作与接收操作是否执行无关 发送操作可以正确返回 进程0 B MPI_SEND MPI BUF MPI BUF A PROC BUF PROC BUF MPI BUF MPI BUF