Rational Krylov subspaces have become a reference tool in dimension reduction procedures for several application problems. When data matrices are symmetric, a short-term recurrence can be used to generate an associated orthonormal basis. In the past this procedure was abandoned because it requires twice the number of linear system solves per iteration compared with the classical long-term method. We propose an implementation that allows one to obtain the rational subspace reduced matrices at lower overall computational costs than proposed in the literature by also conveniently combining the two system solves. Several applications are discussed where the short-term recurrence feature can be exploited to avoid storing the whole orthonormal basis. We illustrate the advantages of the proposed procedure with several examples.
翻译:理性的 Krylov 子空间已成为若干应用问题降低维度程序的参考工具。 当数据矩阵是对称时, 短期重现可以用来产生相关的正态基础。 过去, 这个程序被放弃, 因为它要求线性系统比传统长期方法多出一倍, 解决每个迭代问题。 我们建议实施一个实施方案, 允许人们以比文献中提议的较低的整体计算成本获得合理的子空间减少矩阵。 同时方便地将两个系统合并解决。 讨论了若干应用方案, 可以利用短期重现特征来避免储存整个异常基础。 我们用几个例子来说明拟议程序的优点。