Randomized flexible Krylov methods for $\ell_p$ regularization
By: Malena Sabaté Landman, Yuji Nakatsukasa
Potential Business Impact:
Speeds up solving hard math problems for computers.
The computation of sparse solutions of large-scale linear discrete ill-posed problems remains a computationally demanding task. A powerful framework in this context is the use of iteratively reweighted schemes, which are based on constructing a sequence of quadratic tangent majorants of the $\ell_2$-$\ell_1$ regularization functional (with additional smoothing to ensure differentiability at the origin), and solving them successively. Recently, flexible Krylov-Tikhonov methods have been used to partially solve each problem in the sequence efficiently. However, in order to guarantee convergence, the complexity of the algorithm at each iteration increases with respect to more traditional methods. We propose a randomized flexible Krylov method to alleviate the increase of complexity, which leverages the adaptability of the flexible Krylov subspaces with the efficiency of `sketch-and-solve' methods. A possible caveat of the mentioned methods is their memory requirements. In this case, one needs to rely instead on inner-outer schemes. In these scenarios, we propose a `sketch-to-precondition' method to speed up the convergence of each of the subproblems in the sequence. The performance of these algorithms is shown through a variety of numerical examples.
Similar Papers
Flexible inner-product free Krylov methods for inverse problems
Numerical Analysis
Makes computer math faster and use less memory.
Randomized Krylov methods for inverse problems
Numerical Analysis
Cleans up blurry pictures and earthquake maps.
New flexible and inexact Golub-Kahan algorithms for inverse problems
Numerical Analysis
Improves blurry pictures and scans using math.