pith. sign in

arxiv: 2604.20812 · v1 · submitted 2026-04-22 · 🧮 math.NA · cs.NA

Rigorous High-Order Hausdorff Dimension Estimation of Limit Sets of Continued Fraction Iterated Function Systems via B-Splines

Jacob Brown This is my paper

Pith reviewed 2026-05-09 23:15 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Hausdorff dimensioncontinued fractionsiterated function systemsB-splinesPerron-Frobenius operatorrigorous numerical boundsfinite element methodlimit sets
0
0 comments X

The pith

B-spline finite elements approximate the Perron-Frobenius operator to produce rigorous upper and lower bounds on Hausdorff dimensions of continued fraction limit sets.

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

The paper develops a numerical method to estimate the Hausdorff dimensions of limit sets generated by continued fraction iterated function systems. It uses B-splines in a finite element approximation of the associated Perron-Frobenius operator. The key step is proving that these approximations satisfy a hidden positivity property, which allows bounding the spectral radius that equals the dimension. This provides both rigorous bounds and higher-order convergence rates, as verified in one and two dimensions for quadratic splines.

Core claim

An analogue of the hidden positivity result for B-spline quasi-interpolants permits rigorous enclosure of the leading eigenvalue of the Perron-Frobenius operator, thereby giving computable upper and lower bounds on the Hausdorff dimension of the limit sets.

What carries the argument

B-spline quasi-interpolants approximating the transfer operator of the iterated function system, with the hidden positivity property ensuring bounds on the spectral radius.

If this is right

  • Upper and lower bounds on the dimension become computable with arbitrary precision through refinement.
  • Higher-order convergence rates improve efficiency over lower-order methods.
  • The approach applies directly to both one-dimensional and two-dimensional continued fraction systems.

Where Pith is reading between the lines

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

  • The same positivity argument may apply to other spline families or finite element bases in related dynamical systems.
  • Refinements of this method could yield dimension estimates for limit sets arising from more general iterated function systems in number theory.
  • Such rigorous numerics might help resolve open questions about the dimensions of specific continued fraction attractors.

Load-bearing premise

B-spline quasi-interpolants inherit an analogue of the hidden positivity property that is strong enough to control the spectral radius of the approximated operator.

What would settle it

A concrete continued fraction system whose Hausdorff dimension is known independently, for which the B-spline bounds fail to converge to that value or cross it.

Figures

Figures reproduced from arXiv: 2604.20812 by Jacob Brown.

