pith. sign in

arxiv: 2604.11093 · v2 · submitted 2026-04-13 · 🧮 math.NA · cs.NA

A discontinuous Galerkin method with fractal elements

Sergio G\'omez , David Hewett , Andrea Moiola This is my paper

Pith reviewed 2026-05-10 15:56 UTC · model grok-4.3

classification 🧮 math.NA cs.NA MSC 65N30
keywords discontinuous Galerkin methodfractal boundaryKoch snowflakefinite element methodself-similarityPoisson equationDirichlet eigenvalue problemquasi-optimality
0
0 comments X

The pith

A discontinuous Galerkin method defined directly on fractal elements approximates elliptic problems in the Koch snowflake domain without polygonal prefractal 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 discontinuous Galerkin finite element method that operates on the exact Koch snowflake geometry by tiling the domain with smaller self-similar fractal elements. Fluxes across the fractal inter-element boundaries are handled weakly through integrals over element subdomains rather than direct surface integrals. These integrals are evaluated exactly for polynomial basis functions by exploiting the self-similarity of the elements. The authors establish well-posedness and quasi-optimality of the resulting discrete problem and supply a partial convergence analysis, with numerical tests confirming the approach for linear and quadratic elements on both the Poisson equation and the related eigenvalue problem.

Core claim

We formulate, analyse, and implement a discontinuous Galerkin finite element method on the Koch snowflake domain itself, using a geometry-conforming mesh of fractal elements each similar to the snowflake. Fluxes across inter-element fractal curves are represented weakly by integrals over element subdomains, which for local polynomial bases can be evaluated exactly via self-similarity. The method is shown to be well-posed and quasi-optimal, with a partial convergence theory, and numerical results demonstrate its effectiveness for piecewise linear and quadratic approximations of the Poisson problem and the Dirichlet eigenvalue problem.

What carries the argument

A geometry-conforming fractal tiling of the Koch snowflake domain in which each element is a scaled copy of the snowflake, combined with a weak representation of inter-element fluxes as subdomain integrals that are computed exactly using self-similarity for polynomial basis functions.

If this is right

  • The method produces a well-posed and quasi-optimal discrete problem without first replacing the fractal boundary by a sequence of polygons.
  • Flux integrals over fractal interfaces are replaced by exact subdomain integrals whose evaluation cost is independent of the fractal iteration depth for polynomial bases.
  • The same framework applies directly to the Dirichlet eigenvalue problem on the snowflake.
  • Partial convergence rates follow from the quasi-optimality result once an approximation property for the fractal-element spaces is established.

Where Pith is reading between the lines

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

  • The same subdomain-integral technique could be tested on other self-similar fractals whose tiles admit an exact scaling relation for polynomial integrands.
  • If the exact integral property extends to higher-degree polynomials, the method would immediately support higher-order DG approximations on the same meshes.
  • The approach suggests a route to adaptive refinement that respects the fractal geometry at every level rather than stopping at a fixed prefractal stage.

Load-bearing premise

The weak integrals representing fluxes across fractal inter-element boundaries can be evaluated exactly using self-similarity when the local basis functions are polynomials.

What would settle it

A sequence of numerical solutions on successively refined fractal tilings that fails to converge in the energy norm to the known analytic solution of the Poisson problem on the Koch snowflake.

Figures

Figures reproduced from arXiv: 2604.11093 by Andrea Moiola, David Hewett, Sergio G\'omez.

