Pith. sign in

REVIEW 25 references

Symbolic and Numerical Tools for $L_{\infty}$-Norm Calculation

T0 review · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read Symbolic methods compute L-infinity norms exactly but slowly, and they outperform numerical methods on ill-conditioned or parameterized systems.

arxiv 2505.13980 v1 pith:HI5ZEEJS submitted 2025-05-20 math.OC cs.MScs.SCcs.SYeess.SY

classification math.OCcs.MScs.SCcs.SYeess.SY
keywords computationinftysymbolicnormnumericalparticularlyanalysisapproaches
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 L-infinity norm of a linear system measures its worst-case response to inputs: the largest amplification the system can ever produce. Control engineers use it to judge stability and robustness. Standard numerical methods, like bisection on frequency grids or eigenvalue searches, usually give good answers, but they can fail when the system is nearly singular, has huge or tiny coefficients, or depends on parameters. This paper reviews a different, exact approach. Instead of sampling frequencies, symbolic methods translate the norm into a polynomial equation and then use tools from algebraic geometry: Sturm-Habicht sequences count real roots exactly, rational univariate representations solve polynomial systems, and cylindrical algebraic decomposition handles parameters. The authors demonstrate these tools on small examples and benchmark them against Matlab and Maple. They find that symbolic methods return exact, certified values and remain reliable for extremely small damping values, but are far slower: a 5 by 5 system can take minutes to hours, while numerical methods take seconds. The paper concludes that symbolic methods are best used as a complement to numerical methods, for validation and for parameterized small systems, not as a replacement for large-scale computation.
Extended reading notes

Core claim

The paper's central comparative statement is that symbolic methods are 'robust but computationally expensive' and deliver certified accuracy where numerical methods 'can yield inaccurate results and sometimes fail' (Section 4, Table 1). If true, symbolic computation is a reliable complement to numerical L-infinity norm computation for ill-conditioned or parameterized systems.

Load-bearing premise

The benchmark conclusions rest on the assumption that the reported numerical failures (e.g., Maple NormHinf returning 1.32 in Example 10, Matlab reporting infinity in Example 9) are inherent limitations of numerical methods rather than artifacts of specific implementations or default settings. If those failures are fixable by tuning the solvers, the paper's contrast between symbolic and numerical reliability is overstated.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

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

No free parameters are fitted in the paper; all quantities are either inputs or computed exactly. The paper relies on standard algebraic facts and on propositions from prior work, cited as axioms. No new physical entities are introduced.

assumptions (4)
  • standard math The maximum modulus principle applies to functions in RH_infinity, so the L-infinity norm over the right half-plane equals the supremum over the imaginary axis.
    Invoked in Section 2.1 to reduce the norm to a frequency-domain supremum; cited to standard references, not proven in the paper.
  • domain assumption Proposition 1 (Kanno and Smith 2006): gamma > ||G||_infinity if and only if gamma > sigma-bar(G(i infinity)) and det(Phi_gamma(i omega)) != 0 for all real omega.
    Basis for the symbolic reformulation in Section 3.1; accepted from the cited literature without proof.
  • domain assumption Proposition 2 (Bouzidi et al. 2021): the L-infinity norm is the maximum of the gamma-projections of the real solutions of the critical system and of the leading coefficient roots.
    Transfers L-infinity norm computation to real root finding for a bivariate polynomial system; cited from prior work.
  • standard math The specialization property of subresultants and Sturm-Habicht sequences holds under ring homomorphisms.
    Used in Section 3.2 to specialize the Sturm-Habicht sequence to parameter values; cited from Gonzalez-Vega et al. 1990/1994.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Symbolic and Numerical Tools for $L_{\infty}$-Norm Calculation." pith.science (2026). https://pith.science/paper/HI5ZEEJS

@misc{pith2026250513980,
  author       = {Pith},
  title        = {Pith review of: Symbolic and Numerical Tools for $L_\infty$-Norm Calculation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HI5ZEEJS}},
  note         = {Machine review of arXiv:2505.13980}
}
abstract

The computation of the $L_\infty $-norm is an important issue in $H_{\infty}$ control, particularly for analyzing system stability and robustness. This paper focuses on symbolic computation methods for determining the $L_{\infty} $-norm of finite-dimensional linear systems, highlighting their advantages in achieving exact solutions where numerical methods often encounter limitations. Key techniques such as Sturm-Habicht sequences, Rational Univariate Representations (RUR), and Cylindrical Algebraic Decomposition (CAD) are surveyed, with an emphasis on their theoretical foundations, practical implementations, and specific applicability to $ L_{\infty} $-norm computation. A comparative analysis is conducted between symbolic and conventional numerical approaches, underscoring scenarios in which symbolic computation provides superior accuracy, particularly in parametric cases. Benchmark evaluations reveal the strengths and limitations of both approaches, offering insights into the trade-offs involved. Finally, the discussion addresses the challenges of symbolic computation and explores future opportunities for its integration into control theory, particularly for robust and stable system analysis.

Figures

Figures reproduced from arXiv: 2505.13980 by the authors.