Figure 2
Figure 2. Figure 2: shows an example of a knot sequence. The use of the index [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 2.1
Figure 2.1. Figure 2.1: shows an example of a knot sequence. The use of the index ℓ for the knot intervals comes from the fact that intervals of the form [ξk, ξk+1) can be degenerate in the case when the knot sequence has repeated knots. If the knot sequence has no repeated knots, then the indices ℓ and k coincide, in which case we will write [ξk, ξk+1) to denote the knot intervals. Definition 2.2. Given a knot sequence ξ = {ξk… view at source ↗
Figure 2.2
Figure 2.2. Figure 2.2: The parameter interval D2 ξ = [2, 5] when using quadratic B-splines on the domain [0, 7] with knot sequence ξ = {0, 1, 2, . . . , 7}. the B-splines “relevant” to x, i.e., k ∼ x = {k : bk(x) > 0}, and for a set Ω we write k ∼ Ω to denote the B-splines such that k ∼ x for some x ∈ Ω. The construction of a d-dimensional B-spline is rather straightforward, as it is simply the tensor product of one-dimensiona… view at source ↗
Figure 2.3
Figure 2.3. Figure 2.3: Example of the parameter rectangle for 2-dimensional B-splines. [PITH_FULL_IMAGE:figures/full_fig_p004_2_3.png] view at source ↗
read the original abstract

We develop a method for the rigorous estimation of Hausdorff dimensions of limit sets produced by continued fraction iterated function systems. Our method is based on the approximation of a Perron-Frobenius operator using the finite element method with B-splines as the choice of basis functions. This choice provides key numerical advantages including higher-order convergence and computational flexibility. We prove an analogue of Falk and Nussbaum's result on "hidden positivity" for B-spline quasi-interpolants to give rigorous upper and lower bounds for the Hausdorff dimensions of various limit sets. We provide numerical results to verify both the rigor and higher-order convergence of our method for quadratic B-spline interpolants in one and two dimensions.

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 / 1 minor

Summary. The manuscript develops a finite-element method using B-splines to approximate the Perron-Frobenius operator associated with continued-fraction iterated function systems. It proves an analogue of the Falk-Nussbaum hidden-positivity property for B-spline quasi-interpolants, thereby obtaining rigorous upper and lower bounds on the spectral radius and hence on the Hausdorff dimension of the limit sets. Numerical experiments are presented to verify both the rigor of the bounds and the higher-order convergence of the method for quadratic B-splines in one and two dimensions.

Significance. If the proof of the positivity analogue is correct and the error analysis fully accounts for the approximation, the work supplies a concrete advance: a high-order, rigorously justified computational scheme for Hausdorff dimensions that exploits the flexibility and convergence properties of B-splines. The explicit proof of the key positivity property together with the numerical verification of convergence rates constitute the principal strengths.

minor comments (1)
  1. Abstract: the text refers to 'B-spline quasi-interpolants' when stating the positivity result and to 'quadratic B-spline interpolants' when describing the numerical tests; a brief clarification of whether the same operators are used in both parts, or whether a distinction is intended, would remove potential ambiguity.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript, the recognition of its significance in providing a high-order rigorously justified scheme for Hausdorff dimensions via B-splines, and the recommendation for minor revision. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper's derivation proceeds by approximating the Perron-Frobenius operator via B-spline finite elements, then proving an original analogue of the Falk-Nussbaum hidden positivity property for the resulting quasi-interpolants; this positivity is used to bound the spectral radius that determines the Hausdorff dimension. The proof is presented as new work rather than obtained by fitting parameters to the target dimension or by reducing to prior self-citations. Numerical experiments serve only as verification of convergence order and rigor, not as the source of the bounds themselves. No equation or step is shown to be equivalent to its inputs by construction, and the central claim rests on approximation theory plus the newly established positivity analogue, which is independent of the final dimension values.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The approach rests on standard IFS theory and spline approximation results; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption Limit sets of continued fraction IFS are well-defined and their Hausdorff dimension is given by the logarithm of the spectral radius of the associated Perron-Frobenius operator.
    Invoked as the foundation for using the operator approximation to bound dimension.
  • domain assumption B-spline quasi-interpolants can be constructed to preserve positivity properties analogous to those in Falk-Nussbaum.
    Central to the rigorous bound claim.

pith-pipeline@v0.9.0 · 5414 in / 1217 out tokens · 36399 ms · 2026-05-09T23:15:43.090503+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

15 extracted references · 15 canonical work pages

  1. [1]

    On the dimension spectrum of infinite subsystems of continued fractions

    Vasileios Chousionis, Dmitriy Leykekhman, and Mariusz Urba´ nski. “On the dimension spectrum of infinite subsystems of continued fractions”. In:Transactions of the American Mathematical Society 373.2 (2020), pp. 1009–1042

    work page 2020

  2. [2]

    Rigorous Hausdorff dimension estimates for conformal fractals

    Vasileios Chousionis et al. “Rigorous Hausdorff dimension estimates for conformal fractals”. In:arXiv preprint arXiv:2408.06330(2024)

    work page arXiv 2024

  3. [3]

    Spline approximation by quasiinterpolants

    Carl De Boor and George J Fix. “Spline approximation by quasiinterpolants”. In:Journal of Approx- imation Theory8.1 (1973), pp. 19–45

    work page 1973

  4. [4]

    On polynomial approximation in Sobolev spaces

    Ricardo G Dur´ an. “On polynomial approximation in Sobolev spaces”. In:SIAM journal on numerical analysis20.5 (1983), pp. 985–988

    work page 1983

  5. [5]

    A new approach to numerical computation of Hausdorff di- mension of iterated function systems: applications to complex continued fractions

    Richard S Falk and Roger D Nussbaum. “A new approach to numerical computation of Hausdorff di- mension of iterated function systems: applications to complex continued fractions”. In:Integral Equa- tions and Operator Theory90.5 (2018), p. 61

    work page 2018

  6. [6]

    C m Eigenfunctions of Perron–Frobenius operators and a new approach to numerical computation of Hausdorff dimension: applications inR 1

    Richard S Falk and Roger D Nussbaum. “C m Eigenfunctions of Perron–Frobenius operators and a new approach to numerical computation of Hausdorff dimension: applications inR 1”. In:Journal of Fractal Geometry5.3 (2018), pp. 279–337

    work page 2018

  7. [7]

    Hidden positivity and a new approach to numerical com- putation of Hausdorff dimension: higher order methods

    Richard S Falk and Roger D Nussbaum. “Hidden positivity and a new approach to numerical com- putation of Hausdorff dimension: higher order methods”. In:Journal of Fractal Geometry9.1 (2021), pp. 23–72

    work page 2021

  8. [8]

    SIAM, 2013

    Klaus H¨ ollig and J¨ org H¨ orner.Approximation and modeling with B-splines. SIAM, 2013

    work page 2013

  9. [9]

    Fractals and self similarity

    John E Hutchinson. “Fractals and self similarity”. In:Indiana University Mathematics Journal30.5 (1981), pp. 713–747

    work page 1981

  10. [10]

    On the order of the approximation of functions by linear positive opera- tors

    Pavel Petrovich Korovkin. “On the order of the approximation of functions by linear positive opera- tors”. In:Doklady Akademii Nauk. Vol. 114. 6. Russian Academy of Sciences. 1957, pp. 1158–1161

    work page 1957

  11. [11]

    Lecture notes on Legendre polynomials: their origin and main properties

    F´ abio MS Lima. “Lecture notes on Legendre polynomials: their origin and main properties”. In:arXiv preprint arXiv:2210.10942(2022)

    work page arXiv 2022

  12. [12]

    Dimensions and measures in infinite iterated function systems

    R Daniel Mauldin and Mariusz Urba´ nski. “Dimensions and measures in infinite iterated function systems”. In:Proceedings of the London Mathematical Society3.1 (1996), pp. 105–154

    work page 1996

  13. [13]

    R Daniel Mauldin and Mariusz Urbanski.Graph directed Markov systems: geometry and dynamics of limit sets. Vol. 148. Cambridge University Press, 2003

    work page 2003

  14. [14]

    Hausdorff dimension estimates applied to Lagrange and Markov spectra, Zaremba theory, and limit sets of Fuchsian groups

    Mark Pollicott and Polina Vytnova. “Hausdorff dimension estimates applied to Lagrange and Markov spectra, Zaremba theory, and limit sets of Fuchsian groups”. In:Transactions of the American Math- ematical Society, Series B9.35 (2022), pp. 1102–1159

    work page 2022

  15. [15]

    Polynomial approximation on domains bounded by diffeomorphic images of graphs

    Ulrich Reif. “Polynomial approximation on domains bounded by diffeomorphic images of graphs”. In: Journal of Approximation theory164.7 (2012), pp. 954–970. 27

    work page 2012