Pith. sign in

REVIEW 3 minor 1 cited by

Some new properties of an Active flux type scheme: PamPa

T0 review · 0 major / 3 minor · reviewed 2026-05-18 · grok-4.3

Challenge this review Download PDF Read on arXiv ↗

Pith's one-line read Active flux and PamPa schemes include the discontinuous Galerkin method as a building block while possessing intrinsic bound preservation and summation-by-parts properties.

desk verdict The paper shows active flux schemes embed a discontinuous Galerkin scheme as a building block, gives intrinsic bound preservation for the current pampa variant under standard assumptions, and proves a 1D summation-by-parts property with direct algebra. read the letter →

arxiv 2508.17147 v3 submitted 2025-08-23 math.NA cs.NA

Rémi Abgrall , Philipp Öffner , Yongle Liu This is my paper · ORCID
classification math.NAcs.NA
keywords activefluxpampaschemediscontinuousgalerkinboundpreservingsummationbypartsnumericalschemesconservationlawshigh-ordermethods
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper establishes that Active Flux and PamPa schemes admit a general interpretation in which the discontinuous Galerkin scheme functions as one of their constituent elements. It further shows that the present PamPa variant maintains solution bounds by design, without requiring separate limiters, and verifies this behavior through numerical examples. In one dimension the scheme satisfies the summation-by-parts identity, a property that supports discrete conservation and stability analysis. Readers working on high-order methods for hyperbolic conservation laws would care because these structural links suggest routes to combine existing techniques and to reduce reliance on ad-hoc fixes for maintaining physical solution ranges.

What carries the argument

The general embedding of discontinuous Galerkin schemes inside Active Flux and PamPa formulations, together with the specific PamPa construction that directly enforces bound preservation and the one-dimensional summation-by-parts relation.

What would settle it

A direct calculation of the discrete inner-product matrices for the one-dimensional PamPa scheme that fails to satisfy the summation-by-parts identity, or a computed solution that exits the expected bounds in the absence of any limiter.

Watch

Extended reading notes

Core claim

In full generality the AF/pampa schemes can be interpreted such that the discontinuous Galerkin scheme is one of their building blocks; the current variant of pampa has intrinsic bound preserving properties; and at least in one dimension the pampa scheme has the summation by parts property.

Load-bearing premise

The current specific variant of the pampa scheme admits the stated bound-preserving properties without additional limiting or post-processing steps, and the discontinuous Galerkin interpretation and summation-by-parts property extend directly from the scheme definition.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

0 major / 3 minor

Summary. The manuscript explores new properties of Active Flux (AF) and Point-Average-Moment PolynomiAl-interpreted (PamPa) schemes. It shows that AF/PamPa schemes admit an interpretation in which the discontinuous Galerkin (dG) scheme is one of the building blocks, establishes intrinsic bound-preserving properties for the current PamPa variant (with numerical illustrations), and proves that the PamPa scheme satisfies the summation-by-parts (SBP) property at least in one dimension.

Significance. If the derivations hold, the results supply a structural unification of AF/PamPa schemes with dG methods and confirm stability features (bound preservation via explicit algebraic cancellation under the stated CFL restriction and monotone flux, and SBP via direct summation of discrete operators) without extra limiting steps. These properties are load-bearing for the practical use of such schemes on conservation laws and constitute a clear contribution to the numerical analysis literature.

minor comments (3)
  1. [Introduction] The precise definition of the 'current variant' of PamPa (including the specific reconstruction and flux choices) should be stated explicitly at the beginning of the analysis sections so that the generality claims can be checked without cross-referencing earlier papers.
  2. [Numerical results] In the bound-preservation numerical examples, report the exact CFL numbers employed and include a short table or statement quantifying any observed maximum-principle violations (even if zero) to make the verification fully reproducible.
  3. [SBP property] The 1-D SBP identity is obtained by direct summation; a brief remark on whether the same cancellation pattern extends to the multi-dimensional case or requires additional assumptions would improve the discussion.

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's summary accurately captures our main contributions: the general interpretation of AF/PamPa schemes with the discontinuous Galerkin method as a building block, the intrinsic bound-preserving properties of the current PamPa variant (with numerical illustrations), and the proof of the summation-by-parts property in one dimension. These results provide a structural unification and confirm stability features without additional limiting steps.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; derivations are self-contained from scheme definitions

full rationale

