pith. sign in

arxiv: 2604.18175 · v1 · submitted 2026-04-20 · 🧮 math.NA · cs.NA

Trefftz methods with evanescent plane waves

Andrea Moiola , Nicola Galante , Emile Parolin This is my paper

Pith reviewed 2026-05-10 04:03 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Trefftz methodsevanescent plane wavesHelmholtz equationUltraweak Variational Formulationnumerical instabilitiesplane wave bases
0
0 comments X

The pith

Evanescent plane wave bases substantially mitigate numerical instabilities in classical Trefftz methods for the Helmholtz equation.

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

Classical Trefftz methods approximate solutions to the Helmholtz equation using propagative plane waves, but they often encounter strong numerical instabilities. The authors propose that incorporating evanescent plane waves, chosen with a simple recipe, can address this issue effectively. They demonstrate that applying this approach in the Ultraweak Variational Formulation leads to greatly improved numerical results. This matters because it provides a straightforward enhancement to existing Trefftz frameworks without requiring complex modifications. Readers interested in numerical methods for wave problems would find this a practical step toward more reliable simulations.

Core claim

Evanescent plane wave bases can substantially mitigate numerical instabilities in classical Trefftz methods. A simple recipe is proposed to select such basis functions. Numerical results for the Ultraweak Variational Formulation improve greatly with this choice.

What carries the argument

Evanescent plane wave basis functions selected by a proposed simple recipe, used within the Ultraweak Variational Formulation of Trefftz methods.

Load-bearing premise

The simple recipe for selecting evanescent plane waves will produce consistent improvements for the UWVF and other Trefftz formulations beyond the limited cases referenced.

What would settle it

Numerical experiments on additional Helmholtz problems or other Trefftz variants where the evanescent wave selection fails to reduce instabilities or improve accuracy.

Figures

Figures reproduced from arXiv: 2604.18175 by Andrea Moiola, Emile Parolin, Nicola Galante.

