高斯消去法c++_高斯消去法的算法改进

高斯消去法的过程如图所示

其中括号内的数字表示对该行处理的次数，比如第三列，该列中的第一个元素没有变化，第二个元素处理了一次，第三个元素处理了两次，处理的过程为

现将这个过程写成数组形式 A=A-B*C，于是就有了下列算法：

同传统算法相比较，改进算法只需一重循环，大大提升了效率

算法验证

这个方程组的解为x=[1,2,3]

自编程序计算结果为：

PS: Fortran中的spread函数用法。假定一个二维数组A

A(1, 2:4)是一个一维数组[12 13 14],spread(A(1, 2:4),1,2)就是如下的二维数组

spread(A(2:3, 1),2,3)就是如下的二维数组

spread(A(1, 2:4),1,2)*spread(A(2:3, 1),2,3)的结果就是

该算法的瓶颈就是spread函数的效率究竟如何？当然，任何事情都有其两面性。鱼和熊掌不可兼得。

☆☆☆ 往期相关 ☆☆☆

高斯消去法解线性方程组及MATLAB实现

高斯消去法解方程组及Python实现