地震软件结构4.ppt

地震处理软件系统结构与分析;*;并行计算;并行计算基础;并行计算;并行计算基础;科学研究方法;应 用 需 求;美国HPCC计划公布的重大挑战性应用课题;美国HPCC方面重大挑战性课题的需求;应用需求类型;如何解决大规模科学计算问题;优化设计层次; 计算速度如何提升？;为什麽要采用并行计算?;并行计算基础;;并行计算机的分类(续);并行计算机的发展;并行计算;并行算法的设计;并行算法的定义与分类;数值算法与非数值算法;问题规模;相关性与可并行化;数据相关;数据反相关;数据输出相关;并行计算模型;PRAM模型;PRAM模型;LogP 模型 ;LogP 模型;并行算法的一般设计过程;PCAM设计方法学;PCAM设计过程;划分;四种通讯模式;局部/全局通讯;结构化/非结构化通讯;并行计算运行模式;相并行模式;分治并行模式;流水线并行模式;主-从并行模式;工作池并行模式;并行程序的调试;航天部某所串行程序并行化;计算流程及空间离散;空间离散分块并行;并行计算流程;第二章 MPI简介;什么是MPI?;MPI的历史;MPI的实现;MPI的安装--MPICH;MPI的安装--MPICH;MPI的安装--MPICH;MPI的编译;MPI的运行;第三章 基本MPI编程;第三章 基本MPI编程;MPI 介绍;程序1、第一个FORTRAN77+MPI程序;程序1在一台机器上执行的结果：; 程序2、 简单C+MPI的例子; 程序2、 简单C+MPI的例子;MPI程序的一般结构;MPI程序结构; ;第三章 基本MPI编程;MPI函数——并行环境管理函数;MPI函数——并行环境管理函数;MPI函数——通信子操作;MPI函数——通信子操作;MPI函数——点到点通信函数;MPI函数——点到点通信函数;MPI函数——点到点通信函数;MPI函数——点到点通信函数;阻塞通信;参数说明;参数说明;参数说明;进程号(rank); ; MPI_Send(buf,count,datatype,dest,tag,comm);程序3、 环形的消息传递;程序3、 环形的消息传递;程序3、 环形的消息传递;MPI函数——点到点通信函数;MPI函数——点到点通信函数;MPI函数——点到点通信函数;MPI_Sendrecv示例;MPI函数——点到点通信函数;MPI函数——点到点通信函数;MPI_Sendrecv_replace示例;MPI函数——点到点通信函数;MPI函数——点到点通信函数;程序4、 π的计算;程序4、 π的计算—send(recv);; for (j=1;jnumprocs;j++) { MPI_Send(n,1,MPI_INT,j,99,MPI_COMM_WORLD); } } else { MPI_Recv(n,1,MPI_INT,MPI_ANY_SOURCE,99, MPI_COMM_WORLD,status); } h=1.0/(double) n; sum=0.0; for(i=myid+1;i=n;i+=numprocs) num=n/numprocs for(i=myid*num+1,i=(myid+1)*num,i++) ;;;;MPI程序的一些惯例;MPI的核心函数;时间函数;第三章 基本MPI编程; ;; 程序10 串行表示的Jacobi迭代;用MPI程序实现Jacobi迭代;;;程序11 并行的Jacobi迭代;;;;;; ;二.主从模式的MPI程序设计;发送矩阵A的各行数据 回收各行与B相乘的结果;程序12 矩阵向量乘;;;; (如果还有其他的行没有被计算，则继续发送) if (numsent.le.rows) then do j=1,cols (准备好新一行的数据) buffer(j)=a(numsent+1,j) end do (将该行数据发送出去) call MPI_SEND(buffer,cols,MPI_DOUBLE_PRECISION, * sender,numsent+1,MPI_COMM_WORLD,ierr) numsent=numsent+1;;