Fractional backward spectral approximation theory for weakly singular adjoint integral equations

Mahmoud A. Zaky

arxiv: 2605.15825 · v2 · pith:YVKUW6EVnew · submitted 2026-05-15 · 🧮 math.NA · cs.NA

Fractional backward spectral approximation theory for weakly singular adjoint integral equations

Mahmoud A. Zaky This is my paper

Pith reviewed 2026-05-19 22:20 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords spectral approximationweakly singular integral equationsJacobi polynomialsendpoint singularityfractional differential equationsweighted Sobolev spacesVolterra equationscollocation methods

0 comments

The pith

Composing Jacobi polynomials with an endpoint algebraic mapping yields bases that preserve terminal singularities for high-order spectral approximations.

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

The paper develops a fractional approximation framework for functions whose regularity is limited near a terminal point. The basis arises from composing classical Jacobi polynomials with an algebraic mapping that embeds the terminal singularity directly into the approximation space. Structural properties including orthogonality, derivative identities, and a singular Sturm-Liouville formulation are derived, along with weighted Sobolev spaces and corresponding projection and interpolation error bounds. These results supply the analytic foundation needed for spectral and collocation methods applied to endpoint-singular and weakly regular problems.

Core claim

The central claim is that the fractional polynomials obtained by composing classical Jacobi polynomials with an endpoint algebraic mapping incorporate the terminal singular structure into the approximation space. This composition produces orthogonality relations, derivative identities, and a singular Sturm-Liouville eigenvalue problem. The associated weighted Sobolev spaces then admit explicit projection and Gauss-type interpolation error estimates in weighted norms, together with inverse inequalities and embedding results, thereby justifying high-order convergence for spectral approximations of terminal-value problems, fractional differential equations, and weakly singular Volterra integral

What carries the argument

Fractional polynomials formed by composing classical Jacobi polynomials with an endpoint algebraic mapping that directly incorporates the terminal singular structure into the approximation space.

Load-bearing premise

That composing classical Jacobi polynomials with an endpoint algebraic mapping directly yields the claimed orthogonality, derivative identities, and weighted error estimates for functions with limited regularity near the terminal point.

What would settle it

Numerical computation of the weighted-norm error for a known exact solution of a terminal-value fractional differential equation, checking whether the observed rate of decrease matches the predicted high-order rate as the polynomial degree increases.

Figures

Figures reproduced from arXiv: 2605.15825 by Mahmoud A. Zaky.

**Figure 1.** Figure 1: Errors in the L ∞(J)- and L 2 κ−1/4,−1/4,ρ (J)-norms versus the ρ-polynomial degree N: (a) θ = 0.5, ρ = 1 2 ; (b) θ = 0.5, ρ = 1 4 ; (c) θ = 2 3 , ρ = 1 3 ; (d) θ = 2 3 , ρ = 1 6 . These test problems are used to evaluate the capability of the proposed Zaky–Jacobi spectral method to approximate solutions with reduced regularity near the terminal endpoint t = 1. For arbitrary choices of γ1 and γ2, the trans… view at source ↗

**Figure 2.** Figure 2: Comparison of the weighted L 2 κ0,0,ρ (J)-error and the L ∞(J)-error for the test problem in case (i), with γ1 = √ 2, γ2 = √ 3, and different values of θ and ρ. 7. Conclusion In this paper, we developed a Zaky–Jacobi approximation framework for functions with limited regularity near the terminal point. The proposed basis is constructed by composing classical Jacobi polynomials with an endpoint algebraic ma… view at source ↗

**Figure 3.** Figure 3: Comparison of the weighted L 2 κ0,0,ρ (J)-error and the L ∞(J)-error for the test problem in case (ii), with γ1 = √ 2, γ2 = √ 3, and different values of θ and ρ. smooth polynomial bases are applied to weakly regular functions. A rigorous approximation theory was established for the Zaky–Jacobi polynomials. We derived their main structural properties, including orthogonality relations, recurrence formulas,… view at source ↗

read the original abstract

