矩阵乘法的并行化实验介绍.docVIP

北京科技大学 计算机与通信工程学院 实 验 报 告 实验名称： 学生姓名： 专 业： 班 级： 学 号： 指导教师： 实验成绩：________________________________ 实验地点： 实验时间：2015年05月 一、实验目的与实验要求 1、实验目的 1对比矩阵乘法的串行和并行算法，查看运行时间，得出相应的结论 2观察并行算法不同进程数运行结果，分析得出结论； 1编写矩阵乘法的串行程序，多运行得到汇总； 2矩阵乘法的串行实验 编写串行程序，运行汇总结果。 按照正常的矩阵乘法方法编写的串行，后多次运行得到结果汇总。 乘法的并行化 3个总进程 5个总进程 7个总进程 9个进程 16个进程 （一）乘法并行并行化算法分析 并行策略：行带划分法 描述：将C=A*B的按行划分，分得其中的几行后同时进行计算最后通信将从进程的结果合并的主进程的中B 如图：矩阵行 进程2：矩阵 进程3：矩阵 进程1：矩阵 时间复杂度分析： f(n) =6+2+8+k*n+k*n+k*n+3+10+n+k*n+k*n+n+2 (k从进程分到的行数)n)=(n)； 空间复杂度分析： 存储空间不共用，n)=n; 因此O(n)=(n); 2间隔行带划分法 描述：将C=A*B的按行划分，分得其中的几行后同时进行计算最后通信将从进程的结果合并的主进程的中B 如图：矩阵行 进程2：矩阵 进程3：矩阵 进程3：矩阵 时间复杂度分析： f(n) =6+2+8+k*n+k*n+k*n+3+10+n+k*n+k*n+n+2 (k从进程分到的行数)n)=(n)； 空间复杂度分析： 存储空间不共用，n)=n; 因此T(n)=O(n); 测试环境简介： 《VS2013》 Win7旗舰版性测试结果： 结果 串行结果： 通过比较发现两者相同，因此可以确定运算结果正确。 测试结果： ms 0.009936ms 0.010825ms 0.009207ms 0.008811ms 0.009831ms 5进程 0.007795ms 0.007648ms 0.008588ms 0.008334ms 0.007686ms 0.007731ms 7进程 0.010729ms 0.011336ms 0.011606ms 0.011403ms 0.011621ms 0.014563ms 9进程 0.016547ms 0.014439ms 0.016300ms 0.015455ms 0.015187ms 0.011240ms 16进程 0.024964ms 0.027917ms 0.024721ms 0.028152ms 0.022296ms 0.021321ms 通过对比分析得出，并行算法能够有效提高矩阵乘法运算效率； ）算法分析 复杂度： n)=2+n^2+2*n+3*n+7+n+n+7; 因此T(n)=O(n^2) 空间复杂度： ）O(1) 五：结论（讨论） 实验结论 根据前面实验结果，可以证明并行算法比算法更能有效地运算矩阵乘法，可以提效率数十倍 ，根据结果可以得出并行程序并不是数越多，运行时间越，程序运行还受计算机运器数的 2、讨论 1由于矩阵乘法的运算并不是很复杂，因此运算时间很短，达到了毫秒级，因此准确度有一定的影响；2因为时间紧迫，未能进一步分析代码，只用了一种并行算法。：附录 ：#include stdio.h #include conio.h #include iostream #include time.h #define N 30 #define M 30 void juzhen_mul(int m, int n, int *p1[M], int m1, int n1, int *p2[M]) { int i, j, x = 0; int c[N][M] = { 0 }; 2 for (i = 0; im;) { for (j = 0; jm1; j++) { c[i][x] += *(p1[i] + j)**(p2[j] + x); 30*30 m1*m } printf(%d, c[i][x]); printf( ); 2*