假设我有一个线性系统
在哪里是对称正定的,是对称半正定的。 是一个大小兼容的满秩矩阵。 都是大小合适的向量。
在决定在大系统上使用 GMRES 和使用 Schurr 补码技术之间有哪些相关因素?
我猜想相关的因素是有多贵是反转,条件数,以及是否是零矩阵。
任何与文献的链接都值得赞赏。
编辑
我问这个问题是因为我需要解决以下系统。
现在,我使用 Schurr 补码方法解决了这个系统,但它需要三个迭代的 CG 解决方案。我想知道“蛮力”并使用GMRES是否会更好。
假设我有一个线性系统
在决定在大系统上使用 GMRES 和使用 Schurr 补码技术之间有哪些相关因素?
我猜想相关的因素是有多贵是反转,条件数,以及是否是零矩阵。
任何与文献的链接都值得赞赏。
编辑
我问这个问题是因为我需要解决以下系统。
现在,我使用 Schurr 补码方法解决了这个系统,但它需要三个迭代的 CG 解决方案。我想知道“蛮力”并使用GMRES是否会更好。
您通常可以使用块 LU 分解来操作这样的系统;然后你提到的因素开始发挥作用。然后可以通过利用分解所揭示的块对角结构来设计预条件子,其中块是和舒尔补.
这种方法的一个优点是您还可以利用近似无矩阵逆和 Schur 补码,如果您的矩阵太大而无法形成和操作,这很有用。例如,对于不可压缩流的有限元离散化,可以用几个多重网格迭代来代替。Schur 补码的预条件子也可以使用无矩阵方法来近似。
在构造这些预条件子时应该注意一些(例如,Wathen展示了用迭代方法替换矩阵逆失败的情况,因为它引入了非线性算子而不是线性算子)。
编辑:正如 Wolfgang Bangerth 指出的,这个问题也可以通过使用灵活的 CG 或 GMRES 外部迭代来解决。