We introduce a new class of fractional backward orthogonal functions designed for the spectral approximation of weakly singular adjoint Volterra integral equations. These basis functions generate an approximation space that naturally reflects the terminal-endpoint singular behaviour produced by weakly singular kernels. We develop the basic approximation theory for the proposed backward orthogonal basis, including weighted projection estimates, Gauss-type interpolation estimates, inverse inequalities, and stability bounds for the associated weakly singular adjoint integral operator. The error analysis and numerical results show that the proposed backward Jacobi method is particularly suitable for solutions with terminal-endpoint weak singularities and can recover high-order convergence rates that are typically lost when usual polynomial approximations are applied directly to such weakly regular solutions.

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.

Desk Editor's Note private letter to a colleague

The paper gives a mapped Jacobi basis that bakes endpoint singularities into the approximation space for weakly singular integrals, with supporting weighted error theory that appears new but needs the derivations checked.

read the letter

Hi, the central point is that this constructs fractional polynomials by algebraically mapping Jacobi polynomials to the endpoint so the basis itself reflects the terminal singularity. That setup is then used to prove orthogonality, derivative identities, a singular Sturm-Liouville form, and weighted Sobolev error estimates for projection and interpolation. The result is positioned as a foundation for spectral methods on terminal-value problems, fractional DEs, and weakly singular Volterra equations. If the mapping works as claimed, it supplies a direct way to get high-order convergence without extra grading or special weights chosen by hand. The construction and the move to weighted spaces look like the genuinely new piece rather than a routine extension. The abstract states the proofs are there, including inverse inequalities and embeddings, so the framework is at least internally consistent on paper. The soft spot is whether the specific algebraic map really delivers the full set of weighted estimates for functions whose regularity drops only near the endpoint; the constants and any restrictions on the singularity strength would need to be inspected in the proofs. No circularity or fitted parameters show up. This is written for numerical analysts already working on high-order methods for singular integral equations and fractional models. A reader in that subfield would get concrete tools to try or extend. It is worth a serious referee because the claims are specific enough to verify and the target applications are common. I would send it out for review and ask for the explicit mapping formula plus at least one set of numerical rates to confirm the theory translates to practice.

Referee Report

2 major / 3 minor

Summary. The manuscript introduces a fractional approximation framework for functions with limited regularity near the terminal point. The basis is constructed by composing classical Jacobi polynomials with an endpoint algebraic mapping to incorporate the terminal singular structure. Structural properties including orthogonality relations, derivative identities, and a singular Sturm-Liouville eigenvalue formulation are established. Weighted Sobolev spaces are defined, and projection and Gauss-type interpolation error estimates in weighted norms are proved, together with inverse inequalities and weighted Sobolev embedding estimates. The theory is presented as a foundation for high-order spectral and collocation approximations of endpoint-singular and weakly regular problems such as terminal value problems, fractional differential equations, and weakly singular Volterra integral equations.

Significance. If the central constructions and estimates hold, the work supplies a direct, mapping-based route to high-order spectral methods for a practically important class of singular problems in numerical analysis. The explicit derivation of orthogonality, derivative identities, and weighted error bounds from a classical polynomial family constitutes a clear technical contribution that could support improved convergence analysis for the listed applications.

major comments (2)

[§3] §3 (Basis construction and structural properties): The passage from the classical Jacobi orthogonality to the claimed orthogonality of the mapped fractional polynomials with respect to the singular weight must be shown explicitly. The algebraic mapping (presumably of the form (1+x)^β or equivalent) needs to be stated with its precise exponent, and the proof that this mapping preserves the required weighted inner-product orthogonality and yields the singular Sturm-Liouville formulation should be written out in full; without this step the subsequent projection and interpolation estimates rest on an unverified transfer of properties.
[Theorem 4.2] Theorem 4.2 (projection error estimates): The weighted Sobolev error bounds are derived under the assumption that the mapped basis reproduces the correct singularity. The proof should include an explicit dependence of the constants on the singularity strength parameter; if the constants blow up or become non-uniform as the regularity index approaches the endpoint singularity, the claimed foundation for high-order approximation of weakly singular Volterra equations is weakened.

minor comments (3)

