pith. sign in

arxiv: 1907.08309 · v1 · pith:WG57OFB5new · submitted 2019-07-18 · 🧮 math.NA · cs.NA

A roadmap for Generalized Plane Waves and their interpolation properties

Lise-Marie Imbert-Gerard , Guillaume Sylvand This is my paper

Pith reviewed 2026-05-24 19:17 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Generalized Plane WavesDiscontinuous GalerkinTrefftz methodsvariable coefficientsinterpolation propertieswave equationsnumerical PDEhigh-order approximation
0
0 comments X

The pith

A step-by-step algorithm constructs Generalized Plane Waves that satisfy necessary conditions for high-order interpolation on variable-coefficient wave problems.

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

The paper develops a general roadmap for building Generalized Plane Waves as approximate local solutions to PDEs with variable coefficients. These functions serve as basis elements in Discontinuous Galerkin schemes for wave problems where exact Trefftz solutions do not exist. The authors isolate each construction step, supply an algorithm that assembles the functions, and derive algebraic conditions required for the basis to reproduce polynomials at high order. A reader would care because the conditions make it possible to retain the accuracy advantages of oscillating bases even when material properties change in space.

Core claim

By systematically approximating the governing PDE locally and then building the Generalized Plane Waves from that approximation, one obtains functions whose interpolation properties meet explicit necessary conditions for high-order polynomial reproduction; the same procedure applies to families of wave equations beyond the constant-coefficient Helmholtz case.

What carries the argument

The algorithm that summarizes the construction of Generalized Plane Waves by sequencing local PDE approximation, coefficient expansion, and algebraic enforcement of interpolation conditions.

If this is right

  • The same construction steps extend Trefftz-type bases to a wider class of variable-coefficient wave equations.
  • The derived algebraic conditions directly control the polynomial reproduction order of the resulting DG basis.
  • The roadmap supplies a repeatable procedure that can be applied to new PDEs without starting from scratch.
  • High-order interpolation follows once the local PDE approximation and the algebraic steps are both satisfied.

Where Pith is reading between the lines

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

  • The procedure could be tested on concrete variable-coefficient acoustic or elastic problems to measure actual convergence rates.
  • Similar algebraic conditions might appear in other non-polynomial bases that incorporate PDE information.
  • If the local approximation error can be bounded a priori, the method would give explicit a-priori error estimates for the full DG scheme.

Load-bearing premise

The local approximation to the governing PDE stays accurate enough that the resulting functions still obey the algebraic relations needed for high-order polynomial reproduction.

What would settle it

Numerical tests in which the observed convergence order for the interpolation error drops below the predicted rate as the spatial variation of the coefficients is increased would show the necessary conditions do not hold.

Figures

Figures reproduced from arXiv: 1907.08309 by Guillaume Sylvand, Lise-Marie Imbert-Gerard.

