An efficient spectral method for ordinary differential equations with rational function coefficients
HTML articles powered by AMS MathViewer
- by Evangelos A. Coutsias, Thomas Hagstrom and David Torres PDF
- Math. Comp. 65 (1996), 611-635 Request permission
Abstract:
We present some relations that allow the efficient approximate inversion of linear differential operators with rational function coefficients. We employ expansions in terms of a large class of orthogonal polynomial families, including all the classical orthogonal polynomials. These families obey a simple 3-term recurrence relation for differentiation, which implies that on an appropriately restricted domain the differentiation operator has a unique banded inverse. The inverse is an integration operator for the family, and it is simply the tridiagonal coefficient matrix for the recurrence. Since in these families convolution operators (i.e., matrix representations of multiplication by a function) are banded for polynomials, we are able to obtain a banded representation for linear differential operators with rational coefficients. This leads to a method of solution of initial or boundary value problems that, besides having an operation count that scales linearly with the order of truncation $N$, is computationally well conditioned. Among the applications considered is the use of rational maps for the resolution of sharp interior layers.References
- Milton Abramowitz and Irene A. Stegun, Handbook of mathematical functions with formulas, graphs, and mathematical tables, National Bureau of Standards Applied Mathematics Series, No. 55, U. S. Government Printing Office, Washington, D.C., 1964. For sale by the Superintendent of Documents. MR 0167642
- A. Bayliss, D. Gottlieb, B. Matkowsky and M. Minkoff, An adaptive pseudo-spectral method for reaction diffusion problems, J. Comp. Phys., 81, (1989), 421-443.
- A. Bayliss and E. Turkel, Mappings and accuracy for Chebyshev pseudo-spectral computations, J. Comp. Phys., 101, (1992), 349–359.
- Christine Bernardi and Yvon Maday, Polynomial interpolation results in Sobolev spaces, J. Comput. Appl. Math. 43 (1992), no. 1-2, 53–80. Orthogonal polynomials and numerical methods. MR 1193294, DOI 10.1016/0377-0427(92)90259-Z
- Claudio Canuto, M. Yousuff Hussaini, Alfio Quarteroni, and Thomas A. Zang, Spectral methods in fluid dynamics, Springer Series in Computational Physics, Springer-Verlag, New York, 1988. MR 917480, DOI 10.1007/978-3-642-84108-8
- Cahit Arf, Untersuchungen über reinverzweigte Erweiterungen diskret bewerteter perfekter Körper, J. Reine Angew. Math. 181 (1939), 1–44 (German). MR 18, DOI 10.1515/crll.1940.181.1
- Tadasi Nakayama, On Frobeniusean algebras. I, Ann. of Math. (2) 40 (1939), 611–633. MR 16, DOI 10.2307/1968946
- E.A. Coutsias, T. Hagstrom, J.S. Hesthaven and D. Torres, Integration preconditioners for differential operators in spectral $\tau$-methods, 1995 (preprint).
- E.A. Coutsias, F.R. Hansen, T. Huld, G. Knorr and J.P. Lynov, Spectral Methods for Numerical Plasma Simulations, Phys. Scripta, 40, 270-279, 1989.
- B. Fornberg, A review of pseudospectral methods for solving partial differential equations, Acta Numerica (1994), 203-267.
- L. Fox and I. B. Parker, Chebyshev polynomials in numerical analysis, Oxford University Press, London-New York-Toronto, Ont., 1968. MR 0228149
- David Gottlieb and Steven A. Orszag, Numerical analysis of spectral methods: theory and applications, CBMS-NSF Regional Conference Series in Applied Mathematics, No. 26, Society for Industrial and Applied Mathematics, Philadelphia, Pa., 1977. MR 0520152, DOI 10.1137/1.9781611970425
- L. Greengard, Spectral integration and two-point boundary value problems, SIAM J. Numer. Anal. 28 (1991), no. 4, 1071–1080. MR 1111454, DOI 10.1137/0728057
- Harry Hochstadt, The functions of mathematical physics, 2nd ed., Dover Publications, Inc., New York, 1986. With a foreword by Wilhelm Magnus. MR 883962
- Cahit Arf, Untersuchungen über reinverzweigte Erweiterungen diskret bewerteter perfekter Körper, J. Reine Angew. Math. 181 (1939), 1–44 (German). MR 18, DOI 10.1515/crll.1940.181.1
- A. R. Collar, On the reciprocation of certain matrices, Proc. Roy. Soc. Edinburgh 59 (1939), 195–206. MR 8, DOI 10.1017/S0370164600012281
- Lloyd N. Trefethen and Manfred R. Trummer, An instability phenomenon in spectral methods, SIAM J. Numer. Anal. 24 (1987), no. 5, 1008–1023. MR 909061, DOI 10.1137/0724066
- Laurette S. Tuckerman, Transformations of matrices into banded form, J. Comput. Phys. 84 (1989), no. 2, 360–376. MR 1051527, DOI 10.1016/0021-9991(89)90238-6
- N. N. Lebedev, Special functions and their applications, Dover Publications, Inc., New York, 1972. Revised edition, translated from the Russian and edited by Richard A. Silverman; Unabridged and corrected republication. MR 0350075
Additional Information
- Evangelos A. Coutsias
- Affiliation: Department of Mathematics and Statistics, University of New Mexico, Albuquerque, New Mexico 87131
- Email: vageli@math.unm.edu
- Thomas Hagstrom
- Affiliation: Department of Mathematics and Statistics, University of New Mexico, Albuquerque, New Mexico 87131
- Email: hagstrom@math.unm.edu
- David Torres
- Affiliation: Department of Mathematics and Statistics, University of New Mexico, Albuquerque, New Mexico 87131
- Email: dtorres@math.unm.edu
- Received by editor(s): August 9, 1994
- Received by editor(s) in revised form: February 12, 1995
- Additional Notes: Part of the work of the first author was performed at Risø National Laboratory, DK–4000 Roskilde, Denmark. All authors supported in part by DOE Grant DE-FG03-92ER25128.
The work of the second author was partially supported by NSF Grants DMS-9108072, DMS-9304406 and by ICOMP, NASA Lewis Res. Ctr., Cleveland, OH, USA - © Copyright 1996 American Mathematical Society
- Journal: Math. Comp. 65 (1996), 611-635
- MSC (1991): Primary 65Q05, 65L60, 65P05, 76M25, 33A45, 33C55, 33C45
- DOI: https://doi.org/10.1090/S0025-5718-96-00704-1
- MathSciNet review: 1333309