Sine transform based preconditioners for symmetric Toeplitz systems

Research output: Journal Publications and Reviews (RGC: 21, 22, 62)21_Publication in refereed journalpeer-review

34 Scopus Citations
View graph of relations

Author(s)

Detail(s)

Original languageEnglish
Pages (from-to)237-259
Journal / PublicationLinear Algebra and Its Applications
Volume232
Issue number1-3
Publication statusPublished - 1 Jan 1996
Externally publishedYes

Abstract

The optimal circulant preconditioner for a given matrix A is defined to be the minimizer of ∥C - AF  over the set of all circulant matrices C. Here ∥·∥ is the Frobenius norm. Optimal circulant preconditioners have been proved to be good preconditioners in solving Toeplitz systems with the preconditioned conjugate gradient method. In this paper, we construct an optimal sine transform based preconditioner which is defined to be the minimizer of ∥B - A over the set of matrices B that can be diagonalized by sine transforms. We will prove that for general n-by-n matrices A, these optimal preconditioners can be constructed in O (n2) real operations and in (n) real operations if A is Toeplitz. We will also show that the convergence properties of these optimal sine transform preconditioners are the same as that of the optimal circulant ones when they are employed to solve Toeplitz systems. Numerical examples are given to support our convergence analysis.

Citation Format(s)