在 GMRES 和 schur 补码解决方案之间做出决定时,哪些因素是相关的?

计算科学 线性代数 线性求解器 格瑞斯
2021-12-05 04:59:35

假设我有一个线性系统

[M1SSTM2][XY]=[FG],
在哪里M1Rn×n是对称正定的,M2Rm×m是对称半正定的。 S是一个大小兼容的满秩矩阵。 X,Y,F,G都是大小合适的向量。

在决定在大系统上使用 GMRES 和使用 Schurr 补码技术之间有哪些相关因素?

我猜想相关的因素是有多贵M1是反转,条件数S,以及是否M2是零矩阵。

任何与文献的链接都值得赞赏。

编辑

我问这个问题是因为我需要解决以下系统。

[M1S0ΓTST0ΞT00Ξ00Γ000][x1x2x3x4]=[f1f2f3f4]

现在,我使用 Schurr 补码方法解决了这个系统,但它需要三个迭代的 CG 解决方案。我想知道“蛮力”并使用GMRES是否会更好。

1个回答

您通常可以使用块 LU 分解来操作这样的系统;然后你提到的因素开始发挥作用。然后可以通过利用分解所揭示的块对角结构来设计预条件子,其中块是M1和舒尔补M2STM11S.

这种方法的一个优点是您还可以利用近似无矩阵逆M1和 Schur 补码,如果您的矩阵太大而无法形成和操作,这很有用。例如,对于不可压缩流的有限元离散化,M11可以用几个多重网格迭代来代替。Schur 补码的预条件子也可以使用无矩阵方法来近似。

在构造这些预条件子时应该注意一些(例如,Wathen展示了用迭代方法替换矩阵逆失败的情况,因为它引入了非线性算子而不是线性算子)。

编辑:正如 Wolfgang Bangerth 指出的,这个问题也可以通过使用灵活的 CG 或 GMRES 外部迭代来解决。