![645b20116492fd378c04ab8eac046954.png](https://i-blog.csdnimg.cn/blog_migrate/c4dbfd6ef498d9e7b6a38b9cff152c03.jpeg)
其实关于迭代解法的预调节器还有很多。前面讲迭代解法讲得特别多,多到快忘记了第一节的直接解法了。我们来复习一下,前面提到LU分解的算法
function
矩阵
实际上这个顺序也并不固定,上面是从左上到右下的顺序,也可以是下面这样的从左至右顺序,等等
function
还可以从分块的角度来看
先分解左上角,再分解右下角
最终得到
主元选取
实际上这个版本的LU分解对于一般的矩阵来说数值不稳定,因为即使
来得到解。在LU分解过程中,不断对
在每一步变换前可以从
最终
注意到交换两个元素的排列阵
而
实际上需要数值稳定性只需要
填入
如果
如果
为表示方便令
前面用到了很多的
完全填满
则完全保持稀疏性
数据依赖关系图
当
从列与列之间的依赖关系来看,这个依赖关系与
![2db945ae6c3ec2110e5f002b8c2d73ce.png](https://i-blog.csdnimg.cn/blog_migrate/c048082d9badc724f2387de02835ff54.jpeg)
由于依赖关系的传递性,这个图其实是有一定信息冗余的。当图中同时存在
因此对称剪枝方法又可以用到这里:对于
但是要注意,这个依赖关系不是从原始矩阵