Stable Spectral Methods for Time-Dependent Problems and the Preservation of Structure
Arieh Iserles
DOI: https://doi.org/10.1007/s10208-024-09647-w
IF: 3.4388
2024-02-17
Foundations of Computational Mathematics
Abstract:This paper is concerned with orthonormal systems in real intervals, given with zero Dirichlet boundary conditions. More specifically, our interest is in systems with a skew-symmetric differentiation matrix (this excludes orthonormal polynomials). We consider a simple construction of such systems and pursue its ramifications. In general, given any weight function such that , we can generate an orthonormal system with a skew-symmetric differentiation matrix. Except for the case , , only few powers of that matrix are bounded and we establish a connection between properties of the weight function and boundedness. In particular, we examine in detail two weight functions: the Laguerre weight function for and and the ultraspherical weight function , , , and establish their properties. Both weights share a most welcome feature of separability, which allows for fast computation. The quality of approximation is highly sensitive to the choice of , and we discuss how to choose optimally this parameter, depending on the number of zero boundary conditions.
mathematics, applied,computer science, theory & methods