在我看来,有 3 种方法可以解决非结构化密集方程组:GMRES、QMR 和 LU。
有没有人对比过这三个?据我所知,LU 是首选,它在解决非结构化系统方面相当流行。我正在尝试实施 GMRES 和 QMR 作为替代。
编辑:限制问题的范围
我的问题是:
- 有没有人发表过 GMRES 和 QMR 之间的比较?如果没有,有没有人发表过关于两者融合的很好的调查?有这方面的资源吗?
在我看来,有 3 种方法可以解决非结构化密集方程组:GMRES、QMR 和 LU。
有没有人对比过这三个?据我所知,LU 是首选,它在解决非结构化系统方面相当流行。我正在尝试实施 GMRES 和 QMR 作为替代。
编辑:限制问题的范围
我的问题是:
还有更多方法可以解决非结构化密集方程组(改变一个众所周知的标题,可能至少有 19 种不同的方法)。
特别是,CGS 和 BiCGStab 是有趣的迭代,而 QR 是一个非常重要的直接迭代 - 在数值上比具有列旋转的 LU 更稳定。
比较是没有实际意义的,因为优点和缺点取决于所需的准确性、所需的稳健性水平以及(在随机测试中)假设的矩阵条目的分布。它还取决于实现(LU 有不同的旋转变体),矩阵是否适当缩放等。
Vanilla GMRES 微乎其微地收敛,如下所示步骤,它最小化了整个空间。(如果你限制记忆,那么你有很多不同的 GMRES,比较就更成问题了。)
对于特定条件下的超线性收敛,请参见 http://ta.twi.tudelft.nl/nw/users/vuik/papers/vdV93V.pdf
QMR 的定义论文“QMR:非厄米线性系统的准最小残差方法”在第 6 节中给出了收敛性分析,同样是在某些条件下。