线性方程组的求解-稳定的双共轭梯度法
试图解决线性方程系统X
=biggstab(一种
那B.
)A * x =
为了X
使用Biconjugate梯度稳定方法.当尝试成功时,Bicgstab.
显示一条消息以确认融合。如果Bicgstab.
在最大迭代次数后无法收敛或由于任何原因停止,将显示包含相对残差的诊断消息范数(b-A*x)/范数(b)
以及方法停止的迭代次数。
大多数迭代方法的融合取决于系数矩阵的条件数量,气孔导度(A)
.您可以使用平衡
改善条件数一种
,它本身就使得大多数迭代求解器更容易收敛。然而,使用平衡
当您随后因素为平衡的矩阵时也会导致更好的质量预处理器矩阵b = r * p * a * c
.
您可以使用矩阵重新排序功能,例如解剖
和symrcm
对系数矩阵进行分解以生成预条件时,排列系数矩阵的行和列,并最小化非零的数量。这可以减少后续求解预处理线性系统所需的内存和时间。
[1] 巴雷特,R.,M。Berry,T.F.Chan等人。,用于线性系统解决方案的模板:用于迭代方法的构建块,暹罗,费城,1994年。
[2] Van der Vorst,H.A.,“Bi-CGSTAB:用于非对称线性系统的解决方案的Bi-CG的快速和平滑地融合变体”暹罗J. SCI。统计。计算。,1992年3月,卷。13,2,第2页,第631-644页。