pith. sign in

arxiv: 2604.28103 · v1 · submitted 2026-04-30 · 🧮 math.NA · cs.NA

Bounded, Commuting, Discrete-trace Preserving Projections

Alexandre Ern , Johnny Guzm\'an , Pratyush Potu This is my paper

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

classification 🧮 math.NA cs.NA
keywords de Rham complexcommuting projectionsdiscrete trace preservationfinite element methodsbounded operatorsgraph norm stabilitystable liftingsextension operators
0
0 comments X p. Extension

The pith

Bounded commuting projections for the 3D de Rham complex preserve piecewise polynomial traces.

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

The paper constructs projection operators for finite element spaces associated with the three-dimensional de Rham complex. These operators commute with the exterior derivative, remain bounded, and reproduce the trace exactly when that trace belongs to the piecewise polynomial trace space. The operators are defined element-wise on the mesh and satisfy a graph-norm stability estimate whose derivative term is controlled solely by oscillations inside a thin layer of elements next to the boundary. Such operators matter because they directly support the construction of stable liftings for polynomial data and the proof that discrete extensions of that data achieve the same order as continuous extensions.

Core claim

We construct bounded, commuting projections for the three-dimensional de Rham complex with the additional property that the projections preserve the trace of functions or fields whenever the trace is a piecewise polynomial in the appropriate trace space. The projections are locally defined and stable in the graph norm, with the part involving the exterior derivative controlled by the oscillation of that derivative inside a narrow strip of elements touching the boundary and weighted by local mesh size. Moreover, the projections are L2-stable locally whenever the exterior derivative itself is piecewise polynomial in the matching space. Two direct applications are the construction of stableLift

What carries the argument

Locally defined bounded commuting projections that preserve discrete traces, ensuring exact boundary reproduction for polynomial inputs while commuting with the exterior derivative.

Load-bearing premise

The mesh admits a well-defined narrow strip of boundary-touching elements and the finite element spaces are chosen so that their trace spaces consist of piecewise polynomials of matching degree.

What would settle it

On a single-tetrahedron mesh, take a piecewise polynomial field whose trace is polynomial of the correct degree and check whether the constructed projection reproduces that trace exactly on the boundary; any failure would show the preservation property does not hold.

read the original abstract

We construct bounded, commuting projections for the three-dimensional de Rham complex with the additional property that the projections preserve the trace of functions/fields if the latter is a piecewise polynomial in the appropriate trace space. The projections are locally defined and stable in the graph norm. More precisely, the part of the graph norm involving the exterior derivative only involves the oscillation of this derivative in a narrow strip of elements touching the boundary and weighted by the local mesh size. Moreover, the projections are $L^2$-stable locally when acting on functions/fields whose exterior derivative is a piecewise polynomial in the appropriate space. We present two salient applications of the present bounded, commuting, discrete-trace preserving projections: the construction of stable liftings of piecewise polynomial data and an optimality result on the discrete versus continuous extension of piecewise polynomial data.

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

Summary. The manuscript constructs bounded, commuting projections for the three-dimensional de Rham complex that additionally preserve discrete traces on the boundary whenever the input is piecewise polynomial in the appropriate trace space. The projections are locally defined and stable in the graph norm, with the exterior-derivative contribution to the norm controlled solely by the oscillation of that derivative within a narrow boundary strip of elements (scaled by local mesh size). They are also locally L²-stable when the exterior derivative is piecewise polynomial. Two applications are presented: stable liftings of piecewise polynomial data and an optimality result comparing discrete versus continuous extensions of such data.

Significance. If the construction and estimates hold, the result supplies a practical tool for finite-element analysis of de Rham complexes, especially for boundary-value problems and extension operators. The localized graph-norm stability (confining the exterior-derivative oscillation to an O(h)-width strip) and the discrete-trace preservation are genuine strengths that simplify proofs of stability for liftings and extensions without global mesh dependencies. The explicit construction via interior commuting projections plus boundary-layer averaging, together with the stated hypotheses (Lipschitz polyhedron, compatible polynomial degrees, admissible boundary strip), makes the work directly usable in applications such as Maxwell or mixed Stokes discretizations.

minor comments (4)
  1. [Abstract] Abstract, sentence on graph-norm stability: the phrasing 'the part of the graph norm involving the exterior derivative only involves the oscillation...' is slightly awkward and could be clarified (e.g., 'is controlled by the oscillation...').
  2. [Construction section] Section describing the averaging operator (likely §3): while the reproduction of polynomials on trace spaces is asserted, an explicit statement of the precise polynomial degree reproduced by the averaging kernel relative to the finite-element degree would make the commutativity verification fully transparent.
  3. [Applications] Applications section: the optimality result on discrete versus continuous extensions is stated abstractly; a short remark or numerical illustration of the constant independence would strengthen the practical impact.
  4. [Preliminaries] Preliminaries: ensure the notation for the de Rham sequence (spaces, exterior derivative, trace operators) is introduced with a single consistent reference before the construction begins.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and for recommending minor revision. The referee accurately summarizes the key features of our construction—bounded commuting projections for the 3D de Rham complex that preserve discrete polynomial traces, with local definition and graph-norm stability in which the exterior-derivative contribution is controlled solely by oscillation within an O(h)-width boundary strip. We are pleased that the localized stability, discrete-trace preservation, and the two applications (stable liftings and optimality of discrete versus continuous extensions) are recognized as strengths that simplify stability arguments in applications such as Maxwell or mixed Stokes discretizations. Since the report lists no specific major comments, we see no need for changes to the manuscript at present. We remain ready to address any minor clarifications the editor may request.