Figure 1
Figure 1. Figure 1: A PPW and an EPW on the disc in R 2 , an EPW on the ball in R 3 . Why do we expect linear combinations of the EPWs (3) to give accurate small-coefficient ap￾proximations? The reason is in [11, Thm. 6.7] (2D) and [5, Thm. 3.9] (3D): any Helmholtz so￾lution u ∈ H1 (B), with B the unit disc/ball, is a continuous EPW superposition. The coefficient of this superposition, measured in a weighted L 2 norm, is boun… view at source ↗
Figure 3
Figure 3. Figure 3: Fundamental solution: convergence result for increasing κ using an approximation space of size increasing linearly in κ. mesh contains 64 triangles. Since the solution to this problem is not known analytically, we use as reference solution the EPW approximation with twice as many waves as the largest approxima￾tion space. This test case is more challenging because of the presence of the cavity and the corn… view at source ↗
Figure 2
Figure 2. Figure 2: Fundamental solution: ℜ(u) (top), pointwise error (middle) for PPWs and EPWs, convergence history (bottom). 103 104 Number of degrees of freedom N trial 100 101 102 10−8 10−6 10−4 10−2 100 Wavenumber κ Relative κ weighted H1-error PPW EPW [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: Scattering problem: ℜ(u) (top), pointwise error (middle) for PPWs and EPWs, convergence history (bottom). Helmholtz equation. SINUM 35.1 (1998), 255–299. [4] J. Coyle, and N. Nigam, The whys and hows of conditioning of DG plane wave Tr￾efftz methods: a single element. arXiv:2509.14500 (2025). [5] N. Galante, A. Moiola, and E. Parolin, Stable approximation of Helmholtz solu￾tions in the 3D ball using evanes… view at source ↗
read the original abstract

Classical Trefftz methods approximate Helmholtz solutions using propagative plane waves and are subject to strong numerical instabilities. Evanescent plane wave bases can substantially mitigate this phenomenon. We propose a simple recipe to select such basis functions. We show that the numerical results obtained by the Ultraweak Variational Formulation (UWVF) greatly improve thanks to this choice. More details and examples will soon be available in [Galante, Moiola, Parolin 2026].

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

4 major / 0 minor

Summary. The manuscript claims that classical Trefftz methods for the Helmholtz equation suffer from numerical instabilities when using only propagative plane waves, that evanescent plane wave bases can substantially mitigate this issue, and that a simple recipe for selecting such bases yields greatly improved numerical results when applied to the Ultraweak Variational Formulation (UWVF). Further details and examples are deferred to a forthcoming publication.

Significance. If the asserted stabilization effect holds with supporting analysis and experiments, the work could address a recognized practical limitation of plane-wave Trefftz discretizations at moderate-to-high frequencies. At present, however, the manuscript contains no derivations, no numerical data, no error tables, and no stability estimates, so any potential significance remains entirely speculative.

major comments (4)
  1. [Abstract] The manuscript supplies neither the explicit form of the proposed 'simple recipe' for selecting evanescent plane-wave directions nor any equations defining the basis functions or their incorporation into the UWVF.
  2. [Abstract] No numerical results, convergence tables, condition-number plots, or comparisons against standard propagative bases or other stabilizations are provided to substantiate the claim that the UWVF 'greatly improve[s]' thanks to the choice.
  3. [Abstract] The central assertion that evanescent waves mitigate instabilities rests solely on an unproven statement; the text contains no a-priori bounds, no analysis of conditioning under mesh refinement or frequency increase, and no discussion of the heuristic assumptions underlying the selection recipe.
  4. [Abstract] The manuscript explicitly defers all details and examples to the forthcoming reference [Galante, Moiola, Parolin 2026], rendering the present submission non-self-contained and unsuitable for independent review.

Simulated Author's Rebuttal

4 responses · 1 unresolved

We appreciate the referee's detailed comments on our manuscript. We recognize that the submission is concise and primarily serves to introduce the concept, with full details deferred to a forthcoming publication. We address each major comment below.

read point-by-point responses

  1. Referee: [Abstract] The manuscript supplies neither the explicit form of the proposed 'simple recipe' for selecting evanescent plane-wave directions nor any equations defining the basis functions or their incorporation into the UWVF.

    Authors: We agree that the explicit form of the recipe and the defining equations are not provided in the current manuscript. These are detailed in the forthcoming paper [Galante, Moiola, Parolin 2026]. The present submission is intended as a brief announcement of the stabilization effect observed with evanescent plane waves in Trefftz methods for the Helmholtz equation. revision: no

  2. Referee: [Abstract] No numerical results, convergence tables, condition-number plots, or comparisons against standard propagative bases or other stabilizations are provided to substantiate the claim that the UWVF 'greatly improve[s]' thanks to the choice.

    Authors: The numerical results, tables, plots, and comparisons are indeed not included here, as they form part of the forthcoming publication. The claim in the abstract is based on those results, which demonstrate the improvement in the UWVF performance. revision: no

  3. Referee: [Abstract] The central assertion that evanescent waves mitigate instabilities rests solely on an unproven statement; the text contains no a-priori bounds, no analysis of conditioning under mesh refinement or frequency increase, and no discussion of the heuristic assumptions underlying the selection recipe.

    Authors: The manuscript does not provide a-priori bounds or detailed analysis, as these are developed in [Galante, Moiola, Parolin 2026]. The assertion is based on numerical evidence from that work, and the selection recipe is heuristic as stated. revision: no

  4. Referee: [Abstract] The manuscript explicitly defers all details and examples to the forthcoming reference [Galante, Moiola, Parolin 2026], rendering the present submission non-self-contained and unsuitable for independent review.

    Authors: We acknowledge that the deferral to the forthcoming reference limits the self-containment of this submission. This format is chosen to quickly communicate the main idea and its potential impact on Trefftz methods. We are open to expanding the manuscript if the journal permits, but the core contribution is the proposal of the simple recipe for evanescent waves. revision: partial

standing simulated objections not resolved
  • The explicit recipe, numerical results, convergence data, condition number analysis, and stability estimates, which are all deferred to the forthcoming publication [Galante, Moiola, Parolin 2026] and not present in the current manuscript.

Circularity Check

0 steps flagged

No circularity: proposal and numerical evidence are independent of inputs

full rationale

The paper proposes a simple recipe for choosing evanescent plane-wave directions and reports numerical improvement in UWVF condition numbers and accuracy. No derivation chain, first-principles prediction, or fitted parameter is presented in the abstract or described content. The central claim rests on numerical illustrations rather than any self-definitional relation, renamed known result, or load-bearing self-citation. The single forward self-citation to forthcoming work is not used to justify the recipe or its stability benefit. The manuscript is therefore self-contained against external benchmarks for the purpose of circularity analysis.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no information on free parameters, axioms, or invented entities; ledger is therefore empty.

pith-pipeline@v0.9.0 · 5362 in / 946 out tokens · 29795 ms · 2026-05-10T04:03:29.796179+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

12 extracted references · 12 canonical work pages

  1. [1]

    Adcock, and D

    B. Adcock, and D. Huybrechs, Frames and numerical approximation.SIAM Rev.61 (2019), 443–473

    work page 2019

  2. [2]

    Barucq, A

    H. Barucq, A. Bendali, J. Diaz, and S. Tordeux, Local strategies for improving the conditioning of the plane-wave ultra-weak variational formulation.J. Comput. Phys. 441(2021), 110449, 18

    work page 2021

  3. [3]

    Cessenat, and B

    O. Cessenat, and B. Després, Application of an ultra weak variational formulation of elliptic PDEs to the two-dimensional WAVES 2026, Montreal, Canada 6 Real part of the solution ℜ(u) −2 −1 0 1 2 Pointwise error for PPW 10−8 10−7 10−6 10−5 10−4 10−3 10−2 10−1 Pointwise error for EPW 10−8 10−7 10−6 10−5 10−4 10−3 10−2 10−1 103 104 105 10−3 10−2 10−1 100 Nu...

    work page 2026

  4. [4]

    Coyle, and N

    J. Coyle, and N. Nigam, The whys and hows of conditioning of DG plane wave Tr- efftz methods: a single element. arXiv:2509.14500(2025)

    work page arXiv 2025

  5. [5]

    Galante, A

    N. Galante, A. Moiola, and E. Parolin, Stable approximation of Helmholtz solu- tions in the 3D ball using evanescent plane waves.SMAI J. Comput. Math.11(2025), 435–472

    work page 2025

  6. [6]

    Galante, A

    N. Galante, A. Moiola, and E. Parolin, Trefftz Discontinuous Galerkin method us- ing evanescent plane waves. In preparation (2026)

    work page 2026

  7. [7]

    Hiptmair, A

    R. Hiptmair, A. Moiola, and I. Perugia, Plane wave discontinuous Galerkin meth- ods for the 2D Helmholtz equation: anal- ysis of the p-version.SINUM49(2011), 264–284

    work page 2011

  8. [8]

    Hiptmair, A

    R. Hiptmair, A. Moiola, and I. Perugia, A surveyofTrefftzmethodsfortheHelmholtz equation,Lect. Notes Comput. Sci. Eng. (2016), 237–278

    work page 2016

  9. [9]

    Huttunen, P

    T. Huttunen, P. Monk, and J. P. Kaipio, Computational aspects of the ultra-weak variational formulation.J. Comput. Phys. 182.1(2002), 27–46

    work page 2002

  10. [10]

    Moiola, R

    A. Moiola, R. Hiptmair, and I. Perugia, Plane wave approximation of homogeneous Helmholtz solutions.Z.Angew. Math. Phys. 62(2011), 809–837

    work page 2011

  11. [11]

    Parolin, D

    E. Parolin, D. Huybrechs, and A. Moiola, Stable approximation of Helmholtz solu- tions in the disk by evanescent plane waves. M2AN57(6)(2023), 3499–3536

    work page 2023

  12. [12]

    Trefftz, Ein Gegenstuck zum Ritzschen Verfahren,Proceedings of the 2nd Inter- national Congress for Applied Mechanics, Zurich (1926), 131–137

    E. Trefftz, Ein Gegenstuck zum Ritzschen Verfahren,Proceedings of the 2nd Inter- national Congress for Applied Mechanics, Zurich (1926), 131–137

    work page 1926