We formulate a block-iterative algorithmic scheme for the solution of systems of linear inequalities and/or equations and analyze its convergence. This study provides as special cases proofs of convergence of (i) the recently proposed Component Averaging (CAV) method of Censor, Gordon and Gordon ({\it Parallel Computing}, 27:777–808, 2001), (ii) the recently proposed Block-Iterative CAV (BICAV) method of the same authors ({\it IEEE Transactions on Medical Imaging}, 20:1050–1060, 2001) and (iii) the Simultaneous Algebraic Reconstruction Technique (SART) of Andersen and Kak ({\it Ultrasonic Imaging}, 6:81–94, 1984) and generalizes them to linear inequalities. The first two algorithms are projection algorithms which use certain generalized oblique projections and diagonal weighting matrices which reflect the sparsity of the underlying matrix of the linear system. The previously reported experimental acceleration of the initial behavior of CAV and BICAV is thus complemented here by a mathematical study of the convergence of the algorithms.
Citation
SIAM Journal on Matrix Analysis and Applications, Vol. 24 (2002), pp. 40-58.