pith. sign in

arxiv: 1907.07615 · v1 · pith:MMFW6QVEnew · submitted 2019-07-12 · 🧮 math.NA · cs.NA

Weakly regular Sturm-Liouville problems: a corrected spectral matrix method

Cecilia Magherini This is my paper

Pith reviewed 2026-05-24 22:09 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Sturm-Liouville problemsspectral methodsGalerkin methodeigenvalue correctionconvergence analysisweakly regular problemsunbounded potentialnumerical eigenvalues
0
0 comments X

The pith

A Galerkin spectral matrix method for weakly regular Sturm-Liouville problems with unbounded potentials at both ends yields a low-cost correction formula for numerical eigenvalues.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines Sturm-Liouville eigenvalue problems that are only weakly regular because the potential grows unbounded at each endpoint. It applies a Galerkin spectral matrix discretization and derives error estimates for the resulting eigenvalue approximations. These estimates supply an explicit, inexpensive formula that adjusts the raw numerical values to reduce their error. Numerical experiments on several test cases confirm that the corrected values match the true eigenvalues more closely. Readers care because such singular problems appear in physical models, and the correction avoids the expense of refining the discretization further.

Core claim

For weakly regular Sturm-Liouville eigenproblems with unbounded potential at both endpoints, the convergence analysis of the Galerkin spectral matrix method supplies a low-cost and effective formula for correcting the computed eigenvalues, and numerical experiments confirm that the corrected values are substantially more accurate.

What carries the argument

The correction formula for numerical eigenvalues obtained from the convergence analysis of the Galerkin spectral matrix method.

If this is right

  • The approach computes accurate eigenvalues for problems whose potentials are unbounded at both endpoints.
  • The correction requires only low additional cost once the raw eigenvalues are known.
  • The same convergence analysis that justifies the discretization also justifies the correction step.
  • Numerical experiments on multiple cases show the corrected eigenvalues are reliable.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The correction technique might be applied to other spectral discretizations of singular Sturm-Liouville problems without major changes.
  • Similar error-driven corrections could improve efficiency in related eigenvalue computations that involve endpoint singularities.
  • The method could allow practitioners to use coarser spectral resolutions while still reaching a target accuracy level.

Load-bearing premise

The convergence analysis developed for the Galerkin spectral matrix method applies directly to weakly regular Sturm-Liouville problems with unbounded potential at both endpoints and produces a usable correction formula.

What would settle it

If the correction formula applied to the method's output eigenvalues fails to reduce the error relative to the uncorrected values on standard test problems with singularities at both ends, the central claim would be falsified.

Figures

Figures reproduced from arXiv: 1907.07615 by Cecilia Magherini.

