求解线性方程系统 - 广义最小残余方法
大多数迭代方法的收敛性取决于系数矩阵的条件数,COND(a)
.您可以使用平衡
改善条件数一个
,并自己这使得大多数迭代求解器更容易收敛。但是,使用平衡
也导致更好的质量预处理矩阵时,你随后的因子平衡矩阵B = R * P * * C
.
您可以使用矩阵重新排序功能,例如解剖
和Symrcm.
对系数矩阵进行因子分解以生成预调节器时,对系数矩阵的行和列进行排列并使非零数最小化。这样可以减少后续求解预处理线性系统所需的内存和时间。
[1] Barrett, R., M. Berry, T. F. Chan等,线性系统解的模板:迭代方法的构建块,暹罗,费城,1994年。
[2] Saad,Yousef和Martin H. Schultz,“GMRES:求解非对称线性系统的广义最小的残余算法”暹罗j .科学。Stat。第一版。1986年7月,第7卷第3期,第856-869页。