English

Krylov-Simplex method that minimizes the residual in $\ell_1$-norm or $\ell_\infty$-norm

Numerical Analysis 2021-01-28 v1 Numerical Analysis Optimization and Control

Abstract

The paper presents two variants of a Krylov-Simplex iterative method that combines Krylov and simplex iterations to minimize the residual r=bAxr = b-Ax. The first method minimizes r\|r\|_\infty, i.e. maximum of the absolute residuals. The second minimizes r1\|r\|_1, and finds the solution with the least absolute residuals. Both methods search for an optimal solution xkx_k in a Krylov subspace which results in a small linear programming problem. A specialized simplex algorithm solves this projected problem and finds the optimal linear combination of Krylov basis vectors to approximate the solution. The resulting simplex algorithm requires the solution of a series of small dense linear systems that only differ by rank-one updates. The QRQR factorization of these matrices is updated each iteration. We demonstrate the effectiveness of the methods with numerical experiments.

Keywords

Cite

@article{arxiv.2101.11416,
  title  = {Krylov-Simplex method that minimizes the residual in $\ell_1$-norm or $\ell_\infty$-norm},
  author = {Wim Vanroose and Jeffrey Cornelis},
  journal= {arXiv preprint arXiv:2101.11416},
  year   = {2021}
}

Comments

22 pages