pith. sign in

arxiv: 1906.10622 · v1 · pith:VR5ZIBJJnew · submitted 2019-06-25 · 🌊 nlin.PS

The ampsys tool of pde2path

Hannes Uecker , Daniel Wetzel This is my paper

Pith reviewed 2026-05-25 15:47 UTC · model grok-4.3

classification 🌊 nlin.PS
keywords amplitude equationsTuring bifurcationSwift-Hohenberg equationreaction-diffusion systemspattern formationMATLAB toolnumerical continuation
0
0 comments X

The pith

The ampsys tool automates computation of amplitude system coefficients for Turing bifurcations in Swift-Hohenberg type equations and reaction-diffusion systems.

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

The paper presents the ampsys tool that automates the calculation of coefficients in amplitude systems near Turing bifurcations. This covers scalar equations of Swift-Hohenberg type and generalizations as well as reaction-diffusion systems with any number of components. The tool needs only the PDE specification as input and supports symbolic computations in applicable cases. Examples demonstrate its use for 1D, 2D and 3D wave vector lattices. A reader would care because manual coefficient derivation grows tedious quickly in higher dimensions and for multi-component systems, so automation removes that barrier to studying the resulting amplitude equations.

Core claim

The ampsys tool automates the computation of coefficients of amplitude systems for Turing bifurcations for scalar equations of Swift-Hohenberg type and generalizations, and reaction-diffusion systems with an arbitrary number of components. The tool is designed to require minimal user input, and for a number of cases can also deal with symbolic computations. After a brief review of the setup of amplitude systems the tool is explained by a number of 1D, 2D and 3D examples over various wave vector lattices.

What carries the argument

The ampsys tool, which encodes multiple-scales or center-manifold derivations to extract amplitude equation coefficients from a user-supplied PDE with minimal further input.

If this is right

  • Amplitude equations for 2D and 3D lattices can be obtained directly from the PDE without manual algebra.
  • Systems with an arbitrary number of reaction-diffusion components become accessible for amplitude analysis.
  • Symbolic mode reduces transcription errors when the problem permits closed-form coefficients.
  • The resulting amplitude systems can be fed into continuation software for further bifurcation study.

Where Pith is reading between the lines

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

  • The same automation approach could be adapted to other local bifurcation types such as Hopf or steady-state bifurcations with different normal forms.
  • Extending the tool to handle spatially varying coefficients or additional nonlinear terms would cover a wider range of physical models.
  • Batch processing of families of related PDEs could reveal how coefficients vary with parameters without repeated manual setup.

Load-bearing premise

The multiple-scales and center-manifold derivations have been correctly implemented in the code for the two problem classes covered.

What would settle it

Apply the tool to a standard 1D Swift-Hohenberg equation whose amplitude coefficients are already known from hand calculation and check whether the outputs match.

Figures

Figures reproduced from arXiv: 1906.10622 by Daniel Wetzel, Hannes Uecker.

Figure 1
Figure 1. Figure 1: First three layers (wave vectors generated by terms up to cubic order) for (a) Square lattice; (b) hexagonal lattice; (c) 8-fold quasilattice. The thick black stars are the basic wave vectors (on the critical circle), the blue (red) dots are generated by quadratic (cubic) interactions. Acknowledgment. The work of DW was supported by the DFG under Grant No. 264671738. 2 Some amplitude systems on simple latt… view at source ↗
Figure 2
Figure 2. Figure 2: Critical wave vectors (black stars), and first three ’layers’ of the generated quasi-lattice. 3.3 A damped Kuramoto-Sivashinsky type of equation, demo KS The demo KS deals with damped/driven Kuramoto-Sivashinsky (KS) [KY76, Siv88] type of equations, e.g., in 1D, ∂tu = −(1 + ∂ 2 x ) 2u + λu + c2∂x(u 2 ), (40) i.e., a SH equation with a convective nonlinearity, which gives another important class of pattern … view at source ↗
read the original abstract

The computation of coefficients of amplitude systems for Turing bifurcations is a straightforward but sometimes elaborate task, in particular for 2D or 3D wave vector lattices. The Matlab tool ampsys automates such computations for two classes of problems, namely scalar equations of Swift-Hohenberg type and generalizations, and reaction-diffusion systems with an arbitrary number of components. The tool is designed to require minimal user input, and for a number of cases can also deal with symbolic computations. After a brief review of the setup of amplitude systems we explain the tool by a number of 1D, 2D and 3D examples over various wave vector lattices.

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 describes the ampsys Matlab tool (part of the pde2path package) for automating the calculation of coefficients in amplitude equations for Turing bifurcations. It covers two problem classes: scalar Swift-Hohenberg-type equations (and generalizations) and reaction-diffusion systems with arbitrary numbers of components. The tool requires minimal user input, supports some symbolic computations, and is illustrated via 1D/2D/3D examples on various wave-vector lattices after a brief review of the multiple-scales setup.

