Fast & Accurate Randomized Algorithms for Linear Systems and Eigenvalue Problems

We develop a new class of algorithms for general linear systems and a wide range of eigenvalue problems. These algorithms apply fast randomized sketching to accelerate subspace projection methods.  This approach offers great flexibility in designing the basis for the approximation subspace, which can improve scalability in many computational environments. The resulting algorithms outperform the classic methods with minimal loss of accuracy. For model problems, numerical experiments show large advantages over MATLAB’s optimized routines, including a 100x speedup. 

Joint work with Joel Tropp (Caltech).