【MPI并行程序】矩阵乘法


    具体思路为创建矩阵,给每个进程分配矩阵的不同部分,每个进程根据进程号的不同处理矩阵的不同部分,计算完成后将每个进程的数组发送给一个特定的进程,由该进程合成一个完整的矩阵输出出来。问题是不一定每个进程都能分配到等长的矩阵行数。

首先是创建一个矩阵,这里在创建时使用行数和列数的自定义关系作为随机数种子生成随机数

【MPI并行程序】矩阵乘法_第1张图片

    每个进程根据自己的进程号去取矩阵的不同部分进行计算,不过这里需要考虑不同进程可能取得的行数不相等的情况,所以预留了一个进程用来处理分配给其他进程后剩下的行数,但是这里的预留会有一些问题,但是不影响计算结果,先看计算代码

 【MPI并行程序】矩阵乘法_第2张图片

 

 

 

这里的sta用来判断所有的进程是否可以分配等长的行数,如果等长就会有一个进程什么都不做,如果不等长每个进程都会参与计算,当等长时,输出如下

【MPI并行程序】矩阵乘法_第3张图片

可以看到出现了stack溢出 ** stack smashing detected ***,这里应该是由于进程size-1并没有进行任何计算,但是scatter也把该进程的内容拷贝到最后的进程里面了,所以将一个什么都没有的数组赋给一个数组会出现问题,但是当每个进程的行数不等长时,最后一个进程用上了,输出则没有问题了,不会出现stack溢出 ** stack smashing detected ***

【MPI并行程序】矩阵乘法_第4张图片

你可能感兴趣的:(#MPI并行计算#)