pith. machine review for the scientific record. sign in

arxiv: 2512.16352 · v2 · submitted 2025-12-18 · 🧮 math.NA · cs.NA

Recognition: no theorem link

Conserving mass, momentum, and energy for the Benjamin-Bona-Mahony, Korteweg-de Vries, and nonlinear Schr\"odinger equations

Hendrik Ranocha , David I. Ketcheson
Authors on Pith no claims yet

Pith reviewed 2026-05-16 21:35 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords invariant-preserving discretizationFourier Galerkin methodorthogonal projectionrelaxation methodsBenjamin-Bona-Mahony equationKorteweg-de Vries equationnonlinear Schrödinger equationstructure-preserving integrators
0
0 comments X

The pith

High-order Fourier schemes using projection and relaxation conserve mass, momentum, and energy for the BBM, KdV, and NLS equations.

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

The paper develops arbitrarily high-order numerical methods for three nonlinear dispersive wave equations that keep the physical invariants of mass, momentum, and energy exactly conserved to machine precision. Spatial discretization uses Fourier Galerkin methods, while time integration combines orthogonal projection with a relaxation step that corrects the solution after each step. Because the resulting schemes stay essentially explicit and retain their formal accuracy, they support long-time integrations without the artificial growth of errors that occurs when invariants are allowed to drift.

Core claim

The central claim is that Fourier Galerkin semi-discretizations combined with orthogonal projection and relaxation produce arbitrarily high-order, essentially explicit time integrators that conserve mass, momentum, and energy for the Benjamin-Bona-Mahony, Korteweg-de Vries, and nonlinear Schrödinger equations, as well as a hyperbolic approximation of the latter. The conservation holds up to numerical round-off, is proven for the chosen discretizations, and is verified in numerical tests that also show markedly slower growth of solution errors over long integration intervals.

What carries the argument

Orthogonal projection onto the manifold defined by the invariants followed by a relaxation correction, applied to Fourier-Galerkin spatial semi-discretizations

If this is right

  • The schemes conserve mass, momentum, and energy to machine precision for the BBM, KdV, and NLS equations.
  • Conservation is achieved while the methods remain arbitrarily high-order and essentially explicit.
  • Long-term simulations exhibit slower growth of global errors than non-conserving integrators.
  • The same construction works for a hyperbolic approximation of the nonlinear Schrödinger equation.

Where Pith is reading between the lines

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

  • The technique may extend directly to other Hamiltonian PDEs whose invariants can be written as quadratic or cubic functionals.
  • Because the correction is inexpensive, the methods could be paired with adaptive time-stepping without destroying conservation.
  • Application to problems with non-periodic boundaries would require replacing the Fourier basis while retaining the projection step.

Load-bearing premise

That the projection-relaxation correction preserves the three invariants exactly without lowering the formal order or stability of the underlying time integrator.

What would settle it

A long-time numerical integration in which any of the three invariants drifts beyond round-off error while the scheme is run at its design order.

Figures

Figures reproduced from arXiv: 2512.16352 by David I. Ketcheson, Hendrik Ranocha.

Figure 8
Figure 8. Figure 8: Change of invariants for moving one- and two-gray-soliton solutions of the NLS equation with mass- and energy-conserving relaxation. The time integration is performed with the fifth-order method of [52] with Δ𝑡 = 0.04. one-soliton solution (although this is not guaranteed by relaxation). However, the momentum varies clearly for the two solitons if relaxation is used to conserve only the mass and energy, se… view at source ↗
read the original abstract

We propose and study a class of arbitrarily high-order numerical discretizations that preserve multiple invariants and are essentially explicit (they do not require the solution of any large systems of algebraic equations). In space, we use Fourier Galerkin methods, while in time we use a combination of orthogonal projection and relaxation. We prove and numerically demonstrate the conservation properties of the method by applying it to the Benjamin-Bona-Mahony, Korteweg-de Vries, and nonlinear Schr\"odinger (NLS) PDEs as well as a hyperbolic approximation of NLS. For each of these equations, the proposed schemes conserve mass, momentum, and energy up to numerical precision. We show that this conservation leads to reduced growth of numerical errors for long-term simulations.

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

Summary. The manuscript proposes arbitrarily high-order numerical schemes for the Benjamin-Bona-Mahony, Korteweg-de Vries, nonlinear Schrödinger, and hyperbolic NLS equations. These combine Fourier-Galerkin spatial discretization with a time-stepping approach based on orthogonal projection followed by relaxation. The central claims are that the resulting methods conserve mass, momentum, and energy exactly (up to floating-point precision), that the proofs of conservation hold for each equation, and that the invariant preservation reduces long-term error growth while remaining essentially explicit.

Significance. If the conservation proofs are complete and the schemes retain their formal order and stability, the work supplies a practical family of structure-preserving integrators for dispersive PDEs. The explicit character and the numerical evidence across four model problems are strengths that could influence long-time simulation practice in nonlinear wave equations.