Significance. If the implementation is correct, the tool would be a useful practical aid for researchers computing amplitude systems in pattern-formation problems, especially when algebraic complexity grows with lattice dimension or component count. It could lower the barrier to systematic studies of bifurcations in higher-dimensional or multi-component systems.

major comments (2)
  1. [Examples (throughout)] The central claim that ampsys correctly automates the derivations for arbitrary-component RD systems rests on unverified code paths. The examples section supplies illustrative runs but contains no cross-checks of output coefficients against independent hand derivations, known analytic results from the literature, or test cases with N>3 components where term counts grow combinatorially.
  2. [Tool description and setup review] No error analysis, convergence tests for the linear solves over wave-vector lattices, or discussion of how the symbolic engine scales or fails for the claimed 'arbitrary number of components' is provided. This leaves the weakest assumption (correct automated derivation without user intervention) untested precisely where algebraic load is highest.
minor comments (2)
  1. [Abstract] The abstract and introduction could more explicitly state the precise classes of amplitude equations (e.g., real vs. complex Ginzburg-Landau type) that the tool targets.
  2. [Examples] Figure captions and example listings would benefit from explicit statements of the input PDE, the chosen lattice, and the expected output form of the amplitude system.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the major comments point by point below and will revise the manuscript to incorporate additional verification and discussion as outlined.

read point-by-point responses

  1. Referee: [Examples (throughout)] The central claim that ampsys correctly automates the derivations for arbitrary-component RD systems rests on unverified code paths. The examples section supplies illustrative runs but contains no cross-checks of output coefficients against independent hand derivations, known analytic results from the literature, or test cases with N>3 components where term counts grow combinatorially.

    Authors: We acknowledge that the examples in the current manuscript are primarily for illustration and do not contain explicit cross-checks against hand derivations or literature results, nor tests with N>3. This is a fair observation. We will add a dedicated verification subsection in the revised version, including a direct comparison of tool output against a known analytic result from the literature for a two-component system and an additional example with N=4 components to address the combinatorial growth concern. revision: yes

  2. Referee: [Tool description and setup review] No error analysis, convergence tests for the linear solves over wave-vector lattices, or discussion of how the symbolic engine scales or fails for the claimed 'arbitrary number of components' is provided. This leaves the weakest assumption (correct automated derivation without user intervention) untested precisely where algebraic load is highest.

    Authors: We agree that the manuscript lacks discussion of scaling and potential limitations of the symbolic engine. In revision we will add a short section on computational aspects, reporting observed scaling behavior with component number and lattice size from our tests, and noting practical limits for very large N. Regarding error analysis and convergence tests, the underlying linear systems are solved exactly (symbolically or numerically) within the multiple-scales framework; we will clarify this and include numerical accuracy checks on selected examples. revision: yes

Circularity Check

0 steps flagged

Tool description paper with no derivation chain or fitted predictions

full rationale

The paper presents ampsys as a software tool that automates standard multiple-scales and center-manifold calculations for amplitude equations in known problem classes. It reviews the setup of amplitude systems and demonstrates usage via examples over 1D/2D/3D lattices, but contains no new first-principles derivations, no parameter fitting, and no predictions that reduce to inputs by construction. No self-citation load-bearing steps, uniqueness theorems, or ansatzes are invoked in a manner that creates circularity. The central claim is implementational (the code performs the algebra), which is independent of any self-referential reduction. This is the normal non-circular outcome for a methods/tool paper.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The work describes a software tool rather than a mathematical derivation; no free parameters, axioms, or invented entities are introduced or required by the abstract.

pith-pipeline@v0.9.0 · 5628 in / 985 out tokens · 22326 ms · 2026-05-25T15:47:28.257972+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

