pith. sign in

arxiv: 2605.23501 · v1 · pith:VL46CZ6Dnew · submitted 2026-05-22 · 🧮 math.NA · cs.NA

Spectral distribution of Jacobi weighted histopolation matrices via GLT theory

Allal Guessab , Federico Nudo , Stefano Serra-Capizzano This is my paper

Pith reviewed 2026-05-25 03:42 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords histopolationJacobi weightsGLT theoryspectral distributionmatrix factorizationnumerical stabilitytridiagonal couplingapproximation theory
0
0 comments X

The pith

Jacobi weighted histopolation matrices belong to the GLT class with explicit symbols that determine their spectral distributions.

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

The paper first constructs a reconstruction operator from weighted cell averages using primitives of Jacobi polynomials and derives an exact factorization of the resulting histopolation matrix through a backward-difference operator and a sampling operator. The sampling operator is further decomposed using a tridiagonal coupling matrix obtained from the three-term recurrence of the Jacobi polynomials. In the second part the authors show that the matrix sequences generated on regular meshes belong to the generalized locally Toeplitz class and compute the associated symbols, from which the asymptotic eigenvalue distributions follow directly.

Core claim

The histopolation matrices admit an exact factorization through a backward-difference operator and a sampling operator of Jacobi weighted primitives; the sampling operator further decomposes via a tridiagonal coupling matrix from the three-term recurrence. Consequently, under standard mesh-regularity assumptions, the induced matrix sequences belong to the GLT class with explicitly described symbols, from which the spectral distributions follow and implications for numerical stability are drawn.

What carries the argument

The GLT symbol of each histopolation matrix sequence, built from the weighted primitives and the mesh geometry, which encodes the asymptotic spectral behavior.

If this is right

  • The eigenvalues of the matrices are asymptotically distributed according to the measure defined by the GLT symbol.
  • The conditioning of the linear systems can be read off from the essential range of the symbol.
  • Iterative solver performance for the associated systems can be predicted as the mesh is refined.
  • The tridiagonal factorization provides an explicit route to computing the matrices and analyzing their spectra at finite size.

Where Pith is reading between the lines

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

  • The same GLT analysis could be applied to histopolation problems that use other families of orthogonal polynomials once the corresponding recurrence relations are available.
  • The explicit symbols may suggest choices of quadrature rules or basis modifications that improve conditioning without changing the approximation order.
  • If the mesh-regularity assumption is relaxed, the symbols might still be recovered by treating the irregularity as a low-rank perturbation.

Load-bearing premise

The discretization meshes satisfy the standard regularity conditions that guarantee the matrix sequences are in the GLT class.

What would settle it

For a sequence of successively refined regular meshes, compute the eigenvalues of the histopolation matrices and check whether their empirical distribution converges to the measure induced by the derived GLT symbol.

Figures

Figures reproduced from arXiv: 2605.23501 by Allal Guessab, Federico Nudo, Stefano Serra-Capizzano.

