MPI程序设计PC10.pptVIP

现代密码学理论与实践之五 并 行 计 算 中国科学技术大学计算机科学与技术系 国家高性能计算中心(合肥) 2004年12月 第三篇 并行数值算法 第八章 基本通讯操作 第九章 稠密矩阵运算 第十章 线性方程组的求解 第十一章 快速傅里叶变换  第十章 线性方程组的求解 10.1 三角形方程组的求解 10.2 三对角方程组的求解 10.3 稠密线性方程组的求解 10.4 稀疏线性方程组的求解  10.1 三角形方程组的求解 10.1.1 基本术语 10.1.2 上三角方程组的求解  基本术语 线性方程组的定义和符号 a1,1x1 + a1,2x2 + … + a1,nxn = b1 a2,1x1 + a2,1x2 + … + a2,nxn = b2 an,1x1 + an,1x2 + … + an,nxn = bn 记为 AX=b 10.1 三角形方程组的求解 10.1.1 基本术语 10.1.2 上三角方程组的求解  上三角方程组的求解 上三角方程组的回代解法并行化 (1)SISD上的回代算法 Begin (1)for i=n downto 1 do (1.1)xi=bi/aii (1.2)for j=1 to i-1 do bj=bj-ajixi aji=0 endfor endfor End 上三角方程组的求解 上三角方程组的回代解法并行化 (2)SIMD-CREW上的并行回代算法 - 划分: p个处理器行循环带状划分 - 算法 Begin for i=n downto 1 do xi=bi/aii for all Pj, where 1≤j≤p do for k=j to i-1 step p do bk=bk-akixi aki=0 endfor endfor endfor End // p(n)=n, t(n)=n 第十章 线性方程组的求解 10.1 三角形方程组的求解 10.2 三对角方程组的求解 10.3 稠密线性方程组的求解 10.4 稀疏线性方程组的求解  10.2 三对角方程组的求解 10.2.1 直接求解法 10.2.2 奇偶规约法  三对角方程组的直接求解法 Gauss消去法(难以并行化) ①消元 ②回代 注：由于三对角的 方程组的特殊性， 一次消元或一次 回代，只涉及邻 近一个方程，故 难以并行化。 10.2 三对角方程组的求解 10.2.1 直接求解法 10.2.2 奇偶规约法  三对角方程组的直接求解法 奇偶规约求解法(可并行化) 三对角方程可以写成如下形式 fixi-1+gixi+hixi+1=bi i=1~n f1=hn=0 串行算法描述 ①利用上下相邻方程消去偶序号方程中的奇下标变量: f2i-1x2i-2+g2i-1x2i-1+h2i-1x2i =b2i-1 f2ix2i-1 + g2ix2i + h2ix2i+1 =b2i