pith. sign in

arxiv: 2501.13680 · v4 · submitted 2025-01-23 · 💻 cs.SC · math.AG· math.CA

Projecting dynamical systems via a support bound

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

classification 💻 cs.SC math.AGmath.CA
keywords Newton polytopedifferential eliminationpolynomial dynamical systemsminimal differential equationevaluation-interpolationsymbolic computationprojection
0
0 comments X

The pith

A bound on the Newton polytope of the minimal differential equation for a coordinate in a polynomial dynamical system depends only on dimension and the input degrees d and D.

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

The paper addresses the problem of projecting a polynomial dynamical system onto one chosen coordinate by computing the minimal differential equation it satisfies. It establishes an explicit bound on the Newton polytope of this minimal equation, expressed in terms of the system dimension together with the degrees d and D of the polynomials governing the distinguished coordinate and the remaining variables. The bound is shown to be tight whenever d is at most D or the system is two-dimensional. The authors then convert the bound into an evaluation-interpolation algorithm that computes the minimal equation and solves instances beyond the reach of existing differential-elimination software.

Core claim

For a polynomial dynamical system, a bound is given for the Newton polytope of the minimal differential equation satisfied by a chosen coordinate. The bound depends on the dimension of the model and the degrees d and D of the polynomials defining the dynamics of the chosen coordinate and the remaining coordinates. The bound is sharp if d ≤ D or the model is planar. This bound is used to design an algorithm for computing the minimal equation via the evaluation-interpolation paradigm.

What carries the argument

The explicit bound on the Newton polytope of the minimal equation, which directly determines the monomials needed in the evaluation-interpolation procedure.

If this is right

  • The minimal equation can be recovered by evaluating the system at sufficiently many points and solving a linear algebra problem whose size is controlled by the bound.
  • Instances that exceed the capacity of current differential-elimination tools become tractable.
  • When d ≤ D or the system is planar the computed equation is guaranteed to have the smallest possible support.
  • The same bound supplies an a-priori estimate of the computational cost of the projection step.

Where Pith is reading between the lines

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

  • The support bound may be reusable as a complexity measure for other coordinate-projection tasks in algebraic differential systems.
  • In modeling applications the bound offers a way to decide in advance whether projecting a given coordinate is feasible before running the algorithm.
  • The evaluation-interpolation strategy could be combined with sparse interpolation techniques to further reduce the number of required sample points.

Load-bearing premise

The input consists of polynomial vector fields whose degrees are exactly d for the distinguished coordinate and D for the others.

What would settle it

A concrete polynomial dynamical system with d > D in dimension greater than two whose minimal equation has a monomial outside the predicted Newton polytope.

read the original abstract

For a polynomial dynamical system, we study the problem of computing the minimal differential equation satisfied by a chosen coordinate (in other words, projecting the system on the coordinate). This problem can be viewed as a special case of the general elimination problem for systems of differential equations and appears in applications to modeling and control. We give a bound for the Newton polytope of such minimal equation. Our bound depends on the dimension of the model and the degrees $d$ and $D$ of the polynomials defining the dynamics of the chosen coordinate and the remaining coordinates, respectively. We show that our bound is sharp if $d \leqslant D$ or the model is planar. We further use this bound to design an algorithm for computing the minimal equation following the evaluation-interpolation paradigm. We demonstrate that our implementation of the algorithm can tackle problems which are out of reach for the state-of-the-art software for differential elimination.

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

Summary. The manuscript studies the projection of a polynomial dynamical system onto one coordinate by computing the minimal annihilating differential equation satisfied by that coordinate. It derives an explicit upper bound on the Newton polytope of this minimal equation that depends only on the ambient dimension n and the input degrees d (for the distinguished coordinate) and D (for the remaining coordinates). The bound is shown to be sharp when d ≤ D or when n=2, and is then used to certify an evaluation-interpolation algorithm whose implementation is reported to solve instances beyond the reach of existing differential-elimination software.

Significance. If the stated bound and its sharpness hold, the work supplies a parameter-free support bound that directly enables a certified, practical algorithm in the evaluation-interpolation paradigm for differential elimination. This is a concrete advance for symbolic computation, with clear relevance to modeling and control applications. The explicit dependence solely on n, d, and D, together with the sharpness results and the reported computational gains, constitute the main strengths.

minor comments (2)
  1. [Abstract] The explicit form of the bound (presumably stated in §3 or §4) is described rather than displayed in the abstract; adding the formula would make the central contribution immediately visible to readers.
  2. [§6] Section 6 benchmark tables would benefit from explicit reporting of hardware specifications, timeout values, and memory limits to allow direct reproduction of the claimed performance advantage.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary, significance assessment, and recommendation to accept the manuscript.