24 extracted references · 24 canonical work pages

  1. [1]

    Cross and P.C

    M.C. Cross and P.C. Hohenberg. Pattern formation outside equi\-li\-brium. Rev. Mod. Phys. , 65:854--1190, 1993

    work page 1993

  2. [2]

    T. K. Callahan and E. Knobloch. Symmetry-breaking bifurcations on cubic lattices. Nonlinearity , 10:1179--1216, 1997

    work page 1997

  3. [3]

    T. K. Callahan and E. Knobloch. Pattern formation in three-dimensional reaction-diffusion systems. Phys. D , 132(3):339--362, 1999

    work page 1999

  4. [4]

    de Witt, T

    H. de Witt, T. Dohnal, J.D.M. Rademacher, H. Uecker, and D. Wetzel. pde2path - Quickstart guide and reference card , 2018

    work page 2018

  5. [5]

    Golubitsky and I

    M. Golubitsky and I. Stewart. The symmetry perspective . Birkh\"auser, Basel, 2002

    work page 2002

  6. [6]

    R.B. Hoyle. Pattern formation . Cambridge University Press., 2006

    work page 2006

  7. [7]

    Iooss and A

    G. Iooss and A. M. Rucklidge. On the existence of quasipattern solutions of the S wift- H ohenberg equation. J. Nonlinear Sci. , 20(3):361--394, 2010

    work page 2010

  8. [8]

    Kuramoto and T

    Y. Kuramoto and T. Yamada. Pattern formation in oscillatory chemical reactions. Progress of theoretical physics , 56(3):724--740, sep 1976

    work page 1976

  9. [9]

    L.M. Pismen. Patterns and interfaces in dissipative dynamics . Springer , 2006

    work page 2006

  10. [10]

    Prigogine and R

    I. Prigogine and R. Lefever . Symmetry Breaking Instabilities in Dissipative Systems. II . J. Chem. Phys , 48(4):1695--1700, 1968

    work page 1968

  11. [11]

    Subramanian, A.J

    P. Subramanian, A.J. Archer, E. Knobloch, and A.M. Rucklidge. Three-dimensional icosahedral phase field quasicrystal. Phys. Rev. Lett. , 117:075501, 2016

    work page 2016

  12. [12]

    Swift and P.C

    J. Swift and P.C. Hohenberg. Hydrodynamic fluctuations at the convective instability. Physical Review A , 15(1):319--328, 1977

    work page 1977

  13. [13]

    Sivashinsky

    G. Sivashinsky. Nonlinear analysis of hydrodynamic instability in laminar flames— I . D erivation of basic equations. Acta Astronautica , pages 459 -- 488, 1988

    work page 1988

  14. [14]

    Schneider and H

    G. Schneider and H. Uecker. Nonlinear PDE -- a dynamical systems approach , volume 182 of Graduate Studies Mathematics . AMS, 2017

    work page 2017

  15. [15]

    H. Uecker. User guide on H opf bifurcation and time periodic orbits with pde2path, 2018. Available at p2phome

    work page 2018

  16. [16]

    H. Uecker. Hopf bifurcation and time periodic orbits with pde2path -- algorithms and applications . Comm. in Comp. Phys , 25(3):812--852, 2019

    work page 2019

  17. [17]

    H. Uecker. Pattern formation with pde2path -- a tutorial, 2019

    work page 2019

  18. [18]

    H. Uecker. www.staff.uni-oldenburg.de/hannes.uecker/pde2path , 2019

    work page 2019

  19. [19]

    Uecker and D

    H. Uecker and D. Wetzel. Numerical results for snaking of patterns over patterns in some 2D Selkov-Schnakenberg Reaction-Diffusion systems . SIADS , 13(1):94--128, 2014

    work page 2014

  20. [20]

    Uecker and D

    H. Uecker and D. Wetzel. Snaking branches of planar BCC fronts in the 3D Brusselator . preprint , 2019

    work page 2019

  21. [21]

    Uecker, D

    H. Uecker, D. Wetzel, and J.D.M. Rademacher. pde2path -- a Matlab package for continuation and bifurcation in 2D elliptic systems . NMTMA , 7:58--106, 2014

    work page 2014

  22. [22]

    Verdasca, A

    J. Verdasca, A. de Wit, G. Dewel, and P. Borckmans. Reentrant hexagonal T uring structures. Phys. Lett. A , 168(194):194--198, 1992

    work page 1992

  23. [23]

    D. Wetzel. Tristability between stripes, up-hexagons, and down-hexagons and snaking bifurcation branches of spatial connections between up- and down-hexagons. Phys. Rev. E , 97(062221), 2018

    work page 2018

  24. [24]

    L. Yang, M. Dolnik, A.M. Zhabotinsky, and I.R. Epstein. Pattern formation arising from interactions between T uring and wave instabilities. The Journal of Chemical Physics , 117(15):7259--7265, 2002

    work page 2002