标题： 梯度型子空间迭代方法用于对称特征值问题 显示英文标题

标题： Gradient-type subspace iteration methods for the symmetric eigenvalue problem

摘要： 本文探讨了子空间迭代算法的变体，用于计算近似不变子空间。 标准的子空间迭代方法被重新审视，并开发了利用梯度类型技术结合Grassmann流形观点的新变体。 描述了一种梯度方法以及一种非线性共轭梯度技术。 分析了基于梯度的算法的收敛性，并报告了一些数值实验，表明所提出的算法有时优于标准算法。 这包括基于切比雪夫的子空间迭代和局部最优块共轭梯度方法，在矩阵向量乘积数量和计算时间方面进行比较时，效果更好。 另一方面，新方法不需要估计最优参数。 本文为实现这种良好性能的重要贡献是精确且高效的精确线搜索实现。 此外，还为非加速梯度方法提供了新的收敛证明，如果在主导子空间的$\mathcal{O(\sqrt{\delta})}$邻域内开始，并且具有谱间隙$\delta$，则包括局部指数收敛。 显示英文摘要

摘要： This paper explores variants of the subspace iteration algorithm for computing approximate invariant subspaces. The standard subspace iteration approach is revisited and new variants that exploit gradient-type techniques combined with a Grassmann manifold viewpoint are developed. A gradient method as well as a nonlinear conjugate gradient technique are described. Convergence of the gradient-based algorithm is analyzed and a few numerical experiments are reported, indicating that the proposed algorithms are sometimes superior to standard algorithms. This includes the Chebyshev-based subspace iteration and the locally optimal block conjugate gradient method, when compared in terms of number of matrix vector products and computational time, resp. The new methods, on the other hand, do not require estimating optimal parameters. An important contribution of this paper to achieve this good performance is the accurate and efficient implementation of an exact line search. In addition, new convergence proofs are presented for the non-accelerated gradient method that includes a locally exponential convergence if started in a $\mathcal{O(\sqrt{\delta})}$ neighbourhood of the dominant subspace with spectral gap $\delta$.