The paper's central claims—the dG embedding via moment matching, intrinsic bound preservation via algebraic cancellation under the stated CFL and monotone flux, and the 1-D SBP identity via direct operator summation—are obtained by explicit algebraic and structural arguments applied to the PamPa scheme definition itself. No load-bearing step reduces to a fitted parameter renamed as prediction, a self-referential definition, or an unverified self-citation chain; each property is shown to follow directly from the discrete operators and update rules without hidden assumptions that presuppose the target result. The derivations remain internally verifiable from the given equations and hold independently of external fitted data.

Assumptions & free parameters 0 free parameters · 0 assumptions · 0 invented entities

Based solely on the abstract, the work derives properties from the existing definition of Active Flux and PamPa schemes without introducing new free parameters, ad-hoc axioms, or invented entities.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Some new properties of an Active flux type scheme: PamPa." pith.science (2026). https://pith.science/paper/2508.17147

@misc{pith2026250817147,
  author       = {Pith},
  title        = {Pith review of: Some new properties of an Active flux type scheme: PamPa},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2508.17147}},
  note         = {Machine review of arXiv:2508.17147}
}
read the original abstract

In this paper, we provide a few new properties of Active Flux (AF)/Point-Average-Moment PolynomiAl-interpreted (\pampa) schemes. First, we show, in full generality, that the AF/pampa schemes can be interpreted in such a way that the discontinuous Galerkin (dG) scheme is one of their building blocks. Secondly we provide intrinsic bound preserving properties of the current variant of pampa. This is also illustrated numerically. Last, we show, at least in one dimension, that the pampa scheme has the summation by part (SBP) property.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Construction of entropy satisfying Active Flux-type methods

    math.NA 2026-07 conditional novelty 6.0 of 10

    Entropy stability of Active Flux-type schemes reduces to an inequality on cell averages alone; a Tadmor-style blend of high- and low-order fluxes enforces it while retaining bound preservation and oscillation control.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages · cited by 1 Pith paper

  1. [1]

    R. Abgrall , A combination of residual distribution and the active flux formulations or a new class of schemes that can combine several writings of the same hyperbolic problem: application to the 1d Euler equations , Commun. Appl. Math. Comput., 5 (2023), pp. 370– 402, https://doi.org/10.1007/s42967-021-00175-w

    work page doi:10.1007/s42967-021-00175-w 2023

  2. [2]

    R. Abgrall , A combination of residual distribution and the Active Flux formulations or a new class of schemes that can combine several writings of the same hyperbolic problem: application to the 1d Euler equations , Commun. Appl. Math. Comput., 5 (2023), pp. 370– 402

    work page 2023

  3. [3]

    Abgrall and W

    R. Abgrall and W. Barsukow , Extensions of active flux to arbitrary order of accuracy , ESAIM, Math. Model. Numer. Anal., 57 (2023), pp. 991–1027, https://doi.org/10.1051/ m2an/2023004

    work page arXiv 2023

  4. [4]

    Abgrall, M

    R. Abgrall, M. Jiao, Y. Liu, and K. Wu , Bound preserving Point-Average-Moment PolynomiAl-interpreted (PAMPA) scheme: one-dimensional case , Commun. Comput. Phys., (2025). In press

    work page 2025

  5. [5]

    Abgrall, J

    R. Abgrall, J. Lin, and Y. Liu , Active flux for triangular meshes for compressible flows problems, Beijing Journal of Pure and Applied Mathematics, (2025). in press, also Arxiv preprint 2312.11271

    work page arXiv 2025

  6. [6]

    Abgrall, J

    R. Abgrall, J. Lin, and Y. Liu , Active flux for triangular meshes for compressible flows problems, Beijing J. of Pure and Appl. Math., 2 (2025), pp. 1–33

    work page 2025

  7. [7]

    Abgrall and Y

    R. Abgrall and Y. Liu , A new approach for designing well-balanced schemes for the shallow water equations: a combination of conservative and primitive formulations , SIAM J. Sci. Comput., 46 (2024), pp. A3375–A3400

    work page 2024

  8. [8]

    Abgrall, Y

    R. Abgrall, Y. Liu, and W. Boscheri , Bound preserving Point-Average-Moment PolynomiAl-interpreted (PAMPA) on polygonal meshes, 2024, https://arxiv.org/abs/2502. 10069, https://arxiv.org/abs/2502.10069

    work page arXiv 2024

