Computing Krylov iterates in the time of matrix multiplication
Abstract
Krylov methods rely on iterated matrix-vector products Ak uj for an n× n matrix A and vectors u1,…,um. The space spanned by all iterates Ak uj admits a particular basis -- the maximal Krylov basis -- which consists of iterates of the first vector u1, Au1, A2u1,…, until reaching linear dependency, then iterating similarly the subsequent vectors until a basis is obtained. Finding minimal polynomials and Frobenius normal forms is closely related to computing maximal Krylov bases. The fastest way to produce these bases was, until this paper, Keller-Gehrig's 1985 algorithm whose complexity bound O(nω (n)) comes from repeated squarings of A and logarithmically many Gaussian eliminations. Here ω>2 is a feasible exponent for matrix multiplication over the base field. We present an algorithm computing the maximal Krylov basis in O(nω(n)) field operations when m ∈ O(n), and even O(nω) as soon as m∈ O(n/(n)c) for some fixed real c>0. As a consequence, we show that the Frobenius normal form together with a transformation matrix can be computed deterministically in O(nω ((n))2), and therefore matrix exponentiation~Ak can be performed in the latter complexity if (k) ∈ O(nω-1-) for some fixed >0. A key idea for these improvements is to rely on fast algorithms for m× m polynomial matrices of average degree n/m, involving high-order lifting and minimal kernel bases.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.