Randomized Quasi-Newton Updates are Linearly Convergent Matrix Inversion Algorithms
Iterative Hessian sketch: Fast and accurate solution approximation for constrained least-squares