Figure 1
Figure 1. Figure 1: The Koch snowflake domain and illustrations of two different meshing strategies. In this [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Sequence of prefractals approximating the Koch curve [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Sequence of prefractals approximating the Koch snowflake [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The boundary of Ω is the union of either six copies of Γ, or three copies of Γ scaled by a factor of √ 3. 5 [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The self-similar decomposition of the Koch snowflake into seven smaller snowflakes [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: (a) Each pair of adjacent elements K± is the image of s1(Ω) and s2(Ω) under a similarity S. (b) A zoom of the mesh in [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The meshes Tℓ for ℓ = 1, 2, 3 (a–c), T ′ ℓ for ℓ = 2, 3, 4 (d–f), and T ′ 3,ℓ∗ for ℓ ∗ = 1, 2, 3 (g–i); see §3.5. 10 [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: (a) Elements K± and wedges ▼, ▲ associated with a face F = ∂K− ∩ ∂K+ ∈ FI . The quadrilateral ♦ is the interior of the set ▼ ∪ ▲, and its boundary ♢ is the thick black line. (b) Zoom of an LQU mesh showing division of elements into 3, 4, 5, and 6 wedges. The union of the (3, 4, 5, or 6) red segments contained in K ∈ T is denoted ∗K in §6.3. The angle between the segments S − 1 and S − 2 is π/3, so ▼ takes … view at source ↗
Figure 9
Figure 9. Figure 9: For each boundary face F, the wedge ▲ is contained in the domain complement R 2 \ Ω. To summarise, so far we have shown that, if w ∈ V µ ∩ H1 0 (Ω) and vh ∈ Vh, then aSIP(w, vh) = X K∈T Z K ∇w · ∇vh dx − 1 2 X F ∈FI h I▼(w, v− h − v + h ) − I▲(w, v− h − v + h ) i − X F ∈FB I▼(w, vh). For such w and vh, noting that (4.3) holds also for w ∈ V µ and v ∈ P p (R 2 ) (and not just for w, v ∈ P p (R 2 )), for eac… view at source ↗
Figure 10
Figure 10. Figure 10: The subsets Ξ1, . . . , Ξ6 (in red) of ∂Ω introduced in Assumption 6.7. The above observations motivate the following approximation assumption. 22 [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: h-convergence for the problem with smooth (Gaussian) solution (see §8.1) in the DG norm (left panel) and the L2(Ω) norm (right panel). The numbers in yellow rectangles are the associated empirical convergence rates. where the factor h −d F canceled out after the change of variables. Since (ψ −1 Km ◦ ξF ) : Γ → Fb− and (ψ −1 Kn ◦ ξF ) : Γ → Fb+, only evaluations of the basis functions {ϕbi} on the faces Fb… view at source ↗
Figure 12
Figure 12. Figure 12: (a) Discrete approximation uh with p = 2 on the mesh T ′ 3,5 for the problem with f = 1 (see §8.2). (b) Errors in (8.1) computed on the mesh sequences {T ′ ℓ } 9 ℓ=0 and {T ′ 3,ℓ∗ } 5 ℓ ∗=0. The number of degrees of freedom (N) refers to the approximation uhi . to (2.1) is given by the Gaussian function u(x, y) = exp(−(x 2 + y 2 )/σ2 ), with σ = 10−1 . This solution does not vanish exactly on ∂Ω, but is n… view at source ↗
Figure 13
Figure 13. Figure 13: Condition number of the Galerkin matrix ASIP for the sequence of quasi-uniform meshes {T ′ ℓ } 8 ℓ=0 and the sequence of boundary-refined meshes {T ′ 3,ℓ∗ } 5 ℓ ∗=0 (see §8.3). sequences {T ′ ℓ } 8 ℓ=0 and {T ′ 3,ℓ∗ } 5 ℓ ∗=0, with approximations of degree p = 1 and p = 2. For the quasi￾uniform sequence, the condition number grows like O(N) (which corresponds to O(h −2 )), as for the SIP-DG-FEM on standar… view at source ↗
Figure 14
Figure 14. Figure 14: Discrete eigenfunctions computed with approximations of degree [PITH_FULL_IMAGE:figures/full_fig_p032_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: (a) The shaded region is the quadrilateral [PITH_FULL_IMAGE:figures/full_fig_p035_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Every H1 (▼) function can be extended to an H1 (K) function by even reflections. Proof. Let w ∈ H 2,2 µ (Ω). By the standard Hölder inequality (R Ω |fg| dx ≤ ∥f∥Lr(Ω)∥g∥Lr′ (Ω) with 1 r + 1 r ′ = 1, here with f = |D2w| p δ µp , g = δ −µp , r = 2 p ) we have Z Ω |D2w| p dx 1/p ≤ Z Ω |D2w| 2 δ 2µ dx 1/2 Z Ω δ −qµ dx 1/q , (A.6) when 1 ≤ p < 2 and 2 ≤ q < ∞ with 1 q + 1 2 = 1 p , i.e. q = 2p 2−p . Thus… view at source ↗
Figure 17
Figure 17. Figure 17: (a) Decomposition of the snowflake Ω into six wedges W1, . . . , W6. (b) Decomposition of the wedge W1 into the equilateral triangle T and four wedges W1,1, . . . , W1,4. We start by considering integrals of polynomials over the reference wedges W1, . . . , W6 illustrated in [PITH_FULL_IMAGE:figures/full_fig_p038_17.png] view at source ↗
read the original abstract

We formulate, analyse, and implement a discontinuous Galerkin finite element method (DG-FEM) for the approximation of the solution of an elliptic boundary value problem in a domain with fractal boundary. We consider the case of the Poisson equation in the Koch snowflake domain with zero Dirichlet boundary conditions, but our methodology can be generalised to other cases. Rather than first approximating the snowflake domain by a polygonal "prefractal" and then applying a standard DG-FEM on the prefractal, we define a DG-FEM on the snowflake itself, using a geometry-conforming mesh (a fractal tiling) consisting of fractal elements, each similar to the original snowflake. Fluxes across inter-element boundaries, which are fractal curves, are represented in a weak way by integrals over element subdomains. We show how, for local polynomial basis functions, these integrals can be evaluated exactly using the similarity of the elements. We prove well-posedness and quasi-optimality of the method, and provide a partial convergence analysis. We present numerical results for piecewise linear and piecewise quadratic basis functions, which demonstrate the effectiveness of the method. We also apply our method to the related Dirichlet eigenvalue problem.

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

3 major / 2 minor

Summary. The manuscript formulates, analyzes, and implements a discontinuous Galerkin finite element method for the Poisson equation (and related eigenvalue problem) on the Koch snowflake domain. It employs a geometry-conforming mesh of fractal elements similar to the domain itself, represents fluxes across fractal inter-element boundaries weakly via integrals over element subdomains, and evaluates those integrals exactly for local polynomial bases by exploiting self-similarity. Well-posedness and quasi-optimality are proved, a partial convergence analysis is given, and numerical results are shown for piecewise-linear and piecewise-quadratic elements.

Significance. If the central consistency and exact-evaluation claims hold, the work offers a direct, non-approximating approach to DG methods on self-similar fractal domains that avoids prefractal polygonal approximations; the explicit well-posedness and quasi-optimality proofs together with the exact self-similar quadrature for polynomials constitute a clear technical contribution, while the numerical experiments provide concrete evidence of practical effectiveness.

major comments (3)
  1. [Formulation section (abstract and §2)] The replacement of direct integration over fractal inter-element curves by weak subdomain integrals is load-bearing for the entire analysis (well-posedness, quasi-optimality, and partial convergence). The manuscript must supply a precise statement, in the appropriate trace space on the fractal boundary, showing that this substitution recovers the correct distributional flux; without it the discrete bilinear form analyzed in the proofs is not demonstrably the one implemented.
  2. [Implementation of the bilinear form (abstract and §3)] The claim that the subdomain integrals admit exact closed-form evaluation for polynomial bases via the Koch snowflake self-similarity recursion must be verified to terminate exactly rather than only approximately; if the infinite iterative construction prevents exact termination, the numerical results cannot be taken as confirmation of the theory.
  3. [Convergence analysis (abstract and §4)] The convergence analysis is described as partial and no explicit error bounds or rates are visible. This is a load-bearing gap: quasi-optimality alone does not guarantee that the method converges as the fractal mesh is refined, and the numerical results therefore lack theoretical anchoring.
minor comments (2)
  1. [Abstract] The abstract states that the methodology 'can be generalised to other cases' but gives no indication of which other problems or boundary conditions are intended; a single sentence would clarify scope.
  2. [Mesh construction] Notation for the fractal tiling and the self-similar subdomains should be introduced once with a clear diagram or recursive definition to avoid ambiguity when the integrals are later evaluated.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive major comments. We address each point below and will revise the manuscript to incorporate the requested clarifications and extensions while preserving the core contributions on the geometry-conforming DG method with exact self-similar evaluation.

read point-by-point responses

  1. Referee: [Formulation section (abstract and §2)] The replacement of direct integration over fractal inter-element curves by weak subdomain integrals is load-bearing for the entire analysis (well-posedness, quasi-optimality, and partial convergence). The manuscript must supply a precise statement, in the appropriate trace space on the fractal boundary, showing that this substitution recovers the correct distributional flux; without it the discrete bilinear form analyzed in the proofs is not demonstrably the one implemented.

    Authors: We agree that an explicit equivalence statement is necessary for rigor. In the revised manuscript we will insert a new lemma in Section 2 that, working in the trace space dual to the appropriate fractional Sobolev space on the Koch curve (accounting for its Hausdorff dimension), shows via integration by parts on the self-similar subdomains that the weak subdomain-integral representation of the numerical flux coincides with the distributional interface term. This will confirm that the analyzed bilinear form is identical to the one implemented. revision: yes

  2. Referee: [Implementation of the bilinear form (abstract and §3)] The claim that the subdomain integrals admit exact closed-form evaluation for polynomial bases via the Koch snowflake self-similarity recursion must be verified to terminate exactly rather than only approximately; if the infinite iterative construction prevents exact termination, the numerical results cannot be taken as confirmation of the theory.

    Authors: The recursion terminates exactly for any fixed polynomial degree. Because the basis functions are polynomials of degree at most p, each application of the self-similarity map scales the integrand by a known factor and maps it to a polynomial of the same degree on a smaller copy; after a finite number of steps (at most p+1 per similarity branch) the remaining integrals reduce to explicit moments on the initial equilateral-triangle generators, which are computed in closed form. We will add a short inductive proof and pseudocode in Section 3 demonstrating finite termination and confirming that the reported numerics employ this exact procedure rather than truncation. revision: yes

  3. Referee: [Convergence analysis (abstract and §4)] The convergence analysis is described as partial and no explicit error bounds or rates are visible. This is a load-bearing gap: quasi-optimality alone does not guarantee that the method converges as the fractal mesh is refined, and the numerical results therefore lack theoretical anchoring.

    Authors: We accept that the present analysis is only partial and that quasi-optimality must be supplemented by approximation properties to obtain convergence under refinement. In the revised Section 4 we will derive explicit a priori error estimates by combining the existing quasi-optimality result with interpolation estimates on the fractal elements; these estimates follow from the self-similar subdivision and standard polynomial approximation theory applied level-by-level to the prefractal polygons. The resulting bounds will be stated in terms of the mesh-size parameter (subdivision level) and will directly anchor the observed numerical convergence for both linear and quadratic elements. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper defines a DG-FEM directly on the Koch snowflake domain using a fractal tiling mesh, represents inter-element fluxes weakly via subdomain integrals as a modeling choice to avoid direct fractal-curve integration, and derives exact evaluation formulas for polynomial bases from the independent geometric self-similarity of the Koch snowflake. Well-posedness, quasi-optimality, and partial convergence are then proved by adapting standard DG theory (bilinear form coercivity, continuity, and approximation properties) to this setting. No step reduces a claimed prediction or theorem to a fitted parameter, self-citation chain, or definitional tautology; the exact-integral claim is a constructive computation shown from the similarity recursion, not presupposed in the analysis. The derivation remains self-contained against external DG benchmarks and the known fractal geometry.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the self-similarity property of the Koch snowflake to enable exact integral evaluation; standard elliptic PDE well-posedness theory is assumed without new axioms introduced in the abstract.

axioms (1)
  • domain assumption The Koch snowflake domain admits a self-similar fractal tiling by smaller copies of itself
    Invoked to define the geometry-conforming mesh and to evaluate subdomain integrals exactly for polynomial bases.

pith-pipeline@v0.9.0 · 5507 in / 1186 out tokens · 37966 ms · 2026-05-10T15:56:05.829615+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

69 extracted references · 69 canonical work pages

  1. [1]

    Achdou and T

    Y. Achdou and T. Deheuvels , A transmission problem across a fractal self-similar interface , SIAM J. Multiscale Model. Simul., 14 (2016), pp. 708–736

    work page 2016

  2. [2]

    Achdou, C

    Y. Achdou, C. Sabot, and N. Tchou , Transparent boundary conditions for the Helmholtz equation in some ramified domains with a fractal boundary , J. Comput. Phys., 220 (2007), pp. 712–739

    work page 2007

  3. [3]

    D. N. Arnold , An interior penalty finite element method with discontinuous elements , SIAM J. Numer. Anal., 19 (1982), pp. 742–760

    work page 1982

  4. [4]

    Babuška and B

    I. Babuška and B. Q. Guo , The h-p version of the finite element method for domains with curved boundaries, SIAM J. Numer. Anal., 25 (1988), pp. 837–861

    work page 1988

  5. [5]

    Bagnerini, A

    P. Bagnerini, A. Buffa, and E. V acca , Finite elements for a prefractal transmission problem, C. R. Math., 342 (2006), pp. 211–214

    work page 2006

  6. [6]

    Bagnerini, A

    P. Bagnerini, A. Buffa, and E. V acca , Mesh generation and numerical analysis of a Galerkin method for highly conductive prefractal layers , Appl. Numer. Math., 65 (2013), pp. 63– 78

    work page 2013

  7. [7]

    Bandt , Self-similar sets 5

    C. Bandt , Self-similar sets 5. Integer matrices and fractal tilings of Rn, Proc. Amer. Math. Soc., (1991), pp. 549–562

    work page 1991

  8. [8]

    Banjai , Eigenfrequencies of fractal drums , J

    L. Banjai , Eigenfrequencies of fractal drums , J. Comput. Appl. Math., 198 (2007), pp. 1–18. 41 Γ, Ω Koch curve and snowflake §2.1 dimH(·), d = dim H(Γ) Hausdorff dimension §2.1 H 1 0 (Ω), H1(Ω), W 1 p (Ω) Sobolev spaces §2.2–2.3 C0,α(Ω), H2,2 µ (Ω), D2 Hölder and weighted Sobolev space, second derivatives §2.3 δ, µ Boundary distance function and integra...

    work page 2007

  9. [9]

    Bannister, A

    J. Bannister, A. Gibbs, and D. P. Hewett , Acoustic scattering by impedance screens/ cracks with fractal boundary: well-posedness analysis and boundary element approximation , Math. Models Methods Appl. Sci., 32 (2022), pp. 291–319

    work page 2022

  10. [10]

    Bannister, D

    J. Bannister, D. P. Hewett, and A. Gibbs , Scattering by fractal inhomogeneities via geometry-conforming Galerkin methods for the Lippmann-Schwinger equation . arXiv: 2602.05005

    work page arXiv

  11. [11]

    Barnsley and A

    M. Barnsley and A. Vince , Fractal tilings from iterated function systems , Discrete Comput. Geom., 51 (2014), pp. 729–752

    work page 2014

  12. [12]

    M. F. Barnsley and S. Demko , Iterated function systems and the global construction of fractals, Proc. Roy. Soc. Lond. A., 399 (1985), pp. 243–275

    work page 1985

  13. [13]

    Biegert, On traces of Sobolev functions on the boundary of extension domains , Proc

    M. Biegert, On traces of Sobolev functions on the boundary of extension domains , Proc. Am. Math. Soc., 137 (2009), pp. 4169–4176

    work page 2009

  14. [14]

    S. C. Brenner and L. R. Scott , The mathematical theory of finite element methods , Springer, 2008

    work page 2008

  15. [15]

    A. M. Caetano, S. N. Chandler-Wilde, X. Claeys, A. Gibbs, D. P. Hewett, and A. Moiola , Integral equation methods for acoustic scattering by fractals , Proc. R. Soc. A., (2025), p. 48120230650

    work page 2025

  16. [16]

    A. M. Caetano, S. N. Chandler-Wilde, A. Gibbs, D. P. Hewett, and A. Moiola , A Hausdorff-measure boundary element method for acoustic scattering by fractal screens , Numer. Math., 156 (2024), pp. 463–532

    work page 2024

  17. [17]

    A. M. Caetano, S. N. Chandler-Wilde, and D. P. Hewett , Properties of IFS attractors with non-empty interiors, related rough domains, and associated function spaces and scattering problems. arXiv:2511.15213

    work page arXiv

  18. [18]

    A. M. Caetano, D. P. Hewett, and A. Moiola , Density results for Sobolev, Besov and Triebel-Lizorkin spaces on rough sets , J. Funct. Anal., 281 (2021), p. 109019

    work page 2021

  19. [19]

    Cangiani, Z

    A. Cangiani, Z. Dong, and E. Georgoulis , hp-version discontinuous Galerkin methods on essentially arbitrarily-shaped elements , Math. Comp., 91 (2022), pp. 1–35

    work page 2022

  20. [20]

    Cangiani, Z

    A. Cangiani, Z. Dong, E. H. Georgoulis, and P. Houston , hp-Version Discontinuous Galerkin Methods on Polygonal and Polyhedral Meshes , Springer, 2017

    work page 2017

  21. [21]

    Capitanelli , Robin boundary condition on scale irregular fractals , Commun

    R. Capitanelli , Robin boundary condition on scale irregular fractals , Commun. Pure Appl. Anal., 9 (2010), pp. 1221–1234

    work page 2010

  22. [22]

    Capitanelli and M

    R. Capitanelli and M. A. Vivaldi , Weighted estimates on fractal domains , Mathematika, 61 (2015), pp. 370–384

    work page 2015

  23. [23]

    , Dynamical quasi-filling fractal layers , SIAM J. Math. Anal., 48 (2016), pp. 3931–3961

    work page 2016

  24. [24]

    Carvalho and A

    A. Carvalho and A. Caetano , On the Hausdorff dimension of continuous functions belong- ing to Hölder and Besov spaces on fractal d-sets , J. Fourier Anal. Appl., 18 (2012), pp. 386–409

    work page 2012

  25. [25]

    Castillo , Performance of discontinuous Galerkin methods for elliptic PDEs , SIAM J

    P. Castillo , Performance of discontinuous Galerkin methods for elliptic PDEs , SIAM J. Sci. Comput., 24 (2002), pp. 524–547. 43

    work page 2002

  26. [26]

    Cefalo, S

    M. Cefalo, S. Creo, M. Gallo, M. R. Lancia, and P. Vernole , Approximation of 3D Stokes flows in fractal domains , in Fractals in Engineering: Theoretical Aspects and Numerical Approximations, Springer, 2021, pp. 27–53

    work page 2021

  27. [27]

    Cefalo, S

    M. Cefalo, S. Creo, M. R. Lancia, and J. Rodríguez-Cuadrado , Fractal mixtures for optimal heat draining , Chaos, Solitons & Fractals, 173 (2023), p. 113750

    work page 2023

  28. [28]

    S. N. Chandler-Wilde and D. P. Hewett , Well-posed PDE and integral equation formu- lations for scattering by fractal screens , SIAM J. Math. Anal., 50 (2018), pp. 677–717

    work page 2018

  29. [29]

    S. N. Chandler-Wilde, D. P. Hewett, A. Moiola, and J. Besson , Boundary element methods for acoustic scattering by fractal screens , Numer. Math., 147 (2021), pp. 785–837

    work page 2021

  30. [30]

    P. G. Ciarlet , The Finite Element Method for Elliptic Problems , North-Holland, 1978

    work page 1978

  31. [31]

    Claret , Helmholtz transmission problem and intrinsic impedance scattering problem on extension domains

    G. Claret , Helmholtz transmission problem and intrinsic impedance scattering problem on extension domains . arXiv:2601.14866

    work page arXiv

  32. [32]

    Claret, A

    G. Claret, A. Rozanova-Pierrat, and A. Teplyaev , Convergence of layer potentials and Riemann-Hilbert problem on extension domains . arXiv:2411.03767

    work page arXiv

  33. [33]

    S. Creo, M. R. Lancia, P. Vernole, M. Hinz, and A. Teplyaev , Magnetostatic problems in fractal domains , in Analysis, Probability and Mathematical Physics on Fractals, World Scientific, 2020, pp. 477–502

    work page 2020

  34. [34]

    D. A. Di Pietro and A. Ern , Mathematical Aspects of Discontinuous Galerkin Methods , Springer, 2011

    work page 2011

  35. [35]

    F alconer, Fractal Geometry: Mathematical Foundations and Applications , Wiley, 3rd ed., 2014

    K. F alconer, Fractal Geometry: Mathematical Foundations and Applications , Wiley, 3rd ed., 2014

    work page 2014

  36. [36]

    F arkas and N

    W. F arkas and N. Jacob, Sobolev spaces on non smooth domains and Dirichlet forms related to subordinate reflecting diffusions , Math. Nachr., 224 (2001), pp. 75–104

    work page 2001

  37. [37]

    Gibbs, D

    A. Gibbs, D. Hewett, and A. Moiola , Numerical quadrature for singular integrals on fractals, Numer. Alg., 92 (2023), pp. 2071–2124

    work page 2023

  38. [38]

    Gibbs, D

    A. Gibbs, D. P. Hewett, and B. Major , Numerical evaluation of singular integrals on non-disjoint self-similar fractal sets , Numer. Alg., 97 (2024), pp. 311–343

    work page 2024

  39. [39]

    D. S. Grebenkov and B.-T. Nguyen , Geometrical structure of Laplacian eigenfunctions , SIAM Rev., 55 (2013), pp. 601–667

    work page 2013

  40. [40]

    Grisvard, Elliptic Problems in Nonsmooth Domains , Pitman, Boston, 1985

    P. Grisvard, Elliptic Problems in Nonsmooth Domains , Pitman, Boston, 1985

    work page 1985

  41. [41]

    Grochenig and A

    K. Grochenig and A. Haas , Self-similar lattice tilings , J. Fourier Anal. Appl., 1 (1994), pp. 131–170

    work page 1994

  42. [42]

    D. P. Hewett , Discontinuous piecewise polynomial approximation on non-Lipschitz domains . arXiv: 2511.22628

    work page arXiv

  43. [43]

    M. Hinz, S. N. Chandler-Wilde, and D. P. Hewett , Kernels of trace operators via fine continuity. arXiv:2507.04536. 44

    work page arXiv

  44. [44]

    M. Hinz, A. Rozanova-Pierrat, and A. Teplyaev , Boundary value problems on non- Lipschitz uniform domains: stability, compactness and the existence of optimal shapes , Asymp- totic Anal., 134 (2023), pp. 25–61

    work page 2023

  45. [45]

    J. E. Hutchinson , Fractals and self-similarity, Indiana Univ. Math. J., 30 (1981), pp. 713–747

    work page 1981

  46. [46]

    P. Joly, M. Kachanovska, and Z. Moitier , High-order numerical integration on self-affine sets. arXiv:2410.00637

    work page arXiv

  47. [47]

    P. W. Jones , Quasiconformal mappings and extendability of functions in Sobolev spaces , Acta Math., 147 (1981), pp. 71–88

    work page 1981

  48. [48]

    Keesling , The boundaries of self-similar tiles in Rn, Topol

    J. Keesling , The boundaries of self-similar tiles in Rn, Topol. Appl., 94 (1999), pp. 195–205

    work page 1999

  49. [49]

    Kozlov, V

    V. Kozlov, V. Maz’ya, and J. Rossmann , Elliptic boundary value problems in domains with point singularities , American Mathematical Soc., 1997

    work page 1997

  50. [50]

    M. R. Lancia , A transmission problem with a fractal interface , Z. Anal. Anwend., 21 (2002), pp. 113–133

    work page 2002

  51. [51]

    M. R. Lancia, M. Cefalo, and G. Dell Acqua , Numerical approximation of transmission problems across Koch-type highly conductive layers, Appl. Math. Comput., 218 (2012), pp. 5453– 5473

    work page 2012

  52. [52]

    M. R. Lancia and P. Vernole , Irregular heat flow problems , SIAM J. Math. Anal., 42 (2010), pp. 1539–1567

    work page 2010

  53. [53]

    , The Stokes problem in fractal domains: Asymptotic behaviour of the solutions , Discrete Contin. Dyn. Syst. - Ser. S, 13 (2020)

    work page 2020

  54. [54]

    M. L. Lapidus , Fractals and vibrations: can you hear the shape of a fractal drum? , Fractals, 3 (1995), pp. 725–736

    work page 1995

  55. [55]

    M. L. Lapidus, J. W. Neuberger, R. J. Renka, and C. A. Griffith , Snowflake har- monics and computer graphics: numerical computation of spectra on fractal drums , Int. J. Bifurcation Chaos, 6 (1996), pp. 1185–1210

    work page 1996

  56. [56]

    M. L. Lapidus and E. P. J. Pearse , A tube formula for the Koch snowflake curve, with applications to complex dimensions , J. Lond. Math. Soc., II. Ser., 74 (2006), pp. 397–414

    work page 2006

  57. [57]

    McLean , Strongly Elliptic Systems and Boundary Integral Equations , CUP, 2000

    W. McLean , Strongly Elliptic Systems and Boundary Integral Equations , CUP, 2000

    work page 2000

  58. [58]

    B. S. Mityagin , The zero set of a real analytic function , Math. Notes, 107 (2020), pp. 529–530

    work page 2020

  59. [59]

    J. M. Neuberger, N. Sieben, and J. W. Swift , Computing eigenfunctions on the Koch snowflake: a new grid and symmetry , J. Comput. Appl. Math., 191 (2006), pp. 126–142

    work page 2006

  60. [60]

    Nyström, Smoothness properties of solutions to Dirichlet problems in domains with a fractal boundary, PhD thesis, University of Umeå, 1994

    K. Nyström, Smoothness properties of solutions to Dirichlet problems in domains with a fractal boundary, PhD thesis, University of Umeå, 1994

    work page 1994

  61. [61]

    L. G. Rogers , Degree-independent Sobolev extension on locally uniform domains , J. Funct. Anal., 235 (2006), pp. 619–665

    work page 2006

  62. [62]

    Rösler and A

    F. Rösler and A. Stepanenko , Computing scattering resonances of rough obstacles . arXiv:2402.00846. 45

    work page arXiv

  63. [63]

    Comput., 93 (2024), pp

    , Computing eigenvalues of the Laplacian on rough domains , Math. Comput., 93 (2024), pp. 111–161

    work page 2024

  64. [64]

    Sapoval, T

    B. Sapoval, T. Gobron, and A. Margolina , Vibrations of fractal drums , Phys. Rev. Lett., 67 (1991), p. 2974

    work page 1991

  65. [65]

    Sayas, T

    F.-J. Sayas, T. Brown, and M. Hassell , Variational Techniques for Elliptic Partial Dif- ferential Equations , CRC Press, 2019

    work page 2019

  66. [66]

    R. S. Strichartz , Evaluating integrals using self-similarity , American Math. Monthly, 107 (2000), pp. 316–326

    work page 2000

  67. [67]

    R. S. Strichartz and Y. W ang , Geometry of self-affine tiles I , Indiana U. Math. J., (1999), pp. 1–23

    work page 1999

  68. [68]

    E. R. Vrscay and C. J. Roehrig , Iterated function systems and the inverse problem of fractal construction using moments , in Computers and Mathematics, Springer, 1989, pp. 250–259

    work page 1989

  69. [69]

    Ziemer , Weakly Differentiable Functions, Springer, 1989

    W. Ziemer , Weakly Differentiable Functions, Springer, 1989. 46

    work page 1989