Uniform approximation of spectra of linear second order differential operators via discrete sine transform based discretisation
Abstract
We study various discretisation schemes for the regular Sturm--Liouville operator on a bounded interval and for the Laplace operator on an arbitrary planar domain, in both cases subject to zero Dirichlet boundary conditions. The objective is to identify a discretisation scheme such that the corresponding discretised operator---a matrix of size N × N---produces N eigenvalues that approximate as closely as possible the first N eigenvalues of the corresponding operator at the continuous level. By means of numerical experiments we examine several conventional discretisation schemes, and we corroborate the known fact that the conventional discretisations fail to achieve the objective, with the failure attributable to the poor approximation of high-index eigenvalues, and, as a result, to the non-uniform spectrum approximation. In contrast, the newly proposed discrete sine transform based discretisation scheme is designed in such a way that it replicates the asymptotic behaviour of high-index eigenvalues, thereby providing the sought uniform spectrum approximation. Of equal significance is the fact that the proposed discrete sine transform based scheme can be, unlike many Fourier transform based methods, applied to non-rectangular domains.
Related
- Symmetry breaking differential operators and Discrete Series
- FASS: Discrete Fourier Transform Adapters for Speech Separation With Incremental Learning
- The Taylor Integral and a Generalization of the Discrete Fourier Transform
- Counterexamples to uniqueness for a discrete X-ray transform under quadratic decay
- Riesz transform characterization of Hardy spaces associated with Dunkl--Schrödinger operators
BibTeX
@misc{krenek2026uniform,
title = {{Uniform approximation of spectra of linear second order differential operators via discrete sine transform based discretisation}},
author = {Oliver Křenek and V\'ıt Průša and Rebecca Tozzi and Martin Vejvoda},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.25796},
eprint = {2609.25796},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.25796},
}
RIS
TY - UNPB
TI - Uniform approximation of spectra of linear second order differential operators via discrete sine transform based discretisation
AU - Oliver Křenek
AU - V\'ıt Průša
AU - Rebecca Tozzi
AU - Martin Vejvoda
PY - 2026
DA - 2026/09/22
PB - arXiv
DO - 10.48550/arxiv.2609.25796
UR - https://arxiv.org/abs/2609.25796
AB - We study various discretisation schemes for the regular Sturm--Liouville operator on a bounded interval and for the Laplace operator on an arbitrary planar domain, in both cases subject to zero Dirichlet boundary conditions. The objective is to identify a discretisation scheme such that the corresponding discretised operator---a matrix of size N × N---produces N eigenvalues that approximate as closely as possible the first N eigenvalues of the corresponding operator at the continuous level. By means of numerical experiments we examine several conventional discretisation schemes, and we corroborate the known fact that the conventional discretisations fail to achieve the objective, with the failure attributable to the poor approximation of high-index eigenvalues, and, as a result, to the non-uniform spectrum approximation. In contrast, the newly proposed discrete sine transform based discretisation scheme is designed in such a way that it replicates the asymptotic behaviour of high-index eigenvalues, thereby providing the sought uniform spectrum approximation. Of equal significance is the fact that the proposed discrete sine transform based scheme can be, unlike many Fourier transform based methods, applied to non-rectangular domains.
ER -