pith. sign in

arxiv: 2606.03845 · v1 · pith:ZKPUQS4Bnew · submitted 2026-06-02 · 🧮 math.NA · cs.NA

Embedded Trefftz DG method for reaction-diffusion problems on anisotropic meshes

Sergio G\'omez , Chiara Perinati , Paul Stocker , Igor Voulis This is my paper

Pith reviewed 2026-06-28 08:48 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords discontinuous GalerkinTrefftz methodsanisotropic meshesreaction-diffusionquasi-optimalitya priori error estimatesembedded methodsquadrilateral elements
0
0 comments X

The pith

Embedding a relaxed local Trefftz condition into a tensor-product DG space produces a reduced discretization that remains stable and quasi-optimal for reaction-diffusion problems on anisotropic and curved quadrilateral meshes.

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

The paper constructs an embedded Trefftz discontinuous Galerkin method by imposing a relaxed local Trefftz condition through embedding into a tensor-product DG space. This construction reduces the size of the global linear system while retaining the approximation power of the underlying high-order discretization. The authors establish stability and quasi-optimality of the resulting method on anisotropic, possibly curved, quadrilateral elements and obtain corresponding anisotropic a priori error estimates. Numerical tests for both h- and hp-refinement on curved domains confirm that the theoretical guarantees hold in practice.

Core claim

By embedding a relaxed local Trefftz condition into a tensor-product DG space the method yields a smaller global system for reaction-diffusion problems that is stable and quasi-optimal on anisotropic, possibly curved, quadrilateral elements, together with anisotropic a priori error estimates that track the mesh stretching.

What carries the argument

The embedding into a tensor-product DG space that enforces a relaxed local Trefftz condition, thereby reducing degrees of freedom while preserving approximation properties on stretched and curved elements.

If this is right

  • The global linear system is smaller than that of the full tensor-product DG method while retaining the same approximation order.
  • Stability and quasi-optimality hold uniformly for both h- and hp-refinement on anisotropic and curved quadrilateral elements.
  • Anisotropic a priori error estimates are available that explicitly depend on the mesh aspect ratios.
  • The approach extends directly to curved-domain problems without loss of the theoretical guarantees.

Where Pith is reading between the lines

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

  • The same embedding technique could be tested on other elliptic operators where local solution spaces satisfy similar differential relations.
  • The reduction in system size may translate into lower memory use and faster iterative solves when the mesh anisotropy is extreme.
  • Extending the analysis to time-dependent reaction-diffusion problems would require checking whether the same embedding preserves stability in the time-stepping scheme.

Load-bearing premise

Imposing the relaxed local Trefftz condition via embedding into a tensor-product DG space preserves the approximation properties of the underlying high-order discretization on anisotropic and curved elements.

What would settle it

A sequence of numerical tests on increasingly stretched curved quadrilateral meshes in which the observed convergence rates fall below the predicted anisotropic estimates or the discrete system loses stability would refute the quasi-optimality claim.

Figures

Figures reproduced from arXiv: 2606.03845 by Chiara Perinati, Igor Voulis, Paul Stocker, Sergio G\'omez.