major comments (2)
  1. [§3.2] §3.2 (proof of momentum conservation for KdV): the argument that the relaxation step leaves the discrete momentum invariant unchanged relies on the specific inner-product structure after projection; a short explicit calculation showing that the relaxation parameter does not alter the L2 inner product of the projected solution with itself would strengthen the claim.
  2. [§4.3] §4.3 (numerical order verification): the reported convergence rates for the NLS test are consistent with the design order, but the tables do not separate the contribution of the projection-relaxation step from the underlying integrator; an additional column or remark quantifying any order reduction would confirm that the conservation mechanism does not degrade accuracy.
minor comments (3)
  1. The notation for the relaxation parameter is introduced without a dedicated symbol list; adding a short table of symbols would improve readability.
  2. Figure 2 caption refers to 'energy drift' but the y-axis label is 'relative energy error'; harmonizing the terminology would avoid minor confusion.
  3. The statement that the method is 'arbitrarily high-order' is repeated in the abstract and introduction; a single forward reference to the theorem establishing the order would suffice.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and the recommendation for minor revision. The two major comments are constructive and we address them point by point below, indicating the changes we will make to the manuscript.

read point-by-point responses

  1. Referee: [§3.2] §3.2 (proof of momentum conservation for KdV): the argument that the relaxation step leaves the discrete momentum invariant unchanged relies on the specific inner-product structure after projection; a short explicit calculation showing that the relaxation parameter does not alter the L2 inner product of the projected solution with itself would strengthen the claim.

    Authors: We agree that an explicit calculation will strengthen the presentation. In the revised manuscript we will insert a short calculation immediately after the statement of momentum conservation in §3.2. Let u^{n+1/2} denote the solution after the orthogonal projection step and let λ be the relaxation parameter. We will show directly that the momentum inner product satisfies ⟨u^{n+1},u^{n+1}⟩ = ⟨u^{n+1/2},u^{n+1/2}⟩ because the relaxation update is constructed to be orthogonal to the momentum functional gradient; the explicit algebra is ⟨u^{n+1},u^{n+1}⟩ = ⟨u^{n+1/2} + λ v, u^{n+1/2} + λ v⟩ = ⟨u^{n+1/2},u^{n+1/2}⟩ + 2λ⟨u^{n+1/2},v⟩ + λ²⟨v,v⟩ and the choice of λ forces the cross term to vanish while preserving the quadratic invariant. This addition will make the invariance transparent without altering the existing proof structure. revision: yes

  2. Referee: [§4.3] §4.3 (numerical order verification): the reported convergence rates for the NLS test are consistent with the design order, but the tables do not separate the contribution of the projection-relaxation step from the underlying integrator; an additional column or remark quantifying any order reduction would confirm that the conservation mechanism does not degrade accuracy.

    Authors: We accept the suggestion. In the revised §4.3 we will add a brief remark stating that the projection-relaxation step is a consistent perturbation of order equal to the underlying time integrator and therefore does not reduce the formal order. To make this quantitative we will augment the NLS convergence table with an extra column that reports the observed temporal orders obtained when the relaxation step is omitted (i.e., using only the base integrator followed by projection). The new column will confirm that the orders remain identical to those already reported, thereby demonstrating that the conservation mechanism introduces no order reduction. revision: yes

Circularity Check

0 steps flagged

No significant circularity: conservation proved from explicit method construction

full rationale

The paper constructs Fourier-Galerkin spatial discretizations combined with orthogonal projection and relaxation in time specifically to enforce exact conservation of mass, momentum, and energy for the target PDEs. It then proves these properties hold (up to roundoff) directly from the definitions of the projection and relaxation steps, without any fitted parameters, self-referential definitions, or load-bearing self-citations that reduce the central claim to prior unverified assertions. The numerical demonstrations are consistent with the analytic proofs rather than serving as the sole justification. This is a standard structure-preserving design whose invariants follow by construction from the chosen operators, not from any circular reduction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The method rests on standard Fourier analysis and the existence of an orthogonal projection onto the invariant manifold; no new free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (2)
  • standard math Fourier Galerkin truncation converges for the target PDEs under sufficient smoothness
    Invoked implicitly when claiming high-order accuracy and exact conservation in the discrete space.
  • domain assumption The relaxation step can be performed without destroying the formal order of the time integrator
    Required for the claim that the schemes remain high-order while enforcing conservation.

pith-pipeline@v0.9.0 · 5438 in / 1263 out tokens · 25587 ms · 2026-05-16T21:35:51.497687+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

