矩阵乘法的并行计算及可扩展性分析.pdfVIP

( ) 第 21 卷第 2 期 重庆工商大学学报 自然科学版 2004 年 4 月 Vol . 21 　No . 2 J Chongqing Technol Business Univ. (Nat Sci Ed) Apr . 2004 ( ) 　　文章编号 :1672 - 058X 2004 02 - 0121 - 03 矩阵乘法的并行计算及可扩展性分析 雷　澜1 ,2 (1. 重庆大学 理学院 ,重庆 400044 ; 2 . 重庆工商大学 理学院 ,重庆 400067) 摘 　要 :针对科学与工程计算中的大型矩阵乘积 ,与原有的串行算法相比照 ,利用向量外积法 , 构造了一种适用于多处理机系统的并行算法 ,给出了算法及实现过程 ,并进一步分析了算法的可扩 展性 . 关键词 :并行算法 ;矩阵乘积 ;分块 ;可扩展性 中图分类号 :O 241. 6 　　　　文献标识码 :A 传统的计算机在计算矩阵乘法时采用串行算法 ,这需要占用较多的工作单元和较大的存储单元 . 矩阵 3 乘法的时间复杂度按一般算法为 O ( n ) , 其中 n 为矩阵的阶 . 而随着 n 值的增大, 计算效率将受到很大的 影响 . 随着科学技术的发展, 并行计算机的出现, 使大规模的提高矩阵的运算速度成为可能, 特别是 n 值很 大的矩阵, 采用并行算法, 已经可以实现 . 并行计算机利用多个处理器同时解决某一问题, 是当前提高计算 速度的主要途径 . 矩阵的向量积运算是数值分析领域中十分常用的基本运算 . 1 　矩阵乘法的内积算法 设 A = ( aij) 是 m ×p 矩阵, B = ( bij) 是 p ×n 矩阵, 其乘积 C = AB = ( Cij ) 为 m ×n 矩阵 . 矩阵乘法的 内积算法就是普通的矩阵乘法, 其中 p Cij = aik bkj 　　i = 1,2 , . . . , m ; j = 1,2 , . . . n ; k = 1 其算法如下 : SIMD 机上的矩阵相乘算法 : 输入 :A = ( aij) , B = ( bij) 输出: C m ×n begin 　　for i = 1 to m do 　　　for j = 1 to n do ( ) 　　　　 1 cij ←0 ( ) 　　　　 2 for s = 1 to p do