Figure 1
Figure 1. Figure 1: Construction of the elements in K via composition of an affine map FK and diffeomor￾phism GK. The map GK : Ke → K is a C 1 -diffeomorphism with Jacobian JGK : Ke → R 2×2 satisfying (C ⋄ G) −1 ≤ det JGK ≤ C ⋄ G and ∥JGK ∥L∞(Ke) ≤ C ∗ G ∀K ∈ K, (2.1) for some positive constants C ⋄ G and C ∗ G independent of K. 3 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: h-convergence for the problem with exact solution u in (4.1). Comparison between the standard DG method (Vh) and the embedded Trefftz DG method (Th) for polynomial degrees p = 3, 4, 5. This radially symmetric solution satisfies u = 0 on ∂Ω and exhibits a boundary layer near ∂Ω. To resolve this layer, we use curved anisotropic meshes refined towards ∂Ω. The mesh consists of a fixed triangular mesh of the in… view at source ↗
Figure 3
Figure 3. Figure 3: h-convergence for the problem with exact solution u in (4.2). Comparison between the standard DG method (Vh) and the embedded Trefftz method (Th) for polynomial degrees p = 4, 5 on two initial meshes with maximum element size h ≈ 0.15 (top row) and h ≈ 0.09 (bottom row). Refinement is performed anisotropically towards the boundary to resolve the boundary layers. The final mesh after eight refinement steps … view at source ↗
Figure 4
Figure 4. Figure 4: To resolve the boundary layers, we perform p-refinements while simultaneously adapting the size of the small elements adjacent to the boundary. The length of the small element edges is taken as ℓ = λpε, i.e., the increment of the polynomial degree p is accompanied by a proportional increase in the size of the small elements. We consider λ = 0.9, which yields good convergence properties as reported in [7], … view at source ↗
Figure 5
Figure 5. Figure 5: Numerical solution using the embedded Trefftz DG method with [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: hp-refinement for the problem (4.3) with exact solution u in (4.4). Left: finest mesh in the hp-refinement strategy, with geometric refinement towards the origin. Right: L 2 (Ω)-error versus square root of the number of degrees of freedom. Solid lines correspond to the hp-refinement strategy, while dashed lines correspond to p-refinement on the fixed finest mesh. We consider a sequence of three geometrical… view at source ↗
read the original abstract

We present and analyze an embedded Trefftz discontinuous Galerkin method for reaction-diffusion problems on anisotropic meshes. The method is constructed by imposing a relaxed local Trefftz condition via an embedding into a tensor-product DG space, yielding a reduced global system while preserving the approximation properties of the underlying high-order discretization. We prove stability and quasi-optimality on anisotropic, possibly curved, quadrilateral elements, and derive anisotropic a priori error estimates. Numerical experiments for $h$- and $hp$-refinement, including curved-domain examples, validate the theoretical results.

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 presents an embedded Trefftz discontinuous Galerkin method for reaction-diffusion problems on anisotropic meshes. The method is constructed by imposing a relaxed local Trefftz condition via embedding into a tensor-product DG space on anisotropic, possibly curved quadrilateral elements. This yields a reduced global system while preserving the approximation properties of the underlying high-order discretization. The paper proves stability and quasi-optimality, derives anisotropic a priori error estimates, and validates the results with numerical experiments for h- and hp-refinement, including curved-domain examples.

Significance. If the stability, quasi-optimality, and error estimates hold as claimed, the work provides a practical route to reduced-degree-of-freedom high-order DG discretizations that retain optimal approximation on anisotropic and curved meshes. This is relevant for efficient simulation of reaction-diffusion problems with boundary layers or complex geometries. The explicit construction via tensor-product embedding and the numerical validation on curved domains are positive features.

minor comments (2)
  1. The abstract states that the embedding 'preserves the approximation properties,' but a brief sentence in the introduction or §2 clarifying how the relaxed Trefftz condition interacts with the tensor-product space on curved elements would improve readability.
  2. Numerical experiments are mentioned for curved domains; adding a short remark on how the mesh curvature is handled in the implementation (e.g., via isoparametric mapping) would strengthen the validation section.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary of our work and the recommendation of minor revision. The referee's description accurately reflects the manuscript's contributions regarding the embedded Trefftz DG method, its stability and quasi-optimality analysis on anisotropic meshes, and the numerical validation.

Circularity Check

0 steps flagged

No circularity in derivation chain

full rationale

The paper constructs the embedded Trefftz DG method explicitly via relaxed local Trefftz condition in a tensor-product DG space, then applies standard DG stability and approximation arguments to prove quasi-optimality and anisotropic error estimates on curved quadrilaterals. No step reduces by definition to its own inputs, no fitted parameters are relabeled as predictions, and no load-bearing uniqueness or ansatz is imported solely via self-citation. The central claim (preservation of approximation properties) is analyzed rather than assumed tautologically, making the derivation self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central construction rests on the domain assumption that a relaxed Trefftz condition can be embedded without losing approximation power; no free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption A relaxed local Trefftz condition can be imposed via embedding into a tensor-product DG space while preserving approximation properties on anisotropic curved elements
    This is the key modeling choice enabling the reduced system and the subsequent stability and error analysis.

pith-pipeline@v0.9.1-grok · 5621 in / 1194 out tokens · 21836 ms · 2026-06-28T08:48:38.229458+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

29 extracted references · 3 canonical work pages

  1. [1]

    Apel.Anisotropic finite elements: local estimates and applications

    T. Apel.Anisotropic finite elements: local estimates and applications. Advances in Numerical Mathematics. B. G. Teubner, Stuttgart, 1999

    1999

  2. [2]

    Apel and M

    T. Apel and M. Dobrowolski. Anisotropic interpolation with applications to the finite element method.Computing, 47(3-4):277–293, 1992

    1992

  3. [3]

    Apel and G

    T. Apel and G. Lube. Anisotropic mesh refinement in stabilized Galerkin methods.Numer. Math., 74(3):261–282, 1996

    1996

  4. [4]

    Congreve, P

    S. Congreve, P. Houston, and I. Perugia. Adaptive refinement forhp-version Trefftz discon- tinuous Galerkin methods for the homogeneous Helmholtz problem.Adv. Comput. Math., 45(1):361–393, 2019. 20

    2019

  5. [5]

    Formaggia and S

    L. Formaggia and S. Perotto. Anisotropic error estimates for elliptic problems.Numer. Math., 94(1):67–92, 2003

    2003

  6. [6]

    E. H. Georgoulis.Discontinuous Galerkin methods on shape-regular and anisotropic meshes. PhD thesis, University of Oxford, 2003

    2003

  7. [7]

    E. H. Georgoulis.hp-version interior penalty discontinuous Galerkin finite element methods on anisotropic meshes.Int. J. Numer. Anal. Model., 3(1):52–79, 2006

    2006

  8. [8]

    E. H. Georgoulis, E. Hall, and P. Houston. Discontinuous Galerkin methods onhp-anisotropic meshes. I. A priori error analysis.Int. J. Comput. Sci. Math., 1(2-4):221–244, 2007

    2007

  9. [9]

    E. H. Georgoulis, E. Hall, and P. Houston. Discontinuous Galerkin methods for advection- diffusion-reaction problems on anisotropically refined meshes.SIAM J. Sci. Comput., 30(1):246–271, 2007/08

    2007

  10. [10]

    E. H. Georgoulis, E. Hall, and P. Houston. Discontinuous Galerkin methods onhp-anisotropic meshes. II. A posteriori error analysis and adaptivity.Appl. Numer. Math., 59(9):2179–2194, 2009

    2009

  11. [11]

    G´ omez, C

    S. G´ omez, C. Perinati, and P. Stocker. Inf-sup stable space-time local discontinuous Galerkin method for the heat equation.J. Sci. Comput., 106(1):Paper No. 22, 47, 2026

    2026

  12. [12]

    G´ omez, C

    S. G´ omez, C. Perinati, P. Stocker, and I. Voulis. Replication data for: Embedded Trefftz DG method for reaction-diffusion problems on anisotropic meshes, May 2026

    2026

  13. [13]

    C. Heil. Space-time Trefftz DG methods for parabolic PDEs. Master’s thesis, NAM, University of G¨ ottingen, March 2024

    2024

  14. [14]

    Heimann, C

    F. Heimann, C. Lehrenfeld, P. Stocker, and H. von Wahl. Unfitted Trefftz discontinuous Galerkin methods for elliptic boundary value problems.ESAIM Math. Model. Numer. Anal., 57(5):2803–2833, 2023

    2023

  15. [15]

    Hiptmair, A

    R. Hiptmair, A. Moiola, and I. Perugia. Trefftz discontinuous Galerkin methods for acoustic scattering on locally refined meshes.Appl. Numer. Math., 79:79–91, 2014

    2014

  16. [16]

    Hiptmair, A

    R. Hiptmair, A. Moiola, and I. Perugia. A survey of Trefftz methods for the Helmholtz equation. InBuilding bridges: connections and challenges in modern approaches to numerical partial differential equations, volume 114 ofLect. Notes Comput. Sci. Eng., pages 237–278. Springer, 2016

    2016

  17. [17]

    Imbert-G´ erard, A

    L.-M. Imbert-G´ erard, A. Moiola, C. Perinati, and P. Stocker. Polynomial quasi-Trefftz DG for PDEs with smooth coefficients: elliptic problems.IMA J. Numer. Anal., 45(6):3473–3506, 2025

    2025

  18. [18]

    Kretzschmar, S

    F. Kretzschmar, S. M. Schnepp, I. Tsukerman, and T. Weiland. Discontinuous Galerkin methods with Trefftz approximations.J. Comput. Appl. Math., 270:211–222, 2014

    2014

  19. [19]

    P. L. Lederer, C. Lehrenfeld, and P. Stocker. Trefftz discontinuous Galerkin discretization for the Stokes problem.Numer. Math., 156(3):979–1013, 2024

    2024

  20. [20]

    P. L. Lederer, C. Lehrenfeld, P. Stocker, and I. Voulis. Embedded Trefftz DG framework for the analysis of discretizations with local-global decompositions. arXiv:2412.00806, 2024

    work page arXiv 2024

  21. [21]

    Lehrenfeld and P

    C. Lehrenfeld and P. Stocker. Embedded Trefftz discontinuous Galerkin methods.Internat. J. Numer. Methods Engrg., 124(17):3637–3661, 2023

    2023

  22. [22]

    Lehrenfeld, P

    C. Lehrenfeld, P. Stocker, and M. Zienecker. Sparsity comparison of polytopal finite element methods.PAMM, 24(3):e202400150, 2024. 21

    2024

  23. [23]

    Lozinski

    A. Lozinski. A primal discontinuous Galerkin method with static condensation on very general meshes.Numer. Math., 143(3):583–604, 2019

    2019

  24. [24]

    Perinati

    C. Perinati. A quasi-Trefftz discontinuous Galerkin method for the homogeneous diffusion- advection-reaction equation with piecewise-smooth coefficients. arXiv:2312.09919, 2023

    work page arXiv 2023

  25. [25]

    Schlesinger

    E. Schlesinger. Embedded Trefftz Trace DG Methods for PDEs on unfitted Surfaces. Master’s thesis, NAM, University of G¨ ottingen, April 2023

    2023

  26. [26]

    Sch¨ oberl

    J. Sch¨ oberl. C++ 11 implementation of finite elements in NGSolve.Institute for Analysis and Scientific Computing, Vienna University of Technology, 30, 2014

    2014

  27. [27]

    Schwab and M

    C. Schwab and M. Suri. Thepandhpversions of the finite element method for problems with boundary layers.Math. Comp., 65(216):1403–1429, 1996

    1996

  28. [28]

    P. Stocker. NGSTrefftz: Add-on to NGSolve for Trefftz methods.J. Open Source Softw., 7(71):4135, 2022

    2022

  29. [29]

    Stocker and I

    P. Stocker and I. Voulis. Embedded Trefftz DG method for the Helmholtz equation. arXiv.2603.13034, 2026. A Appendix The following lemma allows us to verify the assumption (3.3) via a perturbation argument. Lemma A.1([20, Lemma 4.2]).LetK∈ Kand letA K,0 :L h(K)→L h(K) ′ be an invertible operator with ∥AK,0uh∥Lh(K) ′ ≥c A0 ∥uh∥Lh(K) ∀uh ∈L h(K), for some co...

    work page arXiv 2026