97 extracted references · 97 canonical work pages · 2 internal anchors

  1. [1]

    A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes

    R. Abgrall. “A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes.” In:Journal of Computational Physics372 (2018), pp. 640–666.doi:10.1016/j.jcp.2018.06.031. arXiv: 1711.10358 [math.NA]

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1016/j.jcp.2018.06.031 2018

  2. [2]

    Relaxation Deferred Correction Methods and their Applications to Residual Distribution Schemes

    R. Abgrall, E. L. Mélédo, P. Öffner, and D. Torlo. “Relaxation Deferred Correction Methods and their Applications to Residual Distribution Schemes.” In:The SMAI Journal of Compu- tational Mathematics8 (2022), pp. 125–160.doi:10.5802/smai- jcm.82. arXiv:2106.05005 [math.NA]

    work page doi:10.5802/smai- 2022

  3. [3]

    Reinterpretation and Extension of Entropy Cor- rection Terms for Residual Distribution and Discontinuous Galerkin Schemes: Application to Structure Preserving Discretization

    R. Abgrall, P. Öffner, and H. Ranocha. “Reinterpretation and Extension of Entropy Cor- rection Terms for Residual Distribution and Discontinuous Galerkin Schemes: Application to Structure Preserving Discretization.” In:Journal of Computational Physics453 (Mar. 2022), p. 110955.doi:10.1016/j.jcp.2022.110955. arXiv:1908.04556 [math.NA]. 20

    work page doi:10.1016/j.jcp.2022.110955 2022

  4. [4]

    M. J. Ablowitz, B. Prinari, and A. D. Trubatch.Discrete and continuous nonlinear Schrödinger systems. Vol. 302. Cambridge University Press, 2004.doi:10.1017/CBO9780511546709

    work page doi:10.1017/cbo9780511546709 2004

  5. [5]

    First-order exact solutions of the nonlinear Schrödinger equation in the normal-dispersion regime

    N. Akhmediev and A. Ankiewicz. “First-order exact solutions of the nonlinear Schrödinger equation in the normal-dispersion regime.” In:Physical Review A47.4 (1993), p. 3213.doi: 10.1103/PhysRevA.47.3213

    work page doi:10.1103/physreva.47.3213 1993

  6. [6]

    High-order mass-, energy-and momentum- conserving methods for the nonlinear Schrödinger equation

    G. Akrivis, B. Li, R. Tang, and H. Zhang. “High-order mass-, energy-and momentum- conserving methods for the nonlinear Schrödinger equation.” In:Journal of Computational Physics532 (2025), p. 113974.doi:10.1016/j.jcp.2025.113974

    work page doi:10.1016/j.jcp.2025.113974 2025

  7. [7]

    On fully discrete Galerkin methods of second-order temporal accuracy for the nonlinear Schrödinger equation

    G. D. Akrivis, V. A. Dougalis, and O. A. Karakashian. “On fully discrete Galerkin methods of second-order temporal accuracy for the nonlinear Schrödinger equation.” In:Numerische Mathematik59.1 (1991), pp. 31–53.doi:10.1007/BF01385769

    work page doi:10.1007/bf01385769 1991

  8. [8]

    Error propagation when approximating multi-solitons: The KdV equation as a case study

    J Álvarez and A. Durán. “Error propagation when approximating multi-solitons: The KdV equation as a case study.” In:Applied Mathematics and Computation217.4 (2010), pp. 1522– 1539.doi:10.1016/j.amc.2009.06.033

    work page doi:10.1016/j.amc.2009.06.033 2010

  9. [9]

    B.D.AndrewsandP.E.Farrell.Conservativeanddissipativediscretisationsofmulti-conservative ODEs and GENERIC systems. Nov. 2025. arXiv:2511.23266 [math.NA]

    work page arXiv 2025

  10. [10]

    Error propagation in the numerical integration of solitary waves. Theregularizedlongwaveequation

    A Araújo and A. Durán. “Error propagation in the numerical integration of solitary waves. Theregularizedlongwaveequation.”In:AppliedNumericalMathematics36.2-3(2001),pp.197– 217.doi:10.1016/S0168-9274(99)00148-8

    work page doi:10.1016/s0168-9274(99)00148-8 2001

  11. [11]

    Implicit-explicit Runge-Kuttamethods fortime- dependent partial differential equations

    U. M.Ascher, S.J. Ruuth,and R.J. Spiteri.“Implicit-explicit Runge-Kuttamethods fortime- dependent partial differential equations.” In:Applied Numerical Mathematics25.2-3 (1997), pp. 151–167.doi:10.1016/S0168-9274(97)00056-1

    work page doi:10.1016/s0168-9274(97)00056-1 1997

  12. [12]

    High-ordermass-andenergy-conservingmethodsforthenonlinear Schrödinger equation

    G.Bai,J.Hu,andB.Li.“High-ordermass-andenergy-conservingmethodsforthenonlinear Schrödinger equation.” In:SIAM Journal on Scientific Computing46.2 (2024), A1026–A1046. doi:10.1137/22M152178X

    work page doi:10.1137/22m152178x 2024

  13. [13]

    Model equations for long waves in nonlinear dispersive systems

    T. B. Benjamin, J. L. Bona, and J. J. Mahony. “Model equations for long waves in nonlinear dispersive systems.” In:Philosophical Transactions of the Royal Society of London. Series A, MathematicalandPhysicalSciences272.1220(1972),pp.47–78.doi:10.1098/rsta.1972.0032

    work page doi:10.1098/rsta.1972.0032 1972

  14. [14]

    ArelaxationschemeforthenonlinearSchrödingerequation

    C.Besse.“ArelaxationschemeforthenonlinearSchrödingerequation.”In:SIAMJournalon Numerical Analysis42.3 (2004), pp. 934–952.doi:10.1137/S0036142901396521

    work page doi:10.1137/s0036142901396521 2004

  15. [15]

    Energy-preserving methods for nonlinearSchrödingerequations

    C. Besse, S. Descombes, G. Dujardin, and I. Lacroix-Violet. “Energy-preserving methods for nonlinearSchrödingerequations.”In:IMAJournalofNumericalAnalysis41.1(2021),pp.618– 653.doi:10.1093/imanum/drz067

    work page doi:10.1093/imanum/drz067 2021

  16. [16]

    Julia:AFreshApproachtoNumerical Computing

    J.Bezanson,A.Edelman,S.Karpinski,andV.B.Shah.“Julia:AFreshApproachtoNumerical Computing.” In:SIAM Review59.1 (2017), pp. 65–98.doi:10.1137/141000671. arXiv:1411. 1607 [cs.MS]

    work page doi:10.1137/141000671 2017

  17. [17]

    A.Biswas,L.S.Busaleh,D.I.Ketcheson,C.Muñoz-Moncayo,andM.Rajvanshi.AHyperbolic Approximation of the Nonlinear Schrödinger Equation. 2025. arXiv:2505.21424 [math.AP]

    work page arXiv 2025

  18. [18]

    AccurateSolutionoftheNonlinearSchrödingerEquationvia Conservative Multiple-Relaxation ImEx Methods

    A.BiswasandD.I.Ketcheson.“AccurateSolutionoftheNonlinearSchrödingerEquationvia Conservative Multiple-Relaxation ImEx Methods.” In:SIAM Journal of Scientific Computing 46.6 (2024), A3827–A3848.doi:10.1137/23M1598118. arXiv:2309.02324 [math.NA]

    work page doi:10.1137/23m1598118 2024

  19. [20]

    Traveling-wave solutions and structure-preservingnumericalmethodsforahyperbolicapproximationoftheKorteweg-de Vriesequation

    A. Biswas, D. I. Ketcheson, H. Ranocha, and J. Schütz. “Traveling-wave solutions and structure-preservingnumericalmethodsforahyperbolicapproximationoftheKorteweg-de Vriesequation.”In:JournalofScientificComputing103(May2025),p.90.doi:10.1007/s10915- 025-02898-x. arXiv:2412.17117 [math.NA]. 21

    work page doi:10.1007/s10915-

  20. [21]

    Bleecke, A

    S. Bleecke, A. Biswas, D. I. Ketcheson, J. Schütz, and H. Ranocha.Asymptotic-preserving and energy-conserving methods for a hyperbolic approximation of the BBM equation. Nov. 2025. arXiv: 2511.10044 [math.NA]

    work page arXiv 2025

  21. [22]

    Projection methods preserving Lya- punov functions

    M. Calvo, M. Laburta, J. I. Montijano, and L. Rández. “Projection methods preserving Lya- punov functions.” In:BIT Numerical Mathematics50.2 (2010), pp. 223–241.doi:10 . 1007 / s10543-010-0259-3

    work page 2010

  22. [23]

    Conserved quantities of some Hamiltonian wave equations after full discretiza- tion

    B. Cano. “Conserved quantities of some Hamiltonian wave equations after full discretiza- tion.” In:NumerischeMathematik103.2 (2006), pp.197–223.doi:10.1007/s00211-006-0680- 3

    work page doi:10.1007/s00211-006-0680- 2006

  23. [24]

    Energy-preserving Runge-Kutta methods

    E. Celledoni, R. I. McLachlan, D. I. McLaren, B. Owren, G. R. W. Quispel, and W. Wright. “Energy-preserving Runge-Kutta methods.” In:ESAIM: Mathematical Modelling and Numeri- cal Analysis (M2AN)43 (2009), pp. 645–649.doi:10.1051/m2an/2009020

    work page doi:10.1051/m2an/2009020 2009

  24. [25]

    A new conservative discontinuous Galerkin method via implicit penalization for the generalized Korteweg-de Vries equation

    Y. Chen, B. Dong, and R. Pereira. “A new conservative discontinuous Galerkin method via implicit penalization for the generalized Korteweg-de Vries equation.” In:SIAM Journal on Numerical Analysis60.6 (2022), pp. 3078–3098.doi:10.1137/22M1470827

    work page doi:10.1137/22m1470827 2022

  25. [26]

    Mass-andenergy-preservingexponentialRunge-Kutta methods for the nonlinear Schrödinger equation

    J.Cui,Z.Xu,Y.Wang,andC.Jiang.“Mass-andenergy-preservingexponentialRunge-Kutta methods for the nonlinear Schrödinger equation.” In:Applied Mathematics Letters112 (2021), p. 106770.doi:10.1016/j.aml.2020.106770

    work page doi:10.1016/j.aml.2020.106770 2021

  26. [27]

    Makie.jl: Flexible high-performance data visualization for Julia

    S. Danisch and J. Krumbiegel. “Makie.jl: Flexible high-performance data visualization for Julia.” In:Journal of Open Source Software6.65 (2021), p. 3349.doi:10.21105/joss.03349

    work page doi:10.21105/joss.03349 2021

  27. [28]

    Accuracy and conservation properties in numerical integration: the case of the Korteweg-de Vries equation

    J De Frutos and J. M. Sanz-Serna. “Accuracy and conservation properties in numerical integration: the case of the Korteweg-de Vries equation.” In:Numerische Mathematik75.4 (1997), pp. 421–445.doi:10.1007/s002110050247

    work page doi:10.1007/s002110050247 1997

  28. [29]

    Dekker and J

    K. Dekker and J. G. Verwer.Stability of Runge-Kutta methods for stiff nonlinear differential equations. Vol. 2. CWI Monographs. Amsterdam: North-Holland, 1984

    work page 1984

  29. [31]

    The(𝑛+1)/2Rule for Dealiasing in the Split-Step Fourier Methods for 𝑛-Wave Interactions

    S. Derevyanko. “The(𝑛+1)/2Rule for Dealiasing in the Split-Step Fourier Methods for 𝑛-Wave Interactions.” In:IEEE Photonics Technology Letters20.23 (2008), pp. 1929–1931.doi: 10.1109/LPT.2008.2005420

    work page doi:10.1109/lpt.2008.2005420 2008

  30. [32]

    ExtendedLagrangianapproachforthedefocusing nonlinearSchrödingerequation

    F.Dhaouadi,N.Favrie,andS.Gavrilyuk.“ExtendedLagrangianapproachforthedefocusing nonlinearSchrödingerequation.”In:StudiesinAppliedMathematics142.3(2019),pp.336–358. doi:10.1111/sapm.12238

    work page doi:10.1111/sapm.12238 2019

  31. [33]

    Doehring, H

    D. Doehring, H. Ranocha, and M. Torrilhon.Paired Explicit Relaxation Runge-Kutta Methods: Entropy-Conservative and Entropy-Stable High-Order Optimized Multirate Time Integration. July

    work page

  32. [34]

    arXiv:2507.04991 [math.NA]

    work page arXiv

  33. [35]

    Ahigh-orderfullydiscreteschemefortheKorteweg-deVries equationwithatime-steppingprocedureofRunge-Kutta-compositiontype

    V.A.DougalisandÁ.Durán.“Ahigh-orderfullydiscreteschemefortheKorteweg-deVries equationwithatime-steppingprocedureofRunge-Kutta-compositiontype.”In:IMAJournal of Numerical Analysis42.4 (2022), pp. 3022–3057.doi:10.1093/imanum/drab060

    work page doi:10.1093/imanum/drab060 2022

  34. [36]

    The numerical integration of relative equilibrium solu- tions.ThenonlinearSchrödingerequation

    A. Durán and J. M. Sanz-Serna. “The numerical integration of relative equilibrium solu- tions.ThenonlinearSchrödingerequation.”In:IMAJournalofNumericalAnalysis20.2(2000), pp. 235–261.doi:10.1093/imanum/20.2.235

    work page doi:10.1093/imanum/20.2.235 2000

  35. [37]

    Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations

    D. C. D. R. Fernández, 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.” In:Computers & Fluids95 (2014), pp. 171–196.doi:10 . 1016 / j . compfluid.2014.02.016. 22

    work page 2014

  36. [38]

    Numerical preservation of multiple local conservation laws

    G. Frasca-Caccia and P. E. Hydon. “Numerical preservation of multiple local conservation laws.” In:Applied Mathematics and Computation403 (2021), p. 126203.doi:10.1016/j.amc. 2021.126203

    work page doi:10.1016/j.amc 2021

  37. [39]

    Simple bespoke preservation of two conservation laws

    G. Frasca-Caccia and P. E. Hydon. “Simple bespoke preservation of two conservation laws.” In:IMAJournalofNumericalAnalysis40.2(2020),pp.1294–1329.doi:10.1093/imanum/dry087

    work page doi:10.1093/imanum/dry087 2020

  38. [40]

    The design and implementation of FFTW3

    M. Frigo and S. G. Johnson. “The design and implementation of FFTW3.” In:Proceedings of the IEEE93.2 (2005), pp. 216–231.doi:10.1109/JPROC.2004.840301

    work page doi:10.1109/jproc.2004.840301 2005

  39. [41]

    CRC Press, 2010.doi:10.1201/b10387

    D.FurihataandT.Matsuo.Discretevariationalderivativemethod:astructure-preservingnumerical method for partial differential equations. CRC Press, 2010.doi:10.1201/b10387

    work page doi:10.1201/b10387 2010

  40. [42]

    Hyperbolic approximation of the BBM equation

    S. Gavrilyuk and K.-M. Shyue. “Hyperbolic approximation of the BBM equation.” In:Non- linearity35.3 (2022), p. 1447.doi:10.1088/1361-6544/ac4c49

    work page doi:10.1088/1361-6544/ac4c49 2022

  41. [43]

    Giesselmann and H

    J. Giesselmann and H. Ranocha.Convergence of hyperbolic approximations to higher-order PDEs for smooth solutions. Aug. 2025. arXiv:2508.04112 [math.NA]

    work page arXiv 2025

  42. [44]

    Bespoke finite difference schemes that preserve multiple conservation laws

    T. J. Grant. “Bespoke finite difference schemes that preserve multiple conservation laws.” In:LMS Journal of Computation and Mathematics18.1 (2015), pp. 372–403.doi:10 . 1112 / S1461157015000078

    work page 2015

  43. [45]

    Energy-preservingvariantofcollocationmethods

    E.Hairer.“Energy-preservingvariantofcollocationmethods.”In:JournalofNumericalAnal- ysis, Industrial and Applied Mathematics5 (2010), pp. 73–84

    work page 2010

  44. [46]

    Hairer, C

    E. Hairer, C. Lubich, and G. Wanner.Geometric Numerical Integration: Structure-Preserving AlgorithmsforOrdinaryDifferentialEquations.Vol.31.SpringerSeriesinComputationalMath- ematics. Berlin Heidelberg: Springer-Verlag, 2006.doi:10.1007/3-540-30666-8

    work page doi:10.1007/3-540-30666-8 2006

  45. [47]

    Crank-NicolsonGalerkinapproximationstononlinearSchrödinger equationswithroughpotentials

    P.HenningandD.Peterseim.“Crank-NicolsonGalerkinapproximationstononlinearSchrödinger equationswithroughpotentials.”In:MathematicalModelsandMethodsinAppliedSciences27.11 (2017), pp. 2147–2184.doi:10.1142/S0218202517500415

    work page doi:10.1142/s0218202517500415 2017

  46. [48]

    Introduction to the Hirota bilinear method

    J. Hietarinta. “Introduction to the Hirota bilinear method.” In:Integrability of Nonlinear Sys- tems:ProceedingsoftheCIMPASchoolPondicherryUniversity,India,8–26January1996.Springer, 2007, pp. 95–103.doi:10.1007/BFb0113694

    work page doi:10.1007/bfb0113694 2007

  47. [49]

    Exact solution of the Korteweg—de Vries equation for multiple collisions of solitons

    R. Hirota. “Exact solution of the Korteweg—de Vries equation for multiple collisions of solitons.” In:Physical Review Letters27.18 (1971), p. 1192.doi:10.1103/PhysRevLett.27. 1192

    work page doi:10.1103/physrevlett.27 1971

  48. [50]

    S. G. Johnson.Notes on FFT-based differentiation.https://math.mit.edu/~stevenj/fft- deriv.pdf. 2011

    work page 2011

  49. [51]

    Pseudo-spectral methods and linear instabilities in reaction- diffusionfronts

    W. B. Jones and J. J. O’Brien. “Pseudo-spectral methods and linear instabilities in reaction- diffusionfronts.”In:Chaos:AnInterdisciplinaryJournalofNonlinearScience6.2(1996),pp.219– 228.doi:10.1063/1.166167

    work page doi:10.1063/1.166167 1996

  50. [52]

    Entropy-PreservingandEntropy-StableRelaxationIMEX and Multirate Time-Stepping Methods

    S.KangandE.M.Constantinescu.“Entropy-PreservingandEntropy-StableRelaxationIMEX and Multirate Time-Stepping Methods.” In:Journal of Scientific Computing93 (2022), p. 23. doi:10.1007/s10915-022-01982-w. arXiv:2108.08908 [math.NA]

    work page doi:10.1007/s10915-022-01982-w 2022

  51. [53]

    Higher-order additive Runge-Kutta schemes for ordinary differential equations

    C. A. Kennedy and M. H. Carpenter. “Higher-order additive Runge-Kutta schemes for ordinary differential equations.” In:Applied Numerical Mathematics136 (2019), pp. 183–205. doi:10.1016/j.apnum.2018.10.007

    work page doi:10.1016/j.apnum.2018.10.007 2019

  52. [54]

    Relaxation Runge-Kutta Methods: Conservation and stability for Inner-Product Norms

    D. I. Ketcheson. “Relaxation Runge-Kutta Methods: Conservation and Stability for Inner- Product Norms.” In:SIAM Journal on Numerical Analysis57.6 (2019), pp. 2850–2870.doi: 10.1137/19M1263662. arXiv:1905.09847 [math.NA]

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1137/19m1263662 2019

  53. [55]

    On Using Quasi-Newton Algorithms of the Broyden Class for Model-to-test Correlation

    J. Klement. “On Using Quasi-Newton Algorithms of the Broyden Class for Model-to-test Correlation.” In:Journal of Aerospace Technology and Management6.4 (2014), pp. 407–414.doi: 10.5028/jatm.v6i4.373. 23

    work page doi:10.5028/jatm.v6i4.373 2014

  54. [56]

    Nonlinear and linear conservative finite difference schemes for regularized long wave equation

    S. Koide and D. Furihata. “Nonlinear and linear conservative finite difference schemes for regularized long wave equation.” In:Japan Journal of Industrial and Applied Mathematics26.1 (2009), p. 15.doi:10.1007/BF03167544

    work page doi:10.1007/bf03167544 2009

  55. [57]

    D. A. Kopriva.Implementing Spectral Methods for Partial Differential Equations: Algorithms for Scientists and Engineers. New York: Springer Science & Business Media, 2009.doi:10.1007/ 978-90-481-2261-5

    work page 2009

  56. [58]

    On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves

    D. J. Korteweg and G. De Vries. “On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves.” In:The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science39.240 (1895), pp. 422–443.doi:10. 1080/14786449508620739

    work page

  57. [59]

    Structure-Preserving Numerical Methods for Two Nonlinear Systems of Dispersive Wave Equations

    J. Lampert and H. Ranocha. “Structure-Preserving Numerical Methods for Two Nonlinear Systems of Dispersive Wave Equations.” In:Computational Science and Engineering2 (Nov. 2025), p. 2.doi:10.1007/s44207-025-00006-3. arXiv:2402.16669 [math.NA]

    work page doi:10.1007/s44207-025-00006-3 2025

  58. [60]

    Relaxation Exponential Rosenbrock-Type Methods for Oscillatory Hamil- tonian Systems

    D. Li and X. Li. “Relaxation Exponential Rosenbrock-Type Methods for Oscillatory Hamil- tonian Systems.” In:SIAM Journal on Scientific Computing45.6 (2023), A2886–A2911.doi: 10.1137/22M1511345

    work page doi:10.1137/22m1511345 2023

  59. [61]

    Implicit-explicit relaxation Runge-Kutta methods: construction, analysis and applications to PDEs

    D. Li, X. Li, and Z. Zhang. “Implicit-explicit relaxation Runge-Kutta methods: construction, analysis and applications to PDEs.” In:Mathematics of Computation(2022).doi:10 . 1090 / mcom/3766

    work page 2022

  60. [62]

    H. Li, X. Li, K. Schratz, and B. Wang.Time-Relaxation Structure-Preserving Explicit Low- RegularityIntegratorsfortheNonlinearSchrödingerEquation.2025.arXiv:2510.02963 [math.NA]

    work page arXiv 2025

  61. [63]

    Linearhigh-orderenergy-preservingschemesforthenonlinear Schrödinger equation with wave operator using the scalar auxiliary variable approach

    X.Li,Y.Gong,andL.Zhang.“Linearhigh-orderenergy-preservingschemesforthenonlinear Schrödinger equation with wave operator using the scalar auxiliary variable approach.” In: Journal of Scientific Computing88.1 (2021), p. 20.doi:10.1007/s10915-021-01533-9

    work page doi:10.1007/s10915-021-01533-9 2021

  62. [64]

    ResolvingEntropyGrowthfromIterativeMethods

    V.Linders,H.Ranocha,andP.Birken.“ResolvingEntropyGrowthfromIterativeMethods.” In:BIT Numerical Mathematics63 (4 Sept. 2023).doi:10.1007/s10543-023-00992-w. arXiv: 2302.13579 [math.NA]

    work page doi:10.1007/s10543-023-00992-w 2023

  63. [65]

    Error analysis for spectral approximation of the Korteweg- de Vries equation

    Y Maday and A Quarteroni. “Error analysis for spectral approximation of the Korteweg- de Vries equation.” In:ESAIM: Mathematical Modelling and Numerical Analysis22.3 (1988), pp. 499–529.doi:10.1051/m2an/1988220304991

    work page doi:10.1051/m2an/1988220304991 1988

  64. [66]

    Geometric integration using discrete gradi- ents

    R. I. McLachlan, G. Quispel, and N. Robidoux. “Geometric integration using discrete gradi- ents.”In:PhilosophicalTransactionsoftheRoyalSocietyofLondon.SeriesA:Mathematical,Physical and Engineering Sciences357.1754 (1999), pp. 1021–1045.doi:10.1098/rsta.1999.0363

    work page doi:10.1098/rsta.1999.0363 1999

  65. [67]

    Aconservativefully-discretenumer- icalmethodfortheregularizedshallowwaterwaveequations

    D.Mitsotakis,H.Ranocha,D.I.Ketcheson,andE.Süli.“Aconservativefully-discretenumer- icalmethodfortheregularizedshallowwaterwaveequations.”In:SIAMJournalonScientific Computing42 (2 Apr. 2021), B508–B537.doi:10 . 1137 / 20M1364606. arXiv:2009 . 09641 [math.NA]

    work page 2021

  66. [68]

    Korteweg-deVriesequationandgeneraliza- tions. II. Existence of conservation laws and constants of motion

    R.M.Miura,C.S.Gardner,andM.D.Kruskal.“Korteweg-deVriesequationandgeneraliza- tions. II. Existence of conservation laws and constants of motion.” In:Journal of Mathematical physics9.8 (1968), pp. 1204–1209.doi:10.1063/1.1664701

    work page doi:10.1063/1.1664701 1968

  67. [69]

    Instead of rewriting foreign code for machine learning, auto- matically synthesize fast gradients

    W. Moses and V. Churavy. “Instead of rewriting foreign code for machine learning, auto- matically synthesize fast gradients.” In:Advances in Neural Information Processing Systems33 (2020), pp. 12472–12485

    work page 2020

  68. [70]

    Reverse-mode automatic differentiation and optimization of GPU kernels via Enzyme

    W. S. Moses, V. Churavy, L. Paehler, J. Hückelheim, S. H. K. Narayanan, M. Schanen, and J. Doerfert. “Reverse-mode automatic differentiation and optimization of GPU kernels via Enzyme.” In:Proceedings of the International Conference for High Performance Computing, Net- working, Storage and Analysis. 2021, pp. 1–16.doi:10.1145/3458817.3476165. 24

    work page doi:10.1145/3458817.3476165 2021

  69. [71]

    Euler operators and conservation laws of the BBM equation

    P. J. Olver. “Euler operators and conservation laws of the BBM equation.” In:Mathematical ProceedingsoftheCambridgePhilosophicalSociety.Vol.85.1.CambridgeUniversityPress.1979, pp. 143–160.doi:10.1017/S0305004100055572

    work page doi:10.1017/s0305004100055572 1979

  70. [72]

    On the elimination of aliasing in finite-difference schemes by filtering high- wavenumbercomponents

    S. A. Orszag. “On the elimination of aliasing in finite-difference schemes by filtering high- wavenumbercomponents.”In:JournalofAtmosphericSciences28.6(1971),pp.1074–1074.doi: 10.1175/1520-0469(1971)028<1074:OTEOAI>2.0.CO;2

    work page doi:10.1175/1520-0469(1971)028 1971

  71. [73]

    NonlinearSolve.jl: High-Performance and Robust Solvers for Systems of Non- linear Equations in Julia

    A. Pal, F. Holtorf, A. Larsson, T. Loman, Utkarsh, F. Schäfer, Q. Qu, A. Edelman, and C. Rackauckas. “NonlinearSolve.jl: High-Performance and Robust Solvers for Systems of Non- linear Equations in Julia.” In: (2024).doi:10.48550/arXiv.2403.16341. arXiv:2403.16341 [math.NA]

    work page doi:10.48550/arxiv.2403.16341 2024

  72. [74]

    A new class of energy-preserving numerical integration methods

    G. Quispel and D. I. McLaren. “A new class of energy-preserving numerical integration methods.” In:Journal of Physics A: Mathematical and Theoretical41.4 (2008), p. 045206.doi: 10.1088/1751-8113/41/4/045206

    work page doi:10.1088/1751-8113/41/4/045206 2008

  73. [75]

    SummationByPartsOperators.jl: A Julia library of provably stable semidis- cretization techniques with mimetic properties

    H. Ranocha. “SummationByPartsOperators.jl: A Julia library of provably stable semidis- cretization techniques with mimetic properties.” In:Journal of Open Source Software6.64 (Aug. 2021), p. 3454.doi:10 . 21105 / joss . 03454.url:https : / / github . com / ranocha / SummationByPartsOperators.jl

    work page 2021

  74. [76]

    Fully-Discrete Explicit Locally Entropy-Stable SchemesfortheCompressibleEulerandNavier-StokesEquations

    H. Ranocha, L. Dalcin, and M. Parsani. “Fully-Discrete Explicit Locally Entropy-Stable SchemesfortheCompressibleEulerandNavier-StokesEquations.”In:ComputersandMathe- maticswithApplications80.5(July2020),pp.1343–1359.doi:10.1016/j.camwa.2020.06.016. arXiv:2003.08831 [math.NA]

    work page doi:10.1016/j.camwa.2020.06.016 2020

  75. [77]

    Energy Stability of Explicit Runge-Kutta Methods for Nonautonomous or Nonlinear Problems

    H. Ranocha and D. I. Ketcheson. “Energy Stability of Explicit Runge-Kutta Methods for Nonautonomous or Nonlinear Problems.” In:SIAM Journal on Numerical Analysis58.6 (Nov. 2020), pp. 3382–3405.doi:10.1137/19M1290346. arXiv:1909.13215 [math.NA]

    work page doi:10.1137/19m1290346 2020

  76. [78]

    Ranocha and D

    H. Ranocha and D. I. Ketcheson.High-order mass- and energy-conserving methods for the non- linear Schrödinger equation and its hyperbolization. Oct. 2025. arXiv:2510.14335 [math.NA]

    work page arXiv 2025

  77. [79]

    Relaxation Runge-Kutta Methods for Hamiltonian Prob- lems

    H. Ranocha and D. I. Ketcheson. “Relaxation Runge-Kutta Methods for Hamiltonian Prob- lems.”In:JournalofScientificComputing84.1(July2020).doi:10.1007/s10915-020-01277-y. arXiv:2001.04826 [math.NA]

    work page doi:10.1007/s10915-020-01277-y 2001

  78. [80]

    Conservingmass,momentum,and energy for the Benjamin-Bona-Mahony, Korteweg-de Vries, and nonlinear Schrödinger equations

    H.RanochaandD.I.Ketcheson.Reproducibilityrepositoryfor"Conservingmass,momentum,and energy for the Benjamin-Bona-Mahony, Korteweg-de Vries, and nonlinear Schrödinger equations". https://github.com/ranocha/2025_BBM_KdV_NLS. 2025.doi:10.5281/zenodo.17936837

    work page doi:10.5281/zenodo.17936837 2025

  79. [81]

    General Relaxation Methods for Initial-Value ProblemswithApplicationtoMultistepSchemes

    H. Ranocha, L. Lóczi, and D. I. Ketcheson. “General Relaxation Methods for Initial-Value ProblemswithApplicationtoMultistepSchemes.”In:NumerischeMathematik146(Oct.2020), pp. 875–906.doi:10.1007/s00211-020-01158-4. arXiv:2003.03012 [math.NA]

    work page doi:10.1007/s00211-020-01158-4 2020

  80. [82]

    A Broad Class of Conservative Numerical Methods for Dispersive Wave Equations

    H. Ranocha, D. Mitsotakis, and D. I. Ketcheson. “A Broad Class of Conservative Numerical Methods for Dispersive Wave Equations.” In:Communications in Computational Physics29.4 (Feb.2021),pp.979–1029.doi:10.4208/cicp.OA-2020-0119.arXiv:2006.14802 [math.NA]

    work page doi:10.4208/cicp.oa-2020-0119.arxiv:2006.14802 2021

Showing first 80 references.