Circularity Check

0 steps flagged

Direct construction of operators; no circularity

full rationale

The manuscript constructs the projections explicitly: a standard interior commuting projection is recalled, then modified only in a boundary strip by a trace-preserving averaging operator whose support is confined to O(h) elements. Commutativity on the trace spaces follows because those spaces are piecewise polynomial of matching degree and the averaging reproduces polynomials exactly; boundedness and localized graph-norm stability follow from inverse inequalities and the strip width. All steps are self-contained algebraic verifications under explicitly listed hypotheses (Lipschitz polyhedron, compatible polynomial degrees, mesh admitting a boundary strip). No parameter fitting, self-referential definitions, or load-bearing self-citations appear in the derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper rests on standard finite-element assumptions for the de Rham complex on shape-regular meshes; no new free parameters, ad-hoc constants, or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption The computational domain is a polyhedral Lipschitz domain admitting a shape-regular simplicial triangulation.
    Required for local definition of the boundary strip and for standard trace theorems in the de Rham sequence.
  • domain assumption The discrete spaces are chosen so that the trace spaces consist of piecewise polynomials of fixed degree.
    Needed for the exact trace-preservation property to hold.

pith-pipeline@v0.9.0 · 5438 in / 1387 out tokens · 94639 ms · 2026-05-07T07:17:55.809414+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

