Abstract
In this report we consider block-tridiagonal systems with Toeplitz blocks. Each block is of sizen×n consisting ofn c×n c matrices as entries, and there arem×m blocks in the system. The solution of those systems consists of 2n c m modified sine transforms and an intermediate solution ofn block-tridiagonal systems. Symmetries in the data vectors are exploited such that one modified sine transform can be computed in terms of one Fourier transform of half the length of the original one, hence requiringO(2.5nlog2 n) operations. Similarly, we only have to solve (n+1)/2 of the intermediate systems due to symmetry.
Similar content being viewed by others
References
D.H. Bailey and P.N. Swarztrauber, The fractional Fourier transform and applications, SIAM Rev. 33 (1991) 389–404.
J. Guerra and B. Gustafsson, A semi-implicit method for hyperbolic problems with different timescales. SIAM J. Numer. Anal. 23 (1986) 734–749.
B. Gustafsson and H. Stoor, Navier-Stokes equations for almost incompressible flow, SIAM J. Number. Anal. 28 (1991) 1523–1547.
L. Hemmingsson, Toepliz preconditions with block structure for first-order PDEs, Report No. 156, Dept. of Scientific Computing, Uppsala University, Uppsala, Sweden (1993).
C. van Loan,Computational Frameworks for the Fast Fourier Transform (SIAM, Philadelphia, 1992).
W. Press, B. Flannery, S. Teukolsky and W. Vetterling,Numerical Recipes (Cambridge University Press, Cambridge, 1986).
M. Pickering,An Introduction to Fast Fourier Transform Methods for Partial Differential Equations, with Applications (Research Studies Press, Letchworth, Hertfordshire, 1986).
G.D. Smith,Numerical Solution of Partial Differential Equations, Finite Difference Methods (Clarendon Press, Oxford, 1978).
P.N. Swarztrauber, Symmetric FFTs, Math. Comp. 47 (1986) 323–346.
Author information
Authors and Affiliations
Additional information
Communicated by Å. Björck
This work was supported by the Swedish National Board for Industrial and Technical Development, NUTEK, under contract No. 89-02539 P.
Rights and permissions
About this article
Cite this article
Hemmingsson, L. A fast modified sine transform for solving block-tridiagonal systems with Toeplitz blocks. Numer Algor 7, 375–389 (1994). https://doi.org/10.1007/BF02140691
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF02140691