A Geometric Algorithm for Solving Linear Systems
Abstract
Based on the geometric {\it Triangle Algorithm} for testing membership of a point in a convex set, we present a novel iterative algorithm for testing the solvability of a real linear system , where is an matrix of arbitrary rank. Let be the ellipsoid determined as the image of the Euclidean ball of radius under the linear map . The basic procedure in our algorithm computes a point in that is either within distance to , or acts as a certificate proving . Each iteration takes operations and when is well-situated in , the number of iterations is proportional to . If is solvable the algorithm computes an approximate solution or the minimum-norm solution. Otherwise, it computes a certificate to unsolvability, or the minimum-norm least-squares solution. It is also applicable to complex input. In a computational comparison with the state-of-the-art algorithm BiCGSTAB ({\it Bi-conjugate gradient method stabilized}), the Triangle Algorithm is very competitive. In fact, when the iterates of BiCGSTAB do not converge, our algorithm can verify is unsolvable and approximate the minimum-norm least-squares solution. The Triangle Algorithm is robust, simple to implement, and requires no preconditioner, making it attractive to practitioners, as well as researchers and educators.
Cite
@article{arxiv.2004.12978,
title = {A Geometric Algorithm for Solving Linear Systems},
author = {Bahman Kalantari and Chun Lau and Yikai Zhang},
journal= {arXiv preprint arXiv:2004.12978},
year = {2020}
}
Comments
12 pages, 6 figures