利用MPI计算矩阵相乘的简单算法.docVIP

利用MPICH2计算矩阵相乘的简单算法 烤鱼片 cleverysm@ MPICH2是用来进行并行运算的平台，而对矩阵算法的分解应该是并行运算应用中很常见的。今天在这里就用MPICH2写一个矩阵乘法的并行计算程序来学习一下MPICH2的使用。 首先要复习一下矩阵乘法的算法，我们在线性代数里都学过。 假设矩阵A为m行，k列。m=4，k=3。B是k行，n列，k=3，n=2，计算一个矩阵与列数据，或者叫做一个向量vector的乘积，结果矩阵C应当是m行，n列。 A: 1,2,3 4,5,6 7,8,9 8,7,6 B: 5,4 3,2 1,0 A,B,C都按照行主序表示。即A(0,0)=1,A(0,1)=2, A(0,3)=3。 则 如C(1,0)=A(1,0)×B(0,0)+A(1,1)×B(1,0)+A(1,2)×B(2,0)=4×5+5×3+6×1=41 。 C: 14,8 41,26 68,44 67,46 那如何将这个算法进行分解，分配到参与计算的各进程上去呢？简单的做法是可以将矩阵A进行划分，每个计算进程分配若干行，每次与B的一列数据计算的时候分别计算其进程所分配的行与B的这一列相乘的结果。 比如有两个进程计算，分别为进程1和进程2。1分配A的0,和2行，2 分配A的1和3行，与B的0列计算的时候，1进程将只计算得到C的一列的0和2位置两个数据，而2进程则算得1和3位置的数据