Instance Optimal Iterative Methods for Matrix Function Approximation
summary
Claim: Lanczos (Krylov) is instance-optimal for rational functions—for a degree-M/Q rational f with poles outside the spectrum, Lanczos output is within a provable factor of the best Krylov-subspace approximation, provided k > max(M, Q−1) and compared to degree k−Q+1 approximants. Tools: extend CG/Lanczos projection tricks, tridiagonal T eigendecomposition, and polynomial-to-rational reductions. Concrete bound: approximation factor scales like Q times the product of condition numbers of the shifted matrices (currently exponential in Q), while experiments show only about √(Q·κ) degradation; bound already beats uniform spectral-range guarantees.
AI-written summary · about
chapters (15)
- 0:00Matrix function approximation
- 1:48Applications of matrix functions
- 3:37Krylov subspace methods
- 5:48Lanczos method
- 8:38Lanczos for polynomials
- 13:13Error bound for Lanczos
- 15:40Limitations of error bounds
- 18:29Instance optimality goal
- 20:36Rational function results
- 23:31Rational functions importance
- 24:55Improving Lanczos error bounds
- 26:20Proof sketch for inverse
- 28:29Generalizing to 1/x²
- 30:58Open questions
- 32:02Conclusion & Q&A
