Many scientific applications require the evaluation of the action of the matrix function over a vector and the most common methods for this task are those based on the Krylov subspace. Since the orthogonalization cost and memory requirement can quickly become overwhelming as the basis grows, the Krylov method is often restarted after a few iterations. This paper proposes a new acceleration technique for restarted Krylov methods based on randomization. The numerical experiments show that the randomized method greatly outperforms the classical approach with the same level of accuracy. In fact, randomization can actually improve the convergence rate of restarted methods in some cases. The paper also compares the performance and stability of the randomized methods proposed so far for solving very large finite element problems, complementing the numerical analyses from previous studies.
翻译:众多科学应用需要计算矩阵函数作用于向量的结果,基于Krylov子空间的方法是完成该任务最常用的技术。由于正交化成本与内存需求会随着基向量的增长而急剧增加,Krylov方法通常在若干次迭代后需要重启。本文提出一种基于随机化的重启Krylov方法加速技术。数值实验表明,在保持相同精度水平的前提下,随机化方法显著优于经典方法。事实上,在某些情况下随机化甚至能够提升重启方法的收敛速率。本文还比较了目前已提出的随机化方法在求解超大规模有限元问题时的性能与稳定性,对先前研究的数值分析形成了补充。