Abstract
In this paper we study the use of the Fourier, Sine and Cosine Transform for solving or preconditioning linear systems, which arise from the discretization of elliptic problems. Recently, R. Chan and T. Chan considered circulant matrices for solving such systems. Instead of using circulant matrices, which are based on the Fourier Transform, we apply the Fourier and the Sine Transform directly. It is shown that tridiagonal matrices arising from the discretization of an onedimensional elliptic PDE are connected with circulant matrices by congruence transformations with the Fourier or the Sine matrix. Therefore, we can solve such linear systems directly, using only Fast Fourier Transforms and the Sherman-Morrison-Woodbury formula. The Fast Fourier Transform is highly parallelizable, and thus such an algorithm is interesting on a parallel computer. Moreover, similar relations hold between block tridiagonal matrices and Block Toeplitz-plus-Hankel matrices of ordern 2×n 2 in the 2D case. This can be used to define in some sense natural approximations to the given matrix which lead to preconditioners for solving such linear systems.
Similar content being viewed by others
References
Axelsson, O., Barker, V. A.:Finite Element Solution of Boundary Value Problems: Theory and Computation, Academic Press, Orlando, Fl., 1983.
Birkhoff, G., Lynch, R. E.:Numerical Solution of Ellitic Problems, SIAM, Philadelphia, 1984.
Chan, T.:An optimal circulant preconditioner for Toeplitz systems, SIAM J. Sci. Stat. Comput. 9 (4), 766–771, 1988.
Chan, R. H., Chan, T. F.:Circulant preconditioners for elliptic problems, to appear in Journal of Numerical Linear Algebra with Appl.
Chan, R., Strang, G.:Toeplitz equations by conjugate gradients with circulant preconditioner, SIAM J. Sci. Stat. Comput. 10, 104–119, 1989.
Davis, P. J.:Circulant Matrices, John Wiley, New York, 1979.
Elliot, D. F., Rao, K. Ramamohan:Fast Transforms Algorithms, Analyses, Applications, Academic Press, New York, 1982.
Horn, R. A., Johnson, C. R.:Topics in Matrix Analysis, Cambridge University Press, Cambridge, 1991.
Huckle, T. K.:Circulant and skewcirculant matrices for solving Toeplitz matrix problems, in SIAM J. Matrix Anal. Appl. 13 (3), 767–777, 1992.
Huckle, T. K.:Circulant/Skewcirculant Matrices as Preconditioners for Hermitian Toeplitz Systems, Proceedings of the IMACS Conference on IterativeMethods in Linear Algebra, Brusseles, April 1991.
Huckle, T. K.:A Note on skewcirculant preconditioners for elliptic problems, in Numerical Algorithms 2 (3–4), 279–286, 1992.
Huckle, T. K.:Some aspects of circulant preconditioners, to appear in SIAM J. Sci. Stat. Comp.
Strang, G.:A proposal for Toeplitz matrix computations, Studies in Applied Mathematics 74, 171–176, 1986.
Tyrtyshnikov, E.:Optimal and super optimal circulant preconditioners, in SIAM J. Matrix Anal. Appl. 13 (2), 459–473, 1992.
Van Loan, C.:Computational Frameworks for the Fast Fourier Transform, SIAM Publications, Philadelphia, 1992.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Huckle, T. Fast transforms for tridiagonal linear equations. BIT 34, 99–112 (1994). https://doi.org/10.1007/BF01935019
Received:
Issue Date:
DOI: https://doi.org/10.1007/BF01935019