21 extracted references · 11 canonical work pages

  1. [1]

    Discrete extension operators for mixed finite element spaces on locally refined meshes.Math

    Ainsworth, M., Guzm´ an, J., and Sayas, F.-J. Discrete extension operators for mixed finite element spaces on locally refined meshes.Math. Comp.85(2016), 2639–2650.https://doi-org.revproxy.brown.edu/10.1090/ mcom/3074

    2016

  2. [2]

    LocalL2-bounded commuting projections in FEEC.ESAIM Math

    Arnold, D., and Guzm´ an, J. LocalL2-bounded commuting projections in FEEC.ESAIM Math. Model. Numer. Anal.55(2021), 2169–2184.https://doi.org/10.1051/m2an/2021054

    work page doi:10.1051/m2an/2021054 2021

  3. [3]

    Boffi, F

    Boffi, D., Brezzi, F., and Fortin, M.Mixed finite element methods and applications, vol.44ofSpringer Series in Computational Mathematics. Springer, Heidelberg, 2013.https://doi.org/10.1007/978-3-642-36519-5

    work page doi:10.1007/978-3-642-36519-5 2013

  4. [4]

    On traces forH(curl,Ω) in Lipschitz domains.J

    Buffa, A., Costabel, M., and Sheen, D. On traces forH(curl,Ω) in Lipschitz domains.J. Math. Anal. Appl. 276(2002), 845–867.https://doi-org.revproxy.brown.edu/10.1016/S0022-247X(02)00455-9

    work page doi:10.1016/s0022-247x(02)00455-9 2002

  5. [5]

    A stable local commuting projector and optimalhpapproximation estimates inH(curl).Numer

    Chaumont-Frelet, T., and Vohral´ ık, M. A stable local commuting projector and optimalhpapproximation estimates inH(curl).Numer. Math.156(2024), 2293–2342.https://doi.org/10.1007/s00211-024-01431-w

    work page doi:10.1007/s00211-024-01431-w 2024

  6. [6]

    H., and Winther, R

    Christiansen, S. H., and Winther, R. Smoothed projections in finite element exterior calculus.Math. Comp.77 (2008), 813–829.http://dx.doi.org/10.1090/S0025-5718-07-02081-9

    work page doi:10.1090/s0025-5718-07-02081-9 2008

  7. [7]

    Equivalence of local- and global-best approximations, a simple stable local commuting projector, and optimalhpapproximation estimates inH(div).IMA J

    Ern, A., Gudi, T., Smears, I., and Vohral´ ık, M. Equivalence of local- and global-best approximations, a simple stable local commuting projector, and optimalhpapproximation estimates inH(div).IMA J. Numer. Anal. 42(2022), 1023–1049.http://dx.doi.org/10.1093/imanum/draa103

    work page doi:10.1093/imanum/draa103 2022

  8. [8]

    Mollification in strongly Lipschitz domains with application to continuous and discrete de Rham complexes.Comput

    Ern, A., and Guermond, J.-L. Mollification in strongly Lipschitz domains with application to continuous and discrete de Rham complexes.Comput. Methods Appl. Math.16(2016), 51–75

    2016

  9. [9]

    Approximation and Interpolation, vol.72ofTexts in Applied Mathematics

    Ern, A., and Guermond, J.-L.Finite Elements I. Approximation and Interpolation, vol.72ofTexts in Applied Mathematics. Springer International Publishing, Springer Nature Switzerland AG, 2021.https://doi-org/10. 1007/978-3-030-56341-7

    2021

  10. [10]

    Discrete Poincar´ e inequalities: a review on proofs, equivalent formulations, and behavior of constants.IMA J

    Ern, A., Guzm´ an, J., Potu, P., and Vohral´ ık, M. Discrete Poincar´ e inequalities: a review on proofs, equivalent formulations, and behavior of constants.IMA J. Numer. Anal.(2026). to appear,https://hal.inria.fr/ hal-04837821

    2026

  11. [11]

    LocalL2-bounded commuting projections using discrete local problems on Alfeld splits.Found

    Ern, A., Guzm´ an, J., Potu, P., and Vohral´ ık, M. LocalL2-bounded commuting projections using discrete local problems on Alfeld splits.Found. Comput. Math.(2026). to appear,https://hal.inria.fr/hal-04931497

    2026

  12. [12]

    S., and Winther, R

    Falk, R. S., and Winther, R. Local bounded cochain projections.Math. Comp.83(2014), 2631–2656.http: //dx.doi.org/10.1090/S0025-5718-2014-02827-5

    work page doi:10.1090/s0025-5718-2014-02827-5 2014

  13. [13]

    J., and Licht, M

    Gawlik, E., Holst, M. J., and Licht, M. W. Local finite element approximation of Sobolev differential forms. ESAIM Math. Model. Numer. Anal.55(2021), 2075–2099.https://doi-org.revproxy.brown.edu/10.1051/ m2an/2021034

    work page arXiv 2021

  14. [14]

    Discrete regular decompositions of tetrahedral discrete 1-forms

    Hiptmair, R., and Pechstein, C. Discrete regular decompositions of tetrahedral discrete 1-forms. InMaxwell’s equations—analysis and numerics, vol.24ofRadon Ser. Comput. Appl. Math.De Gruyter, Berlin, [2019] ©2019, pp. 199–258

    2019

  15. [15]

    Oxford Univ

    Monk, P.Finite element methods for Maxwell’s equations. Numerical Mathematics and Scientific Compu- tation. Oxford University Press, New York, 2003.https://doi-org.revproxy.brown.edu/10.1093/acprof: oso/9780198508885.001.0001

    work page doi:10.1093/acprof: 2003

  16. [16]

    Mixed finite elements inR3.Numer

    N´ ed´ elec, J.-C. Mixed finite elements inR3.Numer. Math.35(1980), 315–341

    1980

  17. [17]

    A mixed finite element method for 2nd order elliptic problems

    Raviart, P.-A., and Thomas, J.-M. A mixed finite element method for 2nd order elliptic problems. InMathe- matical aspects of finite element methods (Proc. Conf., Consiglio Naz. delle Ricerche (C.N.R.), Rome, 1975). Springer, Berlin, 1977, pp. 292–315. Lecture Notes in Math., Vol. 606. DISCRETE-TRACE PRESERVING COMMUTING PROJECTIONS 31

    1975

  18. [18]

    A multilevel decomposition result inH(curl)

    Sch¨ oberl, J. A multilevel decomposition result inH(curl). InMultigrid, multilevel and multiscale methods, Proceedings of the8 th European Multigrid Conference(2005), P. Wesseling, C. Oosterlee, and P. Hemker, Eds., TU Delft

    2005

  19. [19]

    47 [Tem79] R

    Scott, L. R., and Zhang, S. Finite element interpolation of nonsmooth functions satisfying boundary conditions. Math. Comp.54(1990), 483–493.https://doi-org.revproxy.brown.edu/10.2307/2008497

    work page doi:10.2307/2008497 1990

  20. [20]

    Princeton University Press, Princeton, NJ, 1957

    Whitney, H.Geometric integration theory. Princeton University Press, Princeton, NJ, 1957

    1957

  21. [21]

    J., and Farin, G

    Worsey, A. J., and Farin, G. Ann-dimensional Clough-Tocher interpolant.Constr. Approx.3(1987), 99–110. https://doi-org.revproxy.brown.edu/10.1007/BF01890556. CERMICS, CNRS, ENPC, Institut Polytechnique de Paris, 77455 Marne-la-Vall ´ee, France & Inria Paris, 48 rue Barrault, 75647 Paris, France Email address:alexandre.ern@enpc.fr Division of Applied Mathe...

    work page doi:10.1007/bf01890556 1987