用mpi进行矩阵三角分解

这个问题可够麻烦的啊。总体思路：
主要利用主行i对其余各行j（j>i)做初等行变换。首先对矩阵进行行交叉划分；设处理器个数p，矩阵的阶数n，m={n/p}，编号为i的处理器存有矩阵的第i，i+p，。。。。。，i+(m-1)p行。以第0，1，......，n-1，为主行，广播到所有处理器。个处理器再用主行对其部分行向量做行变换。这实际上是个处理器轮流选出主行并广播。若以编号为my-rank的处理器的第i行元素做为主行，对其第i+1，....，m-1行数据做行变换。其他处理器利用主行元素对其第i,。。。。,m-1行数据做行变换。这样基本上就可以了说实话你一分不给，我也够忙活一顿的。