Figure 1
Figure 1. Figure 1: Behaviour of q (i) N (ε) for the scaled matrices H (1) N (left) and H (2) N (right), with N = 1000k, k = 1, . . . , 10, in the case α = β = 2. 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 10-4 10-3 10-2 10-1 100 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 10-4 10-3 10-2 10-1 100 [PITH_FULL_IMAGE:figures/full_fig_p039_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Behaviour of q (i) N (ε) for the scaled matrices H (3) N (left) and H (4) N (right), with N = 1000k, k = 1, . . . , 10, in the case α = β = 2. the scalings and for all the thresholds considered, the fraction of singular values greater than the threshold decreases as N increases. This suggests that the observed behaviour is not tied to the symmetric choice of the Jacobi parameters, but reflects the structur… view at source ↗
Figure 3
Figure 3. Figure 3: Behaviour of q (i) N (ε) for the scaled matrices H (1) N (left) and H (2) N (right), with N = 1000k, k = 1, . . . , 10, in the case α = 3/2 and β = 1. 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 10-4 10-3 10-2 10-1 100 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 10-4 10-3 10-2 10-1 100 [PITH_FULL_IMAGE:figures/full_fig_p040_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Behaviour of q (i) N (ε) for the scaled matrices H (3) N (left) and H (4) N (right), with N = 1000k, k = 1, . . . , 10, in the case α = 3/2 and β = 1. for j/N ∈ [0.8, 1], they rapidly approach zero. This suggests that the unscaled sequence {HN }N may admit a singular value distribution in accordance with Definition 3.1. In particular, from the graph, it can be conjectured that the symbol (if it exists) is … view at source ↗
Figure 5
Figure 5. Figure 5: Unscaled singular values of HN for the nodes xi = −1 + 2g(i/N), with g(y) = (e y − 1)/(e − 1) (left) and g(y) = y 2 (right). 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 10-2 10-1 100 101 102 103 104 105 [PITH_FULL_IMAGE:figures/full_fig_p041_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Singular values of NT (J) N for N = 2000 compared with the sorted samples of the modulus of its GLT symbol. scaled sequence {NT (J) N }N is the function defined in (45). In [PITH_FULL_IMAGE:figures/full_fig_p041_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Singular values of ∆N /N for the nodes xi = −1 + 2g(i/N), with g(y) = (e y − 1)/(e − 1) (left) and g(y) = y 2 (right). In [PITH_FULL_IMAGE:figures/full_fig_p042_7.png] view at source ↗
read the original abstract

In this paper we study a weighted histopolation problem on $[-1,1]$ associated with Jacobi weights. In the first part of the present work we prove results in approximation theory, while in the second we analyze the resulting matrices from an asymptotic linear algebra perspective. More in detail, in the first part, given weighted cell averages, we construct a reconstruction operator based on weighted primitives of Jacobi polynomials and investigate the resulting discretization matrices. At any fixed discretization level, we derive an exact factorization of the histopolation matrix through a backward-difference operator and a sampling operator of Jacobi weighted primitives. Combining a sharp integration by parts identity with the three-term recurrence of Jacobi polynomials, we further show that the primitive sampling operator admits an explicit decomposition involving a tridiagonal coupling matrix in the Jacobi spectral index. This yields a tridiagonal factor representation of the histopolation matrix. In the second part, under standard mesh-regularity assumptions, we show that all the various induced matrix sequences belong to the Generalized Locally Toeplitz (GLT) class, by describing in detail the related GLT symbols. As a consequence, we provide the corresponding spectral distributions and discuss their implications for numerical stability when solving the associated linear systems.

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

1 major / 0 minor

Summary. The paper examines weighted histopolation on [-1,1] using Jacobi weights. In the approximation-theory portion it constructs a reconstruction operator from weighted cell averages via weighted primitives of Jacobi polynomials, derives an exact factorization of the histopolation matrix through a backward-difference operator and a Jacobi-primitive sampling operator, and obtains a tridiagonal factor representation by combining integration-by-parts identities with the three-term recurrence. In the second part, under standard mesh-regularity assumptions, it asserts that the induced matrix sequences belong to the GLT class, supplies the corresponding GLT symbols, deduces the spectral distributions, and discusses consequences for numerical stability of the associated linear systems.

Significance. If the GLT membership and symbol derivations hold, the work supplies an explicit asymptotic spectral analysis for a family of histopolation matrices that arise in weighted approximation schemes. The algebraic factorizations (backward-difference plus tridiagonal coupling) are a concrete strength that could facilitate both theoretical and computational follow-up work on stability and preconditioning.

major comments (1)
  1. [second part] Second part (GLT analysis): the central claim that the matrix sequences generated by the Jacobi-weighted histopolation construction belong to the GLT class rests on 'standard mesh-regularity assumptions,' yet the manuscript provides no explicit verification that these assumptions are satisfied when the underlying mesh interacts with the Jacobi weight singularities at the endpoints. This verification is load-bearing for the passage from the algebraic factorization (valid at each fixed level) to the asymptotic GLT symbol and the resulting spectral distributions.

Simulated Author's Rebuttal

1 responses · 0 unresolved

Thank you for the opportunity to respond to the referee's report. We appreciate the positive assessment of the algebraic factorizations and the constructive feedback on the GLT analysis. We address the major comment as follows.

read point-by-point responses

  1. Referee: [second part] Second part (GLT analysis): the central claim that the matrix sequences generated by the Jacobi-weighted histopolation construction belong to the GLT class rests on 'standard mesh-regularity assumptions,' yet the manuscript provides no explicit verification that these assumptions are satisfied when the underlying mesh interacts with the Jacobi weight singularities at the endpoints. This verification is load-bearing for the passage from the algebraic factorization (valid at each fixed level) to the asymptotic GLT symbol and the resulting spectral distributions.

    Authors: We thank the referee for highlighting this point. The manuscript invokes 'standard mesh-regularity assumptions' as commonly used in the GLT literature for sequences of matrices arising from discretizations on non-uniform meshes (see, e.g., the references on GLT for variable coefficient problems). These typically include that the mesh is regular in the sense that the maximum mesh size tends to zero and the local mesh ratios are bounded independently of the level. Since the Jacobi weight is a fixed function (independent of the discretization parameter n) and belongs to L^1([-1,1]) for the admissible range of parameters, the weighted cell averages and the associated sampling operators remain well-defined and the localization procedure for deriving the GLT symbol carries through without additional restrictions. Nevertheless, we agree that an explicit statement confirming the applicability of these assumptions in the presence of possible endpoint singularities would improve clarity. In the revised manuscript, we will insert a short paragraph immediately after the statement of the mesh assumptions, providing this verification by recalling the precise conditions from the GLT theory and noting their compatibility with the Jacobi weight. revision: yes

Circularity Check

0 steps flagged

No circularity; derivations apply external GLT theory to algebraically derived factors under stated assumptions

full rationale

The paper first establishes exact algebraic factorizations of the histopolation matrix via backward differences, Jacobi primitive sampling, integration by parts, and the three-term recurrence (all at fixed level). It then invokes external GLT theory under 'standard mesh-regularity assumptions' to obtain symbols and spectral distributions. No step reduces a claimed prediction to a fitted input by construction, no self-definitional loop appears, and no load-bearing self-citation chain is present. The derivation chain remains independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on established properties of Jacobi polynomials and the GLT framework without introducing new free parameters, ad-hoc constants, or postulated entities.

axioms (2)
  • standard math Three-term recurrence and integration-by-parts identities for Jacobi polynomials
    Invoked to obtain the tridiagonal decomposition of the primitive sampling operator.
  • domain assumption Standard mesh-regularity assumptions sufficient for GLT membership
    Required to conclude that the matrix sequences belong to the GLT class and possess the stated symbols.

pith-pipeline@v0.9.0 · 5749 in / 1353 out tokens · 36752 ms · 2026-05-25T03:42:32.286241+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

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

  1. [1]

    A. C. Kak and M. Slaney.Principles of computerized tomographic imaging. SIAM, 2001

    work page 2001

  2. [2]

    Natterer.The mathematics of computerized tomography

    F. Natterer.The mathematics of computerized tomography. SIAM, 2001

    work page 2001

  3. [3]

    V. P. Palamodov.Reconstruction from integral data. CRC Press Boca Raton, FL, 2016

    work page 2016

  4. [4]

    Bosner, B

    T. Bosner, B. Crnkovi´ c, and J. ˇSkifi´ c. Application of CCC–Schoenberg operators on image resampling.BIT Numer. Math., 60:129–155, 2020

    work page 2020

  5. [5]

    Bruni Bruno, F

    L. Bruni Bruno, F. Dell’Accio, W. Erb, and F. Nudo. Polynomial Histopolation On Mock- Chebyshev Segments.J. Sci. Comput., 104:65, 2025. 43

    work page 2025

  6. [6]

    Bruni Bruno, F

    L. Bruni Bruno, F. Dell’Accio, W. Erb, and F. Nudo. Bivariate polynomial histopolation techniques on Padua, Fekete and Leja triangles.Adv. Comput. Math., 52(37), 2026

    work page 2026

  7. [7]

    Dell’Accio, A

    F. Dell’Accio, A. Guessab, M. Kbiri Alaoui, and F. Nudo. A general probability density framework for local histopolation and weighted function reconstruction from mesh line inte- grals.Numer. Algorithms, 2026

    work page 2026

  8. [8]

    F. Nudo. Function reconstruction using a Jacobi-weighted quadratic enriched histopolation method.Math. Comput. Simul., 245:512–529, 2026

    work page 2026

  9. [9]

    Dell’Accio, A

    F. Dell’Accio, A. Guessab, G. V. Milovanovi´ c, and F. Nudo. Nonconforming approximation methods for function reconstruction on general polygonal meshes via orthogonal polynomials. IMA J. Numer. Anal., 2026

    work page 2026

  10. [10]

    Demichelis

    V. Demichelis. Graphic applications of some interpolating weighted mean functions.Rocky Mt. J. Math., pages 1277–1286, 1995

    work page 1995

  11. [11]

    A. K. Bose. Functions satisfying a weighted average property.Trans. Amer. Math. Soc., 118:472–487, 1965

    work page 1965

  12. [12]

    T. Diagana. The existence of a weighted mean for almost periodic functions.Nonlinear Anal. Theory Methods Appl., 74:4269–4273, 2011

    work page 2011

  13. [13]

    I. Z. Pesenson. Average sampling and average splines on combinatorial graphs. In2019 13th International conference on Sampling Theory and Applications (SampTA), pages 1–4. IEEE, 2019

    work page 2019

  14. [14]

    I. Z. Pesenson. Weighted sampling and weighted interpolation on combinatorial graphs.arXiv preprint arXiv:1905.02603, 2019

    work page internal anchor Pith review Pith/arXiv arXiv 1905

  15. [15]

    Guessab and F

    A. Guessab and F. Nudo. Quadratic Weighted Histopolation on Tetrahedral Meshes with Probabilistic Degrees of Freedom.BIT Numer. Math. (accepted), 2026

    work page 2026

  16. [16]

    Garoni and S

    C. Garoni and S. Serra-Capizzano.Generalized Locally Toeplitz Sequences: Theory and Ap- plications: Volume I. Springer, Cham, 2018

    work page 2018

  17. [17]

    Garoni and S

    C. Garoni and S. Serra-Capizzano.Generalized Locally Toeplitz Sequences: Theory and Ap- plications: Volume II. Springer, Cham, 2018

    work page 2018

  18. [18]

    Dorostkar, M

    A. Dorostkar, M. Neytcheva, and S. Serra-Capizzano. Spectral analysis of coupled PDEs and of their Schur complements via generalized locally Toeplitz sequences in 2D.Comput. Methods Appl. Mech. Eng., 309:74–105, 2016

    work page 2016

  19. [19]

    Salinelli, S

    E. Salinelli, S. Serra-Capizzano, and D. Sesana. Eigenvalue-eigenvector structure of Schoenmakers–Coffey matrices via Toeplitz technology and applications.Linear Algebra Appl., 491:138–160, 2016

    work page 2016

  20. [20]

    Garoni, H

    C. Garoni, H. Speleers, S. Ekstr¨ om, A. Reali, S. Serra-Capizzano, and T. J.R. Hughes. Symbol- Based Analysis of Finite Element and Isogeometric B-Spline Discretizations of Eigenvalue Problems: Exposition and Review.Arch. Comput. Methods Eng., 26(5):1639–1690, 2019

    work page 2019

  21. [21]

    Benedusi, P

    P. Benedusi, P. Ferrari, M. E. Rognes, and S. Serra-Capizzano. Modeling excitable cells with the EMI equations: spectral analysis and iterative solution strategy.J. Sci. Comput., 98(3):58, 2024. 44

    work page 2024

  22. [22]

    Szeg˝ o.Orthogonal polynomials

    G. Szeg˝ o.Orthogonal polynomials. AMS Colloquium Publications, Vol. 23, 4th ed., 1975

    work page 1975

  23. [23]

    G. V. Milovanovi´ c. Orthogonal polynomial systems and some applications.Inner Product Spaces and Applications, Pitman Res. Notes Math. Ser, 376:115–182, 1997

    work page 1997

  24. [24]

    Erd´ elyi, W

    A. Erd´ elyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi.Higher Transcendental Func- tions. Krieger, Malabar, Fla. Vol. 1, 1981. Reprint of the 1953 edition

    work page 1981

  25. [25]

    Erd´ elyi, A

    T. Erd´ elyi, A. P. Magnus, and P. Nevai. Generalized Jacobi weights, Christoffel functions, and Jacobi polynomials.SIAM J. Math. Anal., 25:602–614, 1994

    work page 1994

  26. [26]

    Bruni Bruno and W

    L. Bruni Bruno and W. Erb. Polynomial interpolation of function averages on interval seg- ments.SIAM J. Numer. Anal., 62:1759–1781, 2024

    work page 2024

  27. [27]

    Bruni Bruno and S

    L. Bruni Bruno and S. Serra-Capizzano. On the conditioning of histopolation.arXiv preprint arXiv:2511.15395, 2025

    work page arXiv 2025

  28. [28]

    Bhatia.Matrix Analysis

    R. Bhatia.Matrix Analysis. Springer, New York, 1997

    work page 1997

  29. [29]

    Barbarino, D

    G. Barbarino, D. Bianchi, and C. Garoni. Constructive approach to the monotone rearrange- ment of functions.Expo. Math., 40:155–175, 2022

    work page 2022

  30. [30]

    Serra-Capizzano

    S-E Ekstr¨ om and S. Serra-Capizzano. Eigenvalues and eigenvectors of banded Toeplitz matrices and the related symbols.Numer. Linear Algebra Appl., 25:e2137, 2018

    work page 2018

  31. [31]

    Mazza, A

    M. Mazza, A. Ratnani, and S. Serra-Capizzano. Spectral analysis and spectral symbol for the 2D curl-curl (stabilized) operator with applications to the related iterative solutions.Math. Comput., 88:1155–1188, 2019

    work page 2019

  32. [32]

    Serra-Capizzano

    S. Serra-Capizzano. Distribution results on the algebra generated by Toeplitz sequences: a finite-dimensional approach.Linear Algebra Appl., 328:121–130, 2001

    work page 2001

  33. [33]

    Serra-Capizzano

    S. Serra-Capizzano. Generalized locally Toeplitz sequences: spectral analysis and applications to discretized partial differential equations.Special issue on structured matrices: analysis, algorithms and applications (Cortona, 2000), Linear Algebra Appl., 366:371–402, 2003

    work page 2000

  34. [34]

    Serra-Capizzano

    S. Serra-Capizzano. The GLT class as a generalized Fourier analysis and applications.Linear Algebra Appl., 419:180–233, 2006

    work page 2006

  35. [35]

    P. Tilli. Locally Toeplitz sequences: spectral properties and applications.Linear Algebra Appl., 278:91–120, 1998

    work page 1998

  36. [36]

    Barbarino, C

    G. Barbarino, C. Garoni, and S. Serra-Capizzano. Block generalized locally Toeplitz sequences: theory and applications in the unidimensional case.Electron. Trans. Numer. Anal., 53:28–112, 2020

    work page 2020

  37. [37]

    Barbarino, C

    G. Barbarino, C. Garoni, and S. Serra-Capizzano. Block generalized locally Toeplitz sequences: theory and applications in the multidimensional case.Electron. Trans. Numer. Anal., 53:113– 216, 2020

    work page 2020

  38. [38]

    Serra-Capizzano and C

    S. Serra-Capizzano and C. Tablino-Possio. Analysis of preconditioning strategies for collocation linear systems.Linear Algebra Appl., 369:41–75, 2003. 45

    work page 2003

  39. [39]

    Benedusi, S

    P. Benedusi, S. Riva, L. Belluzzi, and S. Serra-Capizzano. Analysis of eigenvalue clustering leads to optimal scaling in numerical radiative transfer.arXiv preprint arXiv:2602.21958, 2026

    work page arXiv 2026

  40. [40]

    B. Carl. On a Weyl inequality of operators in Banach spaces.Proc. Am. Math. Soc., 137:155– 159, 2009. 46

    work page 2009