14 假设我有一个线性系统,对于所有使用合适的Krylov方法(例如CG或GMRES)可以快速收敛。如果是具有低秩的矩阵,系统上的相同Krylov方法也会快速收敛吗?(理想情况下,迭代的额外数目大约仅取决于)?Ax=bAx=bbbBBrr(A+B)x=b(A+B)x=brr 这样的系统的一个例子是良好的预调节膜弹性和弯曲度加上未预调节的具有密集外部产品结构的气压项。 注意,问题是有或没有预处理的相同,由于是秩的修改。P(A+B)Q=PAQ+PBQP(A+B)Q=PAQ+PBQrrPAQPAQ krylov-method — 杰弗里·欧文 source
7 如果您的Krylov子空间基于幂,则收敛最多会在更正等级处延迟许多迭代。如果它是基于A T A的幂,则此数字最多为两倍。AAATAATA 我在文献中还没有看到这样的说法。但是要查看第一种情况的有效性,足以证明矩阵A + U S V T的第个Krylov空间(其中U ,V具有r列)包含在相应的空间中,而无需进行低秩校正,但具有相应地,较高的索引k + r。这很容易验证。kkA+USVTA+USVTU,VU,Vrrk+rk+r — 阿诺德·纽迈耶 source 您能否解释“基于幂”的含义?该解算器克雷洛夫给出关于信息甲+ 乙只,而不是甲直接。AAA+BA+BAA — 杰弗里·欧文 没关系:想必您是指所讨论矩阵的幂,因此在这种情况下为。A+BA+B — 杰弗里·欧文 是。该方法具有一个矩阵作为参数,该矩阵通常由表示。AA — 阿诺德·诺伊迈耶 也许对于进一步的兴趣,你可以与一些要求重写等式(或解决方案),以X = ( Ë + Σ ∞ ķ = 1(一个- 1个乙) ķ )一- 1 b可能的帮助,如果乙幂零或A − 1 B小规范。人们还认识到对解决问题的依赖。BBx=(E+∑∞k=1(A−1B)k)A−1bx=(E+∑k=1∞(A−1B)k)A−1bBBA−1BA−1Bundisturbed — 巴斯蒂安·埃伯林