bicgstabl
求解线性方程组稳定双共轭梯度(l)法
语法
描述
试图解线性方程组x
= bicgstabl (一个
,b
)A*x = b
为x
使用双共轭梯度稳定(l)法.当尝试成功时,bicgstabl
显示确认收敛的消息。如果bicgstabl
在最大迭代次数之后未能收敛或由于任何原因停止,它将显示包含相对残差的诊断消息规范(b * x) /规范(b)
以及方法停止的迭代次数。
例子
输入参数
输出参数
更多关于
提示
大多数迭代方法的收敛性取决于系数矩阵的条件数,
气孔导度(A)
.当一个
是方的,你能用吗平衡
为了改善它的条件数,并且就其本身而言,这使得大多数迭代求解器更容易收敛。然而,使用平衡
当你随后分解平衡矩阵时,也会导致更好质量的预处理矩阵B = r * p * a * c
.您可以使用矩阵重新排序函数,例如
解剖
而且symrcm
当对系数矩阵进行因式分解以生成预处理条件时,对系数矩阵的行和列进行排列并使非零的数量最小化。这可以减少随后求解预条件线性系统所需的内存和时间。
参考文献
[1]巴雷特,R., M. Berry, T. F. Chan等,线性系统解的模板:迭代方法的构建块, SIAM,费城,1994年。