🤖 AI Summary
This work addresses the challenge posed by the unknown condition number of Chebyshev-filtered vector sets, which hinders the efficient selection of QR decomposition algorithms. We propose a lightweight algorithm to compute a tight upper bound on the condition number of the filtered subspace, enabling accurate and efficient prediction. Leveraging this estimate, our method dynamically selects the optimal QR decomposition strategy. Integrated into the ChASE library for the first time, this approach establishes an adaptive QR decomposition mechanism that significantly enhances computational performance without compromising numerical accuracy. The key innovation lies in tightly coupling condition number estimation with algorithmic choice, thereby providing intelligent scheduling support for linear algebra operations within subspace iteration methods.
📝 Abstract
Chebyshev filtered subspace iteration is a well-known algorithm for the solution of (symmetric/Hermitian) algebraic eigenproblems which has been implemented in several application codes~\cite{Kronik:2006ff, abinit:2020} or in stand alone libraries~\cite{ChASE}. An essential part of the algorithm is the QR-factorization of the array of vectors spanning the active subspace that have been filtered by the Chebyshev filter. Typically such an array has an a-priori unknown high condition number that directly influences the choice of QR-factorization algorithm. In this work we show how such condition number can be bound from above with precise and inexpensive estimates. We then proceed to use these estimates to implement a mechanism for the choice of QR-factorization in the ChASE library. We show how such mechanism enhance the performance of the library without compromising on its accuracy.