Show all 41 references
  1. [9]

    Abgrall, J

    R. Abgrall, J. Nordstr ¨om, P. ¨Offner, and S. Tokareva, Analysis of the SBP-SAT stabi- lization for finite element methods. I: Linear problems , J. Sci. Comput., 85 (2020), p. 28, https://doi.org/10.1007/s10915-020-01349-z. Id/No 43

    2020 doi

  2. [10]

    Ahmed, A

    B. Ahmed, A. Alsaedi, F. Brezzi, L. Marini, and A. Russo , Projectors for Virtual Element Methods, Comput. Math. Appl., 66 (2013)

    2013

  3. [11]

    Barsukow, The active flux scheme for nonlinear problems , J

    W. Barsukow, The active flux scheme for nonlinear problems , J. Sci. Comput., 86 (2021), pp. Paper No. 3, 34, https://doi.org/10.1007/s10915-020-01381-z, https://doi.org/10. 1007/s10915-020-01381-z

    2021 doi

  4. [12]

    Barsukow, P

    W. Barsukow, P. Chandrashekar, C. Klingenberg, and L. Lechner, A generalized active This manuscript is for review purposes only. 25 flux method of arbitrarily high order in two dimensions , 2025, https://arxiv.org/abs/2502. 05101, https://arxiv.org/abs/2502.05101

    2025

  5. [13]

    Barsukow, J

    W. Barsukow, J. Hohm, C. Klingenberg, and P. L. Roe , The active flux scheme on Cartesian grids and its low Mach number limit , J. Sci. Comput., 81 (2019), pp. 594–622, https://doi.org/10.1007/s10915-019-01031-z

    2019 doi

  6. [14]

    Barsukow, C

    W. Barsukow, C. Klingenberg, L. Lechner, J. Nordstr¨om, S. Ortleb, and H. Ranocha, Stability of the active flux method in the framework of summation-by-parts operators, 2025, https://arxiv.org/abs/2507.11068, https://arxiv.org/abs/2507.11068

    2025

  7. [15]

    Beir ˜ao da Veiga, F

    L. Beir ˜ao da Veiga, F. Brezzi, L. D. Marini, and A. Russo , The Hitchhiker’s guide to the virtual element method , Math. Models Methods Appl. Sci., 24 (2014), pp. 1541–1573, https://doi.org/10.1142/S021820251440003X

    2014 doi

  8. [16]

    Calhoun, E

    D. Calhoun, E. Chudzik, and C. Helzel , The Cartesian grid Active Flux method with adaptive mesh refinement, J. Sci. Comput., 94 (2023), p. 54

    2023

  9. [17]

    D. C. Del Rey Fern´andez, P. D. Boom, and D. W. Zingg, A generalized framework for nodal first derivative summation-by-parts operators, J. Comput. Phys., 266 (2014), pp. 214–239, https://doi.org/10.1016/j.jcp.2014.01.038

    2014 doi

  10. [18]

    D. C. Del Rey Fern ´andez, J. E. Hicken, and D. W. Zingg , Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations, Comput. Fluids, 95 (2014), pp. 171–196, https://doi.org/10.1016/j. compfluid.20...

    2014 doi

  11. [19]

    Eyman, Active flux, PhD thesis, University of Michigan, 2013

    T. Eyman, Active flux, PhD thesis, University of Michigan, 2013

    2013

  12. [20]

    Eyman and P

    T. Eyman and P. Roe , Active flux. 49th AIAA Aerospace Science Meeting, 2011

    2011

  13. [21]

    Eyman and P

    T. Eyman and P. Roe , Active flux for systems . 20 th AIAA Computationa Fluid Dynamics Conference, 2011

    2011

  14. [22]

    T. A. Eymann and P. L. Roe , Multidimensional Active Flux schemes , in 21st AIAA Compu- tational Fluid Dynamics Conference, A. I. of Aeronautics and Astronautics, eds., 2013

    2013

  15. [23]

    Glaubitz, J

    J. Glaubitz, J. Nordstr¨om, and P. ¨Offner, Summation-by-parts operators for general func- tion spaces , SIAM J. Numer. Anal., 61 (2023), pp. 733–754, https://doi.org/10.1137/ 22M1470141

    2023

  16. [24]

    Glaubitz, H

    J. Glaubitz, H. Ranocha, A. R. Winters, M. Schlottke-Lakemper, P. ¨Offner, and G. Gassner, Generalized upwind summation-by-parts operators and their application to nodal discontinuous Galerkin methods , J. Comput. Phys., 529 (2025), p. 26, https://doi. org/10.1016/j.jcp.2025.1...

    2025 doi

  17. [25]

    F. He , Towards a New-generation Numerical Scheme for the Com- press- ible Navier-Stokes Equations with the Active Flux Method , PhD thesis, Applied and Interdisciplinary Mathematics, University of Michigan, 2021. https://deepblue.lib.umich.edu/handle/2027.42/169687

    2021

  18. [26]

    Helzel, D

    C. Helzel, D. Kerkmann, and L. Scandurra , A new ADER method inspired by the active flux method, Journal of Scientific Computing, 80 (2019), pp. 35–61

    2019

  19. [27]

    Liu and P

    X.-D. Liu and P. D. Lax , Solution of two-dimensional riemann problems of gas dynamics by positive schemes, SIAM J. Sci. Comput., 19 (1998), pp. 319–340

    1998

  20. [28]

    J. Maeng , On the Advective Component of Active Flux Schemes for Nonlinear Hyperbolic Conservation Laws, PhD thesis, Applied and Interdisciplinary Mathematics, University of Michigan, 2017. https://deepblue.lib.umich.edu/handle/2027.42/138695

    2017

  21. [29]

    Martinez, F

    A. Martinez, F. Piazzon, A. Sommariva, and M. Vianello , Quadrature-based polyno- mial optimization , Optim. Lett., 14 (2020), pp. 1027–1036, https://doi.org/10.1007/ s11590-019-01416-x, hdl.handle.net/11577/3292611

    2020

  22. [30]

    Mattsson, Diagonal-norm upwind SBP operators , J

    K. Mattsson, Diagonal-norm upwind SBP operators , J. Comput. Phys., 335 (2017), pp. 283– 310, https://doi.org/10.1016/j.jcp.2017.01.042

    2017 doi

  23. [31]

    Olsson , Summation by parts, projections, and stability

    P. Olsson , Summation by parts, projections, and stability. I , Math. Comput., 64 (1995), pp. 1035–1065, s23–s26, https://doi.org/10.2307/2153482

    1995 doi

  24. [32]

    Olsson , Summation by parts, projections, and stability

    P. Olsson , Summation by parts, projections, and stability. II , Math. Comput., 64 (1995), pp. 1473–1493, https://doi.org/10.2307/2153366

    1995 doi

  25. [33]

    Sommariva and M

    A. Sommariva and M. Vianello , Gauss-green cubature and moment computation over arbi- trary geometries, J. Comput. Appl. Math., 231 (2009), pp. 886–896, https://doi.org/10. 1016/j.cam.2009.05.014

    2009

  26. [34]

    Sommariva and M

    A. Sommariva and M. Vianello , Cheap and stable quadrature on polyhedral elements . Pre- print, arXiv:2502.03446 [math.NA] (2025), 2025, https://arxiv.org/abs/2502.03446

    2025

  27. [35]

    Sv¨ard and J

    M. Sv¨ard and J. Nordstr ¨om, Review of summation-by-parts schemes for initial-boundary- value problems , J. Comput. Phys., 268 (2014), pp. 17–38, https://doi.org/10.1016/j.jcp. 2014.02.031. Appendix A. VEM approximation: basic facts. This manuscript is for review purposes onl...

    2014 doi

  28. [36]

    wh is a polynomial of degree k on each edge e of P , that is pwhq|e P Pkpeq,

  29. [37]

    wh is continuous on BP ,

  30. [38]

    ş P whmα dx “ ş P π∇whmα dx for |α| “ k ´ 1, k

    ∆ wh P Pk´2pP q, 4. ş P whmα dx “ ş P π∇whmα dx for |α| “ k ´ 1, k. The degrees of freedom are the same as in VkpP q:

  31. [39]

    The value of wh on the vertices of P ,

  32. [40]

    On each edge of P , the value of vk at the k ´ 1 internal points of the k ` 1 Gauss–Lobatto points on this edge,

  33. [41]

    The L2 projection of wh is computable

    The moments up to order k ´ 2 of wh in P , mαpwhq :“ 1 |P | ż P whmα dx, |α| ď k ´ 2. The L2 projection of wh is computable. For any β, if π0pwhq “ ř α,|α|ďk sαmα, we have xπ0pwhq, mβy “ ÿ α,|α|ďk sαxmβ, mαy “ x wh, mβy. The left hand side is computable since the inner product...

Pith tools

Reviewed May 18, 2026 · model on record in the stance chip above.