Decremental $(1+ε)$-Approximate Maximum Eigenvector: Dynamic Power Method
Deeksha Adil, Thatchaphol Saranurak
TL;DR
This work develops a dynamic algorithm to maintain a $(1+\epsilon)$-approximate maximum eigenvalue and eigenvector of a PSD matrix undergoing decremental rank-one updates, achieving polylog-update efficiency under an oblivious adversary. The core method adapts the classical power method to a dynamic setting and is supported by a novel analysis that uses subspace-based potentials to bound the number of power-method executions. A key contribution is a conditional lower bound showing that an efficient adaptive-adversary algorithm with a natural property would imply breakthroughs in PSD-certification hardness, suggesting a separation between oblivious and adaptive models. The paper also connects these results to dynamic PSDPs and outlines open questions, including incremental updates and broader dynamic SDP frameworks.
Abstract
We present a dynamic algorithm for maintaining $(1+ε)$-approximate maximum eigenvector and eigenvalue of a positive semi-definite matrix $A$ undergoing \emph{decreasing} updates, i.e., updates which may only decrease eigenvalues. Given a vector $v$ updating $A\gets A-vv^{\top}$, our algorithm takes $\tilde{O}(\mathrm{nnz}(v))$ amortized update time, i.e., polylogarithmic per non-zeros in the update vector. Our technique is based on a novel analysis of the influential power method in the dynamic setting. The two previous sets of techniques have the following drawbacks (1) algebraic techniques can maintain exact solutions but their update time is at least polynomial per non-zeros, and (2) sketching techniques admit polylogarithmic update time but suffer from a crude additive approximation. Our algorithm exploits an oblivious adversary. Interestingly, we show that any algorithm with polylogarithmic update time per non-zeros that works against an adaptive adversary and satisfies an additional natural property would imply a breakthrough for checking psd-ness of matrices in $\tilde{O}(n^{2})$ time, instead of $O(n^ω)$ time.