Circularity Check

0 steps flagged

No significant circularity; bound derived from input parameters

full rationale

The central result is an explicit upper bound on the Newton polytope of the minimal annihilating equation, stated to depend only on ambient dimension n and the exact input degrees d (distinguished coordinate) and D (remaining coordinates). The paper proves containment, establishes sharpness when d ≤ D or n=2, and applies the bound to an evaluation-interpolation algorithm whose correctness follows directly from the support containment. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations appear; the derivation is self-contained against the stated polynomial vector-field inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No free parameters, axioms, or invented entities are introduced or required by the abstract; the work operates within standard algebraic-geometry notions of Newton polytopes and differential elimination.

pith-pipeline@v0.9.0 · 5681 in / 1061 out tokens · 38323 ms · 2026-05-23T05:24:18.585741+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

65 extracted references · 65 canonical work pages

  1. [1]

    M. F. Atiyah and I. G. Macdonald. Introduction to commutative algebra . Addi- son–Wesley, 1969

  2. [2]

    B¨ achler, V

    T. B¨ achler, V. Gerdt, M. Lange-Hegermann, and D. Robertz. Thomas decomposition of algebraic and differential systems. In Computer Algebra in Scientific Computing , pages 31–54. 2010. URL https://doi.org/10.1007/978-3-642-15274-0_4

  3. [3]

    B¨ achler, V

    T. B¨ achler, V. Gerdt, M. Lange-Hegermann, and D. Robertz. Algorithmic Thomas decomposition of algebraic and differential systems.Journal of Symbolic Computation, 47(10):1233–1266, 2012. URL https://doi.org/10.1016/j.jsc.2011.12.043

  4. [4]

    Bezanson, A

    J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah. Julia: A fresh approach to numerical computing. SIAM review, 59(1):65–98, 2017. URL https://doi.org/10. 1137/141000671

  5. [5]

    Bostan, F

    A. Bostan, F. Chyzak, F. Ollivier, B. Salvy, E. Schost, and A. Sedoglavic. Fast com- putation of power series solutions of systems of differential equations. In Proceedings of the Eighteenth Annual ACM-SIAM Symposium on Discrete Algorithms, SODA ’07, page 1012–1021, 2007. URL https://dl.acm.org/doi/10.5555/1283383.1283492

  6. [6]

    F. Boulier. BLAD: Biblioth` eques Lilloises d’Alg` ebre Diff´ erentielle. URL https: //pro.univ-lille.fr/francois-boulier/logiciels/blad/

  7. [7]

    F. Boulier. Differential elimination and biological modelling. In Gr¨ obner Bases in Symbolic Analysis, page 109–138. De Gruyter, 2007. URL http://dx.doi.org/10. 1515/9783110922752.109

  8. [8]

    Boulier, D

    F. Boulier, D. Lazard, F. Ollivier, and M. Petitot. Representation for the radical of a finitely generated differential ideal. In Proceedings of the 1995 International Symposium on Symbolic and Algebraic Computation - ISSAC’95 , 1995. URL https: //doi.org/10.1145/220346.220367

  9. [9]

    Boulier, D

    F. Boulier, D. Lazard, F. Ollivier, and M. Petitot. Computing representations for radicals of finitely generated differential ideals. Applicable Algebra in Engineering, Communication and Computing , 20(1):73–121, 2009. URL https://doi.org/10. 1007/s00200-009-0091-7. – 34 – Projecting dynamical systems via a support bound Y. Mukhina, G. Pogudin

  10. [10]

    Bourbaki

    N. Bourbaki. ´El´ ements de math´ ematique: Les structures fondamentales de l’analyse; Alg` ebre; Polynomes et fractions rationnelles; Chapitre 5: Corps commutatifs . Her- mann, 1950

  11. [11]

    Buchberger

    B. Buchberger. A criterion for detecting unnecessary reductions in the construction of gr¨ obner-bases. InInternational Symposium on Symbolic and Algebraic Manipulation , pages 3–21. Springer, 1979

  12. [12]

    S. Cabay. Exact solution of linear equations. In Proceedings of the Second ACM Sym- posium on Symbolic and Algebraic Manipulation , SYMSAC ’71, page 392–398. Asso- ciation for Computing Machinery, 1971. URL https://doi.org/10.1145/800204. 806310

  13. [13]

    Carra-Ferro

    G. Carra-Ferro. A resultant theory for the systems of two ordinary algebraic differen- tial equations. Applicable Algebra in Engineering, Communication and Computing , 8 (6):539–560, 1997. URL https://doi.org/10.1007/s002000050090

  14. [14]

    Carr` a Ferro

    G. Carr` a Ferro. A survey on differential Gr¨ obner bases. InGr¨ obner Bases in Sym- bolic Analysis, page 77–108. De Gruyter, 2007. URL http://dx.doi.org/10.1515/ 9783110922752.77

  15. [15]

    D. Cox, J. Little, D. O’Shea, and M. Sweedler. Ideals, varieties, and algorithms , volume 3. Springer, 1997

  16. [16]

    D. A. Cox, J. Little, and D. O’Shea. Using algebraic geometry, volume 185. Springer Science & Business Media, 2005

  17. [17]

    Dickenstein, M

    A. Dickenstein, M. I. Herrero, and B. Mourrain. Curve valuations and mixed volumes in the implicitization of rational varieties.Journal of Algebra, 612:691–721, 2022. URL https://doi.org/10.1016/j.jalgebra.2022.08.026

  18. [18]

    R. Dong, C. Goodbrake, H. A. Harrington, and G. Pogudin. Differential elimination for dynamical models via projections with applications to structural identifiability. SIAM Journal on Applied Algebra and Geometry , 7(1):194–235, Mar. 2023. URL http://dx.doi.org/10.1137/22M1469067

  19. [19]

    D’Alfonso, G

    L. D’Alfonso, G. Jeronimo, and P. Solern´ o. On the complexity of the resolvent rep- resentation of some prime differential ideals. Journal of Complexity , 22(3):396–430,

  20. [20]

    URL http://dx.doi.org/10.1016/j.jco.2005.10.002

  21. [21]

    Ehrenborg and G.-C

    R. Ehrenborg and G.-C. Rota. Apolarity and canonical forms for homogeneous poly- nomials. European Journal of Combinatorics , 14(3):157–181, 1993. URL https: //doi.org/10.1006/eujc.1993.1022

  22. [22]

    I. Z. Emiris, T. Kalinka, and C. Konaxis. Implicitization of curves and surfaces using predicted support. In Proceedings of the 2011 International Workshop on Symbolic- Numeric Computation, volume 4 of ISSAC ’11, page 137–146. ACM, June 2012. URL http://dx.doi.org/10.1145/2331684.2331705

  23. [23]

    I. Z. Emiris, T. Kalinka, C. Konaxis, and T. Luu Ba. Sparse implicitization by interpo- lation: Characterizing non-exactness and an application to computing discriminants. Computer-Aided Design, 45(2):252–261, Feb. 2013. URL http://dx.doi.org/10. 1016/j.cad.2012.10.008. – 35 – Projecting dynamical systems via a support bound Y. Mukhina, G. Pogudin

  24. [24]

    A. Esterov. Engineered complete intersections: slightly degenerate Bernstein– Kouchnirenko–Khovanskii, 2024. URL https://arxiv.org/abs/2401.12099

  25. [25]

    Esterov and A

    A. Esterov and A. Khovanskii. Elimination theory and Newton polytopes. Func- tional Analysis and Other Mathematics , 2(1):45–71, 2008. URL https://doi.org/ 10.1007/s11853-008-0015-2

  26. [26]

    Fieker, W

    C. Fieker, W. Hart, T. Hofmann, and F. Johansson. Nemo/Hecke: Computer algebra and number theory packages for the Julia programming language. In Proceedings of the 2017 ACM International Symposium on Symbolic and Algebraic Computation , ISSAC ’17, page 157–164, 2017. URL https://doi.org/10.1145/3087604.3087611

  27. [27]

    Gawrilow and M

    E. Gawrilow and M. Joswig. Polymake: a framework for analyzing convex polytopes. In Polytopes—combinatorics and computation , pages 43–73. Springer, 2000. URL https://doi.org/10.1007/978-3-0348-8438-9_2

  28. [28]

    I. M. Gelfand, M. M. Kapranov, and A. V. Zelevinsky. A-Discriminants. Birkh¨ auser Boston, 1994. URL http://dx.doi.org/10.1007/978-0-8176-4771-1_10

  29. [29]

    V. P. Gerdt and D. Robertz. Algorithmic approach to strong consistency analysis of finite difference approximations to PDE systems. In Proceedings of the 2019 Inter- national Symposium on Symbolic and Algebraic Computation , ISSAC ’19, July 2019. URL http://dx.doi.org/10.1145/3326229.3326255

  30. [30]

    V. P. Gerdt, M. Lange-Hegermann, and D. Robertz. The MAPLE package TDDS for computing Thomas decompositions of systems of nonlinear PDEs. Computer Physics Communications, 234:202–215, 2019. URL https://doi.org/10.1016/j.cpc.2018. 07.025

  31. [31]

    Grigor’ev

    D. Grigor’ev. Complexity of quantifier elimination in the theory of ordinary dif- ferential equations , page 11–25. Springer Berlin Heidelberg, 1989. URL http: //dx.doi.org/10.1007/3-540-51517-8_81

  32. [32]

    Gustavson, A

    R. Gustavson, A. Ovchinnikov, and G. Pogudin. New order bounds in differential elimination algorithms. Journal of Symbolic Computation , 85:128–147, Mar. 2018. URL http://dx.doi.org/10.1016/j.jsc.2017.07.006

  33. [33]

    H. A. Harrington and R. A. Van Gorder. Reduction of dimension for nonlinear dynamical systems. Nonlinear Dynamics , 88(1):715–734, Dec. 2016. URL http: //dx.doi.org/10.1007/s11071-016-3272-5

  34. [34]

    Hartshorne

    R. Hartshorne. Graduate texts in mathematics. Algebraic geometry, 52, 1977

  35. [35]

    J. Heintz. Definability and fast quantifier elimination in algebraically closed fields. Theoretical Computer Science , 24(3):239–277, 1983. URL https://doi.org/10. 1016/0304-3975(83)90002-6

  36. [36]

    H. Hong, A. Ovchinnikov, G. Pogudin, and C. Yap. Global identifiability of differential models. Communications on Pure and Applied Mathematics , 73(9):1831–1879, 2020. URL https://doi.org/10.1002/cpa.21921

  37. [37]

    E. Hubert. Factorization-free decomposition algorithms in differential algebra. Jour- nal of Symbolic Computation , 29(4–5):641–662, May 2000. URL http://dx.doi. org/10.1006/jsco.1999.0344. – 36 – Projecting dynamical systems via a support bound Y. Mukhina, G. Pogudin

  38. [38]

    E. Hubert. Notes on triangular sets and triangulation-decomposition algorithms II: Differential systems. In Lecture Notes in Computer Science , pages 40–87. Springer Berlin Heidelberg, 2003. URL https://doi.org/10.1007/3-540-45084-x_2

  39. [39]

    Li, C.-M

    W. Li, C.-M. Yuan, and X.-S. Gao. Sparse differential resultant for Laurent differential polynomials. Foundations of Computational Mathematics , 15(2):451–517, Feb. 2015. URL http://dx.doi.org/10.1007/s10208-015-9249-9

  40. [40]

    Marco and J

    A. Marco and J. Martinez. Implicitization of rational surfaces by means of polynomial interpolation. Computer Aided Geometric Design , 19(5):327–344, May 2002. URL http://dx.doi.org/10.1016/S0167-8396(02)00094-8

  41. [41]

    Marker, M

    D. Marker, M. Messmer, and A. Pillay. Model theory of fields , volume 5. Cambridge University Press, 2017

  42. [42]

    McCallum and F

    S. McCallum and F. Winkler. Differential resultants. ITM Web of Conferences , 20: 01005, 2018. URL https://doi.org/10.1051/itmconf/20182001005

  43. [43]

    Meshkat, C

    N. Meshkat, C. Anderson, and J. J. DiStefano III. Alternative to Ritt’s pseudodi- vision for finding the input-output equations of multi-output models. Mathematical Biosciences, 239(1):117–123, Sept. 2012. URL http://dx.doi.org/10.1016/j.mbs. 2012.04.008

  44. [44]

    C. Moog, J. Perraud, P. Bentz, and Q. Vo. Prime differential ideals in nonlinear rational controls systems , page 17–21. 1990. URL http://dx.doi.org/10.1016/ B978-0-08-037022-4.50009-0

  45. [45]

    Ollivier

    F. Ollivier. Standard bases of differential ideals , page 304–321. Springer Berlin Hei- delberg, 1991. URL http://dx.doi.org/10.1007/3-540-54195-0_60

  46. [46]

    Oscar – open source computer algebra research system, version 1.0.0, 2024

    OSCAR. Oscar – open source computer algebra research system, version 1.0.0, 2024. URL https://www.oscar-system.org

  47. [47]

    Ovchinnikov, G

    A. Ovchinnikov, G. Pogudin, and T. N. Vo. Bounds for elimination of unknowns in systems of differential-algebraic equations. International Mathematics Research No- tices, 2022(16):12342–12377, Apr. 2021. URL http://dx.doi.org/10.1093/imrn/ rnaa302

  48. [48]

    A. Pascadi. Computer-assisted proofs of congruences for multipartitions and di- visor function convolutions, based on methods of differential algebra. The Ra- manujan Journal , 57(1):1–36, Nov. 2021. URL http://dx.doi.org/10.1007/ s11139-021-00506-8

  49. [49]

    S. L. Pimm and J. C. Rice. The dynamics of multispecies, multi-life-stage models of aquatic food webs. Theoretical Population Biology, 32(3):303–325, Dec. 1987. URL http://dx.doi.org/10.1016/0040-5809(87)90052-9

  50. [50]

    G. Pogudin. Lecture notes on differential algebra, 2023. URL http://www.lix. polytechnique.fr/Labo/Gleb.POGUDIN/files/da_notes.pdf

  51. [51]

    J. F. Ritt. Differential equations from the algebraic standpoint. American Mathemat- ical Society, 1932. ISBN 978-0821846056. URL https://archive.org/details/ differentialequa033050mbp. – 37 – Projecting dynamical systems via a support bound Y. Mukhina, G. Pogudin

  52. [52]

    D. Robertz. Formal Algorithmic Elimination for PDEs . Springer International Publishing, 2014. ISBN 9783319114453. URL http://dx.doi.org/10.1007/ 978-3-319-11445-3

  53. [53]

    K. Rose, B. Sturmfels, and S. Telen. Tropical implicitization revisited, 2023. URL https://arxiv.org/abs/2306.13015

  54. [54]

    S. L. Rueda. Linear sparse differential resultant formulas. Linear Algebra and its Ap- plications, 438(11):4296–4321, 2013. URL https://doi.org/10.1016/j.laa.2013. 01.016

  55. [55]

    S. L. Rueda. Differential elimination by differential specialization of Sylvester style matrices. Advances in Applied Mathematics , 72:4–37, Jan. 2016. URL http://dx. doi.org/10.1016/j.aam.2015.07.002

  56. [56]

    Samardzija and L

    N. Samardzija and L. Greller. Explosive route to chaos through a fractal torus in a generalized Lotka-Volterra model. Bulletin of Mathematical Biology , 50(5):465–491,

  57. [57]

    URL http://dx.doi.org/10.1016/S0092-8240(88)80003-X

  58. [58]

    Sturmfels, J

    B. Sturmfels, J. Tevelev, and J. Yu. The Newton polytope of the implicit equa- tion. Moscow Mathematical Journal , 7(2):327–346, 2007. URL https://doi.org/ 10.17323/1609-4514-2007-7-2-327-346

  59. [59]

    X. Tang, T. de Wolff, and R. Zhao. A new method for computing elimination ideals of likelihood equations. In Proceedings of the 2019 International Symposium on Symbolic and Algebraic Computation , ISSAC ’19, July 2019. URL http://dx.doi.org/10. 1145/3326229.3326241

  60. [60]

    van der Hoeven

    J. van der Hoeven. Newton’s method and FFT trading. Journal of Symbolic Compu- tation, 45(8):857–878, 2010. URL https://doi.org/10.1016/j.jsc.2010.03.005

  61. [61]

    R. A. Van Gorder. Triple mode alignment in a canonical model of the blue-sky catastrophe. Nonlinear Dynamics , 73:397–403, 2013. URL https://doi.org/10. 1007/s11071-013-0794-y

  62. [62]

    D. Wang. EPSILON: A library of software tools for polynomial elimination. In Math- ematical Software, 2002. URL https://doi.org/10.1142/9789812777171_0040

  63. [63]

    P. J. Wangersky. Lotka-Volterra population models. Annual Review of Ecology and Systematics, 9:189–218, 1978. URL http://www.jstor.org/stable/2096748

  64. [64]

    R. Zippel. Effective Polynomial Computation. Springer US, Boston, MA, 1993. ISBN 978-1-4613-6398-9 978-1-4615-3188-3. doi: 10.1007/978-1-4615-3188-3. URL http: //link.springer.com/10.1007/978-1-4615-3188-3

  65. [65]

    A. I. Zobnin. On standard bases in rings of differential polynomials. Journal of Mathematical Sciences, 135(5):3327–3335, June 2006. URL http://dx.doi.org/10. 1007/s10958-006-0161-3. – 38 –