Some projection-based direct solvers for general linear systems of equations

The Direct Projection-based solvers are direct methods for solving linear systems of equations in which an initial set of vectors are projected onto the hyperplanes of the system by using projections parallel with some specific directions also constructed during the development of the algorithm. They have been initially designed for square nonsingular systems, but further developments were produced also for more general ones (see an overview of these algorithms in [1], [6] and references therein). The main scope of our paper is to introduce and theoretically analyse a new class of such kind of direct solvers. Starting from some preliminary methods and results proposed by one of the authors in [15] we extend the theoretical analysis for one of them and design new block row and column projection versions. In the last section of the paper we use these algorithms and compare them with some classical direct projection-based ones for some numerical experiments on linear systems arising from rigid body dynamics problems.

[1]  Claude Brezinski,et al.  Projection methods for systems of equations , 1997 .

[2]  E. W. Purcell,et al.  The Vector Method of Solving Simultaneous Linear Equations , 1953 .

[3]  J. Navarro-Pedreño Numerical Methods for Least Squares Problems , 1996 .

[4]  C. Popa Block-Projections Algorithms with Blocks Containing Mutually Orthogonal Rows and Columns , 1999 .

[5]  Louis A. Hageman,et al.  Iterative Solution of Large Linear Systems. , 1971 .

[6]  T. Elfving Block-iterative methods for consistent and inconsistent linear equations , 1980 .

[7]  Gene H. Golub,et al.  Matrix computations , 1983 .

[8]  Andrzej Stachurski,et al.  Parallel Optimization: Theory, Algorithms and Applications , 2000, Parallel Distributed Comput. Pract..

[9]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

[10]  K. Tanabe Projection method for solving a singular system of linear equations and its applications , 1971 .

[11]  C. Popa DIRECT AND ITERATIVE KACZMARZ-LIKE SOLVERS , 2006 .

[12]  Michele Benzi,et al.  A Direct Projection Method for Sparse Linear Systems , 1995, SIAM J. Sci. Comput..

[13]  E. Spedicato,et al.  A class of direct methods for linear systems , 1984 .

[14]  Fridrich Sloboda A parallel projection method for linear algebraic systems , 1978 .

[15]  Michele Benzi,et al.  A Sparse Approximate Inverse Preconditioner for the Conjugate Gradient Method , 1996, SIAM J. Sci. Comput..

[16]  C. Popa Extensions of block-projections methods with relaxation parameters to inconsistent and rank-deficient least-squares problems , 1998 .