Block-Projections Algorithms with Blocks Containing Mutually Orthogonal Rows and Columns

In this paper we present an algorithm which, for a given (sparse) matrix, constructs a partition of its set of row-indices, such that each subset of this partition (except the last one obtained) contains indices which correspond to mutually orthogonal rows. We then use such decompositions in some classes of block-projections methods, previously extended by the author to general inconsistent linear least-squares problems. Numerical experiments on an inconsistent and rank-deficient least-squares model problem are described in the last section of the paper.

[1]  Å. Björck,et al.  Accelerated projection methods for computing pseudoinverse solutions of systems of linear equations , 1979 .

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

[3]  A. Lent,et al.  Iterative algorithms for large partitioned linear systems, with applications to image reconstruction , 1981 .

[4]  Y. Censor,et al.  Strong underrelaxation in Kaczmarz's method for inconsistent systems , 1983 .

[5]  J. J. Moré,et al.  Estimation of sparse jacobian matrices and graph coloring problems , 1983 .

[6]  Thomas F. Coleman,et al.  Large Sparse Numerical Optimization , 1984, Lecture Notes in Computer Science.

[7]  Trond Steihaug,et al.  On the successive projections approach to least-squares problems , 1986 .

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

[9]  Y. Censor,et al.  Block-iterative projection methods for parallel computation of solutions to convex feasibility problems , 1989 .

[10]  Constantin Popa,et al.  Least-squares solution of overdetermined inconsistent linear systems using kaczmarz's relaxation , 1995, Int. J. Comput. Math..

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

[12]  Åke Björck,et al.  Numerical methods for least square problems , 1996 .

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

[14]  C. Popa.  Characterization of the solutions set of least-squares problems by an extension of Kaczmarz's projec , 1997 .