English

Orthogonal systems for time-dependent spectral methods

Numerical Analysis 2023-02-09 v1 Numerical Analysis

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 C1(a,b)\mathrm{C}^1(a,b) weight function such that w(a)=w(b)=0w(a)=w(b)=0, we can generate an orthonormal system with a skew-symmetric differentiation matrix. Except for the case a=a=-\infty, b=+b=+\infty, only a limited number of powers of that matrix is 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 xαexx^\alpha \mathrm{e}^{-x} for x>0x>0 and α>0\alpha>0 and the ultraspherical weight function (1x2)α(1-x^2)^\alpha, x(1,1)x\in(-1,1), α>0\alpha>0, and establish their properties. Both weights share a most welcome feature of {\em separability,\/} which allows for fast computation. The quality of approximation is highly sensitive to the choice of α\alpha and we discuss how to choose optimally this parameter, depending on the number of zero boundary conditions.

Keywords

Cite

@article{arxiv.2302.04217,
  title  = {Orthogonal systems for time-dependent spectral methods},
  author = {Arieh Iserles},
  journal= {arXiv preprint arXiv:2302.04217},
  year   = {2023}
}
R2 v1 2026-06-28T08:35:17.128Z