A preconditioner for systems with symmetric Toeplitz blocks
Srinivasa M. Salapaka, Anthony P. Peirce, Marie Dillon Dahleh · Proceedings of the 40th IEEE Conference on Decision and Control (Cat. No.01CH37228) · 2003
Proposes and studies the performance of a preconditioner used in the preconditioned conjugate gradient method for solving a class of symmetric positive definite systems, A/sub p/x = b, which we call lower rank extracted systems (LRES). These systems correspond to integral equations with convolution kernels defined on a union of many line segments in contrast to only one line segment in the case of Toeplitz systems. The p /spl times/ p matrix, A/sub p/, is shown to be a principal submatrix of a larger N /spl times/ N Toeplitz matrix, AN. The preconditioner is provided in terms of the inverse of a 2N /spl times/ 2N circulant matrix constructed from the elements of AN. The preconditioner is shown to yield clustering in the spectrum of preconditioned matrix similar to the clustering results in iterative algorithms used to solve Toeplitz systems. The analysis further demonstrates that the computational expense to solve LRE systems is reduced to O(N log N).