Figure 1
Figure 1. Figure 1: Relative errors for problems with potential (46). [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Relative errors for the symmetric problems with potential [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Relative errors for problems with potentials (50) subject [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
read the original abstract

In this paper, we consider weakly regular Sturm-Liouville eigenproblems with unbounded potential at both endpoints of the domain. We propose a Galerkin spectral matrix method for its solution and we study the error in the eigenvalue approximations it provides. The result of the convergence analysis is then used to derive a low-cost and very effective formula for the computation of corrected numerical eigenvalues. Finally, we present and discuss the results of several numerical experiments which confirm the validity of the approach.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript develops a Galerkin spectral matrix method for weakly regular Sturm-Liouville eigenproblems featuring unbounded potentials at both endpoints. It conducts a convergence analysis of the eigenvalue approximations produced by the method, derives a low-cost correction formula directly from the error analysis, and validates the overall approach through several numerical experiments.

Significance. If the convergence analysis and the subsequent derivation of the correction formula hold, the work supplies a practical, low-cost enhancement to spectral approximations for a difficult class of singular eigenvalue problems. The explicit linkage between the error analysis and the correction formula, together with the numerical confirmation, represents a clear strength for applications in numerical analysis of Sturm-Liouville problems.

minor comments (2)
  1. [Section 2] The description of the Galerkin basis functions and the precise definition of the matrix entries in the spectral method could be expanded for readers unfamiliar with the weakly regular setting.
  2. [Section 4] A brief remark on the computational cost of assembling the correction formula relative to the original eigenvalue solve would strengthen the claim of 'low-cost'.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. The referee's summary correctly captures the main contributions: the Galerkin spectral matrix method, the convergence analysis, the derivation of the correction formula, and the numerical validation. Since the report lists no major comments, we have no point-by-point responses to provide.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper performs a convergence analysis of the Galerkin spectral matrix method for weakly regular Sturm-Liouville problems and derives a correction formula from those results. The abstract and description provide no evidence that the correction formula reduces to a fitted input, self-definition, or self-citation chain; the analysis is presented as supplying independent content that yields the formula. This is the most common honest outcome for a paper whose central claim rests on an explicit error analysis rather than tautological renaming or load-bearing self-reference.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No information available from the abstract to identify free parameters, axioms, or invented entities.

pith-pipeline@v0.9.0 · 5592 in / 984 out tokens · 25775 ms · 2026-05-24T22:09:15.612103+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages · 1 internal anchor

  1. [1]

    P. B. Bailey, W. N. Everitt, A. Zettl. The SLEIGN2 Sturm-Liouville C ode, ACM Trans. Math. Software , 21 (2001), 143-192. Available at: http://www.math.niu.edu/SL2/

    work page 2001

  2. [2]

    T. A. Driscoll, N. Hale, L. N. Trefethen, editors, Chebfun Guide, Pafnuty Publications, Oxford, 2014. Version used: 5.7.0

    work page 2014

  3. [3]

    Erd´ elyi, W

    A. Erd´ elyi, W. Magnus, F. Oberhettinger, F. G. Tricomi, Tables of Inte- gral Transforms, Vol. II . Based, in part, on notes left by Harry Bateman. McGraw-Hill Book Company, Inc., New York-Toronto-London, 195 4

    work page

  4. [4]

    Ghelardoni, C

    P. Ghelardoni, C. Magherini, A matrix method for fractional Stur m- Liouville problems on bounded domain, Adv. Comput. Math. 43 (2017) 1377-1401. 16 65 130 260 520 N -15 -13 -11 -9 -7 -5 relative error = 1/2, k = 25 65 130 260 520 N -15 -13 -11 -9 -7 -5 relative error = 3/4, k = 25 0 10 20 30 40 50 k -15 -13 -11 -9 -7 -5 relative error = 1/2, N = 149 0 10...

    work page 2017

  5. [5]

    Ledoux, M

    V. Ledoux, M. Van Daele, Matslise 2.0: a Matlab toolbox of Sturm-L iouville computations, ACM Trans. Math. Software 42 (2016), no. 4, Art. 29, 18 pp. Available at https://sourceforge.net/projects/matslise/

    work page 2016

  6. [6]

    Ledoux, Study of special algorithms for solving Sturm-Liouville an d Schr¨ odinger equations.Ph.D

    V. Ledoux, Study of special algorithms for solving Sturm-Liouville an d Schr¨ odinger equations.Ph.D. Thesis, Universiteit Gent, 2007

    work page 2007

  7. [7]

    A corrected spectral method for Sturm-Liouville problems with unbounded potential at one endpoint

    C. Magherini, A corrected spectral method for Sturm-Liouville p roblems with unbounded potential at one endpoint, J. Comput. Appl. Math (ac- cepted). Available at arXiv: 1812.02090

    work page internal anchor Pith review Pith/arXiv arXiv

  8. [8]

    Pruess, C

    S. Pruess, C. T. Fulton, Mathematical software for Sturm-Lio uville prob- lems, ACM Trans. on Math. Software , 19 (1993) 360-376. Release 2.2 (1993- 12-04). Available at http://www.netlib.org/misc/sledge

    work page 1993

  9. [9]

    Pryce, Numerical Solution of Sturm-Liouville Problems, Oxford: Univ

    J.D. Pryce, Numerical Solution of Sturm-Liouville Problems, Oxford: Univ. Press, London, 1993

    work page 1993

  10. [10]

    Shen, Efficient Spectral-Galerkin Method I

    J. Shen, Efficient Spectral-Galerkin Method I. Direct Solvers fo r the Second and Fourth Order Equations Using Legendre Polynomials, SIAM J. Sci. Comput. 15 (1994) no.6 1489-1505

    work page 1994

  11. [11]

    Sidi, Asymptotic expansion of Legendre series coefficients fo r functions with endpoint singularities, Asymptotic Analysis 65 (2009) 175-190

    A. Sidi, Asymptotic expansion of Legendre series coefficients fo r functions with endpoint singularities, Asymptotic Analysis 65 (2009) 175-190

    work page 2009

  12. [12]

    Townsend, M

    A. Townsend, M. Werb, S. Olver, Fast polynomial transforms b ased on Toeplitz and Hankel matrices, Math. Comp. 87, (2018) 1913-1934. 17 65 130 260 520 N -15 -13 -11 -9 -7 -5 relative error = 2/5, k = 25 65 130 260 520 N -15 -13 -11 -9 -7 -5 relative error = 4/5, k = 25 0 10 20 30 40 50 k -15 -13 -11 -9 -7 -5 relative error = 2/5, N = 149 0 10 20 30 40 50 ...

    work page 2018