[§4] Notation for the weighted norms in §4 should be introduced with a short comparison table to the classical Sobolev norms to prevent reader confusion.
[§4.3] The statement of the Gauss-type interpolation operator in §4.3 would benefit from an explicit formula for the quadrature nodes induced by the mapped basis.
[Introduction] A brief remark in the introduction on how the present construction differs from existing mapped polynomial or fractional polynomial approaches (e.g., those based on Jacobi or generalized Jacobi weights) would improve context.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and the positive evaluation of our manuscript. We address the major comments point by point below, and we will make the necessary revisions to strengthen the presentation.

read point-by-point responses

Referee: [§3] §3 (Basis construction and structural properties): The passage from the classical Jacobi orthogonality to the claimed orthogonality of the mapped fractional polynomials with respect to the singular weight must be shown explicitly. The algebraic mapping (presumably of the form (1+x)^β or equivalent) needs to be stated with its precise exponent, and the proof that this mapping preserves the required weighted inner-product orthogonality and yields the singular Sturm-Liouville formulation should be written out in full; without this step the subsequent projection and interpolation estimates rest on an unverified transfer of properties.

Authors: We agree that an explicit derivation is essential for clarity. In the revised version, we will state the algebraic mapping precisely (of the form used in the manuscript to incorporate the terminal singularity) and provide a step-by-step proof of the orthogonality by change of variables, showing how the weighted inner product transforms to the standard Jacobi inner product. The singular Sturm-Liouville problem will be derived by applying the chain rule and transformation to the classical Jacobi differential equation. These details will be inserted in Section 3. revision: yes
Referee: [Theorem 4.2] Theorem 4.2 (projection error estimates): The weighted Sobolev error bounds are derived under the assumption that the mapped basis reproduces the correct singularity. The proof should include an explicit dependence of the constants on the singularity strength parameter; if the constants blow up or become non-uniform as the regularity index approaches the endpoint singularity, the claimed foundation for high-order approximation of weakly singular Volterra equations is weakened.

Authors: We will revise Theorem 4.2 and its proof to explicitly display the dependence of the error constants on the singularity parameter β. Our analysis shows that the constants remain bounded independently of the polynomial degree N for fixed β > 0, and we will include the explicit form C(β) in the estimates. This ensures the high-order convergence for the applications to weakly singular Volterra equations, as the singularity strength is typically fixed in those problems. We do not observe blow-up in the relevant regime. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is a direct construction with independent proofs

full rationale

The paper defines a fractional polynomial basis explicitly by composing classical Jacobi polynomials with an endpoint algebraic mapping, then states that the main structural properties (orthogonality relations, derivative identities, singular Sturm-Liouville formulation) are established and that projection/interpolation error estimates in weighted Sobolev spaces are proved. These steps are presented as standard mathematical derivations from the given construction rather than reductions to fitted parameters, self-definitions, or self-citation chains. No equations or claims in the provided abstract reduce any central result to its inputs by construction; the framework remains self-contained against external benchmarks of polynomial approximation theory.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim rests on the algebraic mapping construction, the transfer of Jacobi properties, and the subsequent analytic proofs in weighted spaces; no numerical fitting or external data are invoked.

axioms (1)

standard math Standard algebraic and analytic properties of classical Jacobi polynomials on [-1,1]
Invoked as the base objects that are composed with the endpoint mapping.

invented entities (1)

Fractional polynomials obtained by endpoint algebraic mapping no independent evidence
purpose: To embed the terminal singularity directly into the approximation space
New objects constructed in the paper; no independent experimental evidence supplied.

pith-pipeline@v0.9.0 · 5648 in / 1284 out tokens · 31215 ms · 2026-05-19T22:20:42.210057+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The proposed basis is constructed by composing classical Jacobi polynomials with an endpoint algebraic mapping... orthogonality relations, derivative identities, and a singular Sturm–Liouville eigenvalue formulation... weighted Sobolev spaces... projection and Gauss-type interpolation error estimates
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Zaky–Jacobi polynomial of degree r is defined by P_{μ,υ,ρ}^r(t) = P_{μ,υ}^r(1−2(1−t)^ρ)

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.