Figure 1
Figure 1. Regions and root functions for ξ > 0. JuliaTime column reports execution time of the sym￾bolic method via PACE.jl on Julia. The matrices and benchmark codes are available upon request. Note that the JuliaTime column reports only execution times; the computed H∞-norm values coincide exactly with those of MapSymNorm column. Size – Degree MatNumNorm MapNumNorm MapSymNorm (Time (s)) JuliaTime (s) 2×2 – 2 3.0116 3.0116 3… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

25 extracted references · 22 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter doi edition editor eid howpublished institution journal key month note number organization pages publisher school series title type url volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sent...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

  3. [3]

    (2006 a )

    Basu, S., Pollack, R., and Roy, M.F. (2006 a ). Algorithms in Real Algebraic Geometry. Springer Verlag

  4. [4]

    (2006 b )

    Basu, S., Pollack, R., and Roy, M.F. (2006 b ). Existential theory of the reals, volume 10 of algorithms and computation in mathematics

  5. [5]

    and Praagman, C

    Belur, M. and Praagman, C. (2011). An efficient algorithm for computing the h_ norm. IEEE Transactions on Automatic Control, (56), 1656--1660

  6. [6]

    Bouzidi, Y., Lazard, S., Moroz, G., Pouget, M., Rouillier, F., and Sagraloff, M. (2016). Solving bivariate systems using rational univariate representations. Journal of Complexity, 37, 34--75

  7. [7]

    Bouzidi, Y., Lazard, S., Pouget, M., and Rouillier, F. (2015). Separating linear forms and rational univariate representations of bivariate systems. Journal of Symbolic Computation, 68, 84--119

  8. [8]

    Bouzidi, Y., Quadrat, A., Rouillier, F., and Younes, G. (2021). Computation of the L _ -norm of finite-dimensional linear systems . Communications in Computer and Information Science , 1414, 119--136

Show all 25 references
  1. [9]

    Boyd, S., Balakrishnan, V., and Kabamba, P. (1989). A bisection method for computing the H _ norm of a transfer matrix and related problems. Math. Control, Signals, and Systems, (2), 207--220

  2. [10]

    Chen, C., Maza, M.M., and Xie, Y. (1990). A fast algorithm to compute the h_ -norm of a transfer function matrix. Systems and Control Letters, (14), 287--293

  3. [11]

    Chen, C., Moreno Maza, M., and Xie, Y. (2013). Computing the supremum of the real roots of a parametric univariate polynomial. In Proceedings of the 2013 International Symposium on Symbolic and Algebraic Computation (ISSAC), 141--148. ACM, Philadelphia, PA, USA

  4. [12]

    Cox, D.A., Little, J., and O'Shea, D. (2013). Ideals, varieties, and algorithms: an introduction to computational algebraic geometry and commutative algebra. Springer Science & Business Media

  5. [13]

    Cox, D.A., Little, J., and O'shea, D. (2005). Using algebraic geometry. 185

  6. [14]

    Demin, A., Rouillier, F., and Ruiz, J. (2024). Reading rational univariate representations on lexicographic groebner bases. arXiv preprint arXiv:2402.07141

  7. [15]

    Dominguez, L.F., Narciso, D.A., and Pistikopoulos, E.N. (2010). Recent advances in multiparametric nonlinear programming. Computers & Chemical Engineering, 34(5), 707--716

  8. [16]

    Doyle, J.C., Francis, B.A., and Tannenbaum, A.R. (2013). Feedback Control Theory. Courier Corporation

  9. [17]

    Gonz \'a lez-Vega, L., Lombardi, H., Recio, T., and Roy, M.F. (1990). Sp \'e cialisation de la suite de sturm et sous-r \'e sultants (i). RAIRO-Theoretical Informatics and Applications, 24(6), 561--588

  10. [18]

    Gonz \'a lez-Vega, L., Lombardi, H., Recio, T., and Roy, M.F. (1994). Sp \'e cialisation de la suite de sturm. RAIRO-Theoretical Informatics and Applications, 28(1), 1--24

  11. [19]

    Gonz \'a lez-Vega, L., Recio, T., Lombardi, H., and Roy, M.F. (1998). Sturm- Habicht sequences, determinants and real roots of univariate polynomials. In Quantifier Elimination and Cylindrical Algebraic Decomposition, 300--316. Springer

  12. [20]

    Kahoui, M.E. (2003). An elementary approach to subresultants theory. Journal of Symbolic Computation, 35(3), 281--292

  13. [21]

    and Smith, M.C

    Kanno, M. and Smith, M.C. (2006). Validated numerical computation of the L _ -norm for linear dynamical systems. Journal of Symbolic Computation, 41(6), 697--707

  14. [22]

    Quadrat, A., Rouillier, F., and Younes, G. (2024). L ^ -norm computation for linear time-invariant systems depending on parameters . Maple Transactions , 4(1), 18

  15. [23]

    Rouillier, F. (1999). Solving zero-dimensional systems through the rational univariate representation. Applicable Algebra in Engineering, Communication and Computing, 9(5), 433--461

  16. [24]

    and Zimmermann, P

    Rouillier, F. and Zimmermann, P. (2004). Efficient isolation of polynomial's real roots. Journal of Computational and Applied Mathematics, 162(1), 33--50

  17. [25]

    Zhou, K., Doyle, J.C., and Glover, K. (1996). Robust and Optimal Control. Prentice Hall

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.