Figure 1
Figure 1. Figure 1: Representation of the indices involved in the nonlinear system (10), for [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Representation of the indices of equations and unknowns from the initial nonlinear system (10) [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Representation of the indices of unknowns involved in two equations ( [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: GPW approximation of uAd by ua ∈ V 0 α,p,q with p = 2n + 1. We represent the L∞ error max(x0,y0)∈Ω kuAd − uakL∞({(x,y)∈R2,|(x,y)−(x0,y0)|<h}) , for 50 random points (x0, y0) ∈ ΩAd. We compare results for parameters satisfying Theorem 1 hypotheses i.e. q = max(1, n − 1) (Left panel), and for varying q with fixed n (Right panel). 100 10−2 10−4 10−6 10−8 10−16 10−12 10−8 10−4 100 h max error on disks of radiu… view at source ↗
Figure 5
Figure 5. Figure 5: GPW approximation of uJc by ua ∈ V 0 α,p,q with p = 2n + 1. We represent the L∞ error max(x0,y0)∈Ω kuJc − uakL∞({(x,y)∈R2,|(x,y)−(x0,y0)|<h}) , for 50 random points (x0, y0) ∈ ΩJc. We compare results for parameters satisfying Theorem 1 hypotheses i.e. q = max(1, n − 1) (Left panel), and for varying q with fixed n (Right panel). 33 [PITH_FULL_IMAGE:figures/full_fig_p033_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: GPW approximation of uJJ by ua ∈ V 0 α,p,q with p = 2n + 1. We represent the L∞ error max(x0,y0)∈Ω kuJJ − uakL∞({(x,y)∈R2,|(x,y)−(x0,y0)|<h}) , for 50 random points (x0, y0) ∈ ΩJJ . We compare results for parameters satisfying Theorem 1 hypotheses i.e. q = max(1, n − 1) (Left panel), and for varying q with fixed n (Right panel). 100 10−2 10−4 10−6 10−8 10−16 10−12 10−8 10−4 100 h max error on disks of radi… view at source ↗
Figure 7
Figure 7. Figure 7: GPW approximation of ucs by ua ∈ V 0 α,p,q with p = 2n + 1. We represent the L∞ error max(x0,y0)∈Ω kucs−uakL∞({(x,y)∈R2,|(x,y)−(x0,y0)|<h}) , for 50 random points (x0, y0) ∈ Ωcs. We compare results for parameters satisfying Theorem 1 hypotheses i.e. q = max(1, n − 1) (Left panel), and for varying q with fixed n (Right panel). 34 [PITH_FULL_IMAGE:figures/full_fig_p034_7.png] view at source ↗
read the original abstract

This work focuses on the study of partial differential equation (PDE) based basis function for Discontinuous Galerkin methods to solve numerically wave-related boundary value problems with variable coefficients. To tackle problems with constant coefficients, wave-based methods have been widely studied in the literature: they rely on the concept of Trefftz functions, i.e. local solutions to the governing PDE, using oscillating basis functions rather than polynomial functions to represent the numerical solution. Generalized Plane Waves (GPWs) are an alternative developed to tackle problems with variable coefficients, in which case Trefftz functions are not available. In a similar way, they incorporate information on the PDE, however they are only approximate Trefftz functions since they don't solve the governing PDE exactly, but only an approximated PDE. Considering a new set of PDEs beyond the Helmholtz equation, we propose to set a roadmap for the construction and study of local interpolation properties of GPWs. Identifying carefully the various steps of the process, we provide an algorithm to summarize the construction of these functions, and establish necessary conditions to obtain high order interpolation properties of the corresponding basis.

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

2 major / 2 minor

Summary. The paper develops a roadmap for constructing Generalized Plane Waves (GPWs) as approximate Trefftz functions for Discontinuous Galerkin discretizations of wave problems with variable coefficients. It supplies an algorithm that organizes the construction steps for PDEs beyond the Helmholtz equation and states necessary conditions under which the resulting basis functions achieve high-order interpolation properties.

Significance. If the stated necessary conditions remain valid once the local PDE approximation error is taken into account, the algorithm would give a reproducible procedure for building high-order, PDE-informed bases in the variable-coefficient setting where exact Trefftz functions do not exist. The explicit algorithmic summary is a concrete strength that could be directly implemented and tested.

major comments (2)
  1. [paragraph after the algorithm for GPW construction] The necessary conditions for high-order polynomial reproduction are derived under the assumption that the GPWs satisfy an exactly approximated PDE; the manuscript does not supply an a-priori bound showing that the perturbation introduced by the variable-coefficient approximation preserves those algebraic relations (see the paragraph following the algorithm statement).
  2. [section establishing necessary conditions for interpolation properties] The claim that the constructed basis attains a prescribed interpolation order therefore rests on an unverified continuity argument between the constant-coefficient case and the variable-coefficient approximation; without this step the necessary conditions do not yet imply the stated interpolation result.
minor comments (2)
  1. Notation for the local approximating operator and the precise definition of the residual PDE should be introduced once and used consistently.
  2. The abstract lists “a new set of PDEs” but does not name them; an explicit list in the introduction would help readers assess the scope.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for highlighting the distinction between the approximated PDE and the original variable-coefficient problem. We address each major comment below and will revise the manuscript accordingly to clarify the scope of the necessary conditions.

read point-by-point responses

  1. Referee: [paragraph after the algorithm for GPW construction] The necessary conditions for high-order polynomial reproduction are derived under the assumption that the GPWs satisfy an exactly approximated PDE; the manuscript does not supply an a-priori bound showing that the perturbation introduced by the variable-coefficient approximation preserves those algebraic relations (see the paragraph following the algorithm statement).

    Authors: The construction algorithm is designed so that each GPW satisfies the chosen approximated PDE exactly; the necessary conditions for polynomial reproduction are therefore algebraic identities that hold exactly in that setting. The variable-coefficient approximation enters only through the definition of the approximated PDE itself. We agree that the manuscript would be strengthened by an explicit remark on how this choice of approximation affects the relations when the original PDE is considered. In the revised version we will expand the paragraph following the algorithm to include such a remark, making clear that the algebraic conditions apply to the approximated PDE and that controlling the perturbation for the original PDE belongs to the subsequent global error analysis. revision: yes

  2. Referee: [section establishing necessary conditions for interpolation properties] The claim that the constructed basis attains a prescribed interpolation order therefore rests on an unverified continuity argument between the constant-coefficient case and the variable-coefficient approximation; without this step the necessary conditions do not yet imply the stated interpolation result.

    Authors: The manuscript states necessary conditions under which high-order interpolation properties hold when the GPWs solve the approximated PDE exactly. It does not contain a proof that the basis attains the order for the original variable-coefficient PDE, nor does it invoke a continuity argument between the constant- and variable-coefficient cases. If any phrasing in the section on necessary conditions could be read as asserting attainment without that step, we will revise the wording to emphasize that the stated conditions are necessary (and sufficient) only for the approximated PDE, leaving the passage from the approximation to the original PDE for future work or for the DG convergence analysis. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation supplies algorithm and necessary conditions without reducing to fitted inputs or self-citations

full rationale

The paper presents a roadmap and algorithm for constructing Generalized Plane Waves as approximate Trefftz functions for variable-coefficient PDEs, then states necessary conditions for high-order interpolation properties. No equations, fitted parameters, or predictions appear that reduce by construction to the inputs; the central steps involve identifying construction steps and algebraic conditions rather than renaming or self-referentially justifying results. The provided abstract and context contain no self-citation load-bearing claims or ansatz smuggling, making the derivation self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no explicit free parameters, axioms, or invented entities; the ledger is therefore empty.

pith-pipeline@v0.9.0 · 5722 in / 1103 out tokens · 17403 ms · 2026-05-24T19:17:09.747845+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

34 extracted references · 34 canonical work pages

  1. [1]

    Internat

    Babuska I., Melenk J.M., The partition of unity method. Internat. J. Numer. Methods Engrg. 40 (1997), no. 4, pp. 727–758

    work page 1997

  2. [2]

    Symposium on Advances in Computational Mechanics, Vol

    Babuska I., Zhang Z., The partition of unity method for the elastically supported beam. Symposium on Advances in Computational Mechanics, Vol. 5. Comput. Methods Appl. Mech. Engrg. 152 (1998), no. 1-2, pp. 1–18

    work page 1998

  3. [3]

    Buet C., Despres B., Morel G., Trefftz Discontinuous Galerkin basis functions for a class of Friedrichs systems coming from linear transport, hal-01964528. 37

    work page

  4. [4]

    Probl` emes de Helmholtz 2D et de Maxwell 3D, Universit´ e Paris 9 Dauphine, (1996)

    Cessenat O., Application d’une nouvelle formulation variationnelle aux ´ equations d’ondes harmoniques. Probl` emes de Helmholtz 2D et de Maxwell 3D, Universit´ e Paris 9 Dauphine, (1996)

    work page 1996

  5. [5]

    M.; Savits, T

    Constantine, G. M.; Savits, T. H. A multivariate Fa` a di Bruno formula with applications.Trans. Amer. Math. Soc. 348 no. 2 (1996), pp. 503–520

    work page 1996

  6. [6]

    Craik, A. D. D. Prehistory of Faa di Bruno’s formula. Amer. Math. Monthly 112 no. 2 (2005), pp. 119– 130

    work page 2005

  7. [7]

    , The propagation of gravity waves from deep to shallow water, Circular 20, National Bureau of Standards (1952), pp

    Eckart C. , The propagation of gravity waves from deep to shallow water, Circular 20, National Bureau of Standards (1952), pp. 165–173

    work page 1952

  8. [8]

    P., The discontinuous enrichment method

    Farhat C., Harari I., Franca L. P., The discontinuous enrichment method. Comput. Methods Appl. Mech. Engrg. 190 (2001), no. 48, pp. 6455–6479

    work page 2001

  9. [9]

    Fix G. J. , Gulati S. Wakoff G. I., On the use of singular functions with finite element approximations. J. Computational Phys. 13 (1973), 209–228

    work page 1973

  10. [10]

    J., Hiptmair R., Dispersion analysis of plane wave discontinuous Galerkin methods, Inter- national Journal for Numerical Methods in Engineering, 98 (2014) 5, pp

    Gittelson C. J., Hiptmair R., Dispersion analysis of plane wave discontinuous Galerkin methods, Inter- national Journal for Numerical Methods in Engineering, 98 (2014) 5, pp. 313–323

    work page 2014

  11. [11]

    Electron

    Hardy M., Combinatorics of partial derivatives. Electron. J. Comb. 13(1), 13 (2006)

    work page 2006

  12. [12]

    P., Computational Aspects of the Ultra-Weak Variational Formulation, Journal of Computational Physics, 182 (2002) 1, pp

    Huttunen T., Monk P., Kaipio J. P., Computational Aspects of the Ultra-Weak Variational Formulation, Journal of Computational Physics, 182 (2002) 1, pp. 27–46

    work page 2002

  13. [13]

    Imbert-G´ erard L.-M., Interpolation properties of generalized plane waves, Numer. Math. (2015) 131:683–711

    work page 2015

  14. [14]

    Imbert-G´ erard, L.-M., Generalized plane waves for varying coefficients, Proceedings of Waves 2015, Karslruhe, Germany

    work page 2015

  15. [15]

    Imbert-Gerard, L.-M., Despres, B.: A generalized plane-wave numerical method for smooth noncon- stant coefficients. IMA J. Numer. Anal. (2013). doi:10.1093/imanum/drt030

    work page doi:10.1093/imanum/drt030 2013

  16. [16]

    1387–1406

    Imbert-G´ erard, L.-M., Monk P., Numerical simulation of wave propagation in inhomogeneous media using Generalized Plane Waves, ESAIM: M2AN, 51 4 (2017), pp. 1387–1406

    work page 2017

  17. [17]

    Imbert-G´ erard L.-M., Well-posedness and generalized plane waves simulations of a 2D mode conversion model, Journal of Computational Physics, 303 (2015) pp. 105–124

    work page 2015

  18. [18]

    Kita E., Kamiya N., Trefftz method: an overview, Advances in Engineering Software 24 (1995) pp. 3–12

    work page 1995

  19. [19]

    Kretzschmar F., Moiola A., Perugia I., Schnepp S.M., A priori error analysis of space?time Trefftz discontinuous Galerkin methods for wave problems, IMA Journal of NumericalAnalysis, 36 (2016), p. 1599

    work page 2016

  20. [20]

    Kretzschmar F., Schnepp S.M., Tsukerman I., Weiland T., Discontinuous Galerkin methods with Trefftz approximations, Journal of Computational and Applied Mathematics, 270 (2014), pp. 211–222

    work page 2014

  21. [21]

    Kupradze V.D., Gegelia T.G., Basheleishvili M.O., Burchuladze T.V., Three?Dimensional Problems of the Mathematical Theory of Elasticity and Thermoelasticity, North?Holland, Amsterdam, New York, Oxford (1979)

    work page 1979

  22. [22]

    237–278, Lect

    Hiptmair R., Moiola A., Perugia I., A Survey of Trefftz Methods for the Helmholtz Equation, Building bridges: connections and challenges in modern approaches to numerical partial differential equations, pp. 237–278, Lect. Notes Comput. Sci. Eng., 114, Springer, 2016. 38

    work page 2016

  23. [23]

    Gabard G., B´ eriot H., A comparison of high-order polynomial and wave-based methods for Helmholtz problems

    Lieu A. Gabard G., B´ eriot H., A comparison of high-order polynomial and wave-based methods for Helmholtz problems. J. Comput. Phys. 321 (2016), pp. 105–125

    work page 2016

  24. [24]

    Luostari T., Huttunen T., Monk P., The ultra weak variational formulation using Bessel basis functions. Commun. Comput. Phys. 11 no. 2 (2012), pp. 400–414

    work page 2012

  25. [25]

    Higher chain formula proved by combinatorics

    Ma, T.-W. Higher chain formula proved by combinatorics. Electron. J. Combin. 16 no. 1 (2009), 7 pp

    work page 2009

  26. [26]

    Maunder E.A.W., Trefftz in translation, Computer Assisted Mechanics and Engineering Sciences 10 (2003)

    work page 2003

  27. [27]

    M., On generalized finite element methods

    Melenk J. M., On generalized finite element methods. PhD thesis, The University of Maryland, 1995

    work page 1995

  28. [28]

    M., Babuska I., The partition of unity finite element method: basic theory and applications

    Melenk J. M., Babuska I., The partition of unity finite element method: basic theory and applications. Comput. Methods Appl. Mech. Engrg. 139 (1996), no. 1-4, pp. 289–314

    work page 1996

  29. [29]

    Pergamon Press; distributed by Macmillan, New York (1964)

    Mikhlin S.G., Variational methods in mathematical physics. Pergamon Press; distributed by Macmillan, New York (1964)

    work page 1964

  30. [30]

    , Despres B., Trefftz Discontinuous Galerkin Method for Friedrichs Systems with Linear Relaxation: Application to the P 1 Model

    Morel G., Buet C. , Despres B., Trefftz Discontinuous Galerkin Method for Friedrichs Systems with Linear Relaxation: Application to the P 1 Model. Computational Methods in Applied Mathematics, 18 no. 3 (2018), pp. 521–557

    work page 2018

  31. [31]

    , Perugia I., A space-time Trefftz discontinuous Galerkin method for the acoustic wave equation in first-order formulation, Numer

    Moiola A. , Perugia I., A space-time Trefftz discontinuous Galerkin method for the acoustic wave equation in first-order formulation, Numer. Math., 138 no. 2 (2018), pp. 389–435

    work page 2018

  32. [32]

    Springer Science & Business Media (2012)

    Rektorys K., Variational methods in mathematics, science and engineering. Springer Science & Business Media (2012)

    work page 2012

  33. [33]

    J., An analysis of the finite element method

    Strang G., Fix G. J., An analysis of the finite element method. Prentice-Hall Series in Automatic Computation. Prentice-Hall, Inc., Englewood Cliffs, N. J., 1973

    work page 1973

  34. [34]

    In Proceedings of the 2nd International Congress of Applied Mechanics, pp

    Trefftz E., Ein gegenstuck zum ritzschen verfahren. In Proceedings of the 2nd International Congress of Applied Mechanics, pp. 131–137, Zurich, 1926. Orell Fussli Verlag. 39

    work page 1926