Local Linearizations of Rational Matrices with Application to Rational Approximations of Nonlinear Eigenvalue Problems
Pith reviewed 2026-05-24 16:17 UTC · model grok-4.3
The pith
A unified definition of local linearizations for rational matrices recovers every prior pencil as a special case while preserving zero-pole structure in arbitrary subsets of any algebraically closed field and at infinity.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
A definition of local linearization is introduced that produces matrix pencils associated to a given rational matrix; these pencils preserve the zero-pole structure of the rational matrix inside any prescribed subset of an algebraically closed field together with the point at infinity, and every pencil previously studied in the literature arises as a special case of this definition.
What carries the argument
The definition of a local linearization, which is a matrix pencil tied to a rational matrix such that the pencil and the rational matrix share the same zeros and poles inside any chosen subset of an algebraically closed field and at infinity.
If this is right
- Every pencil used since the 1970s for computing zeros and poles of rational matrices is now a special case inside one theory.
- Pencils arising from rational approximations of nonlinear eigenvalue problems inherit the zero-pole preservation property automatically.
- The same pencil can be analyzed for structure preservation at infinity and inside any finite subset of the field without switching definitions.
- Numerical algorithms that rely on these pencils gain a uniform justification for their correctness across different choices of linearization.
Where Pith is reading between the lines
- Implementations of nonlinear eigenvalue solvers could switch among historically different pencils inside a single code base without loss of theoretical guarantees.
- The framework suggests a systematic way to construct new pencils that preserve additional structure, such as symmetry or palindromicity, on chosen subsets.
- Questions about the minimal size or conditioning of pencils for a given rational matrix can now be posed uniformly rather than case by case.
Load-bearing premise
It is possible to write one definition of local linearization that simultaneously preserves zero-pole locations for every possible subset of an algebraically closed field and recovers every earlier pencil construction as a special case.
What would settle it
Existence of a pencil from the 1970-2019 literature that cannot be obtained as a special case of the new definition, or a rational matrix for which a local linearization fails to match the zeros or poles of the original matrix inside some subset.
read the original abstract
This paper presents a definition for local linearizations of rational matrices and studies their properties. This definition allows us to introduce matrix pencils associated to a rational matrix that preserve its structure of zeros and poles in subsets of any algebraically closed field and also at infinity. Moreover, such definition includes, as particular cases, other definitions that have been used previously in the literature. In this way, this new theory of local linearizations captures and explains rigorously the properties of all the different pencils that have been used from the 1970's until 2019 for computing zeros, poles and eigenvalues of rational matrices. Particular attention is paid to those pencils that have appeared recently in the numerical solution of nonlinear eigenvalue problems through rational approximation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces a definition of local linearizations for rational matrices. This definition constructs associated matrix pencils that preserve the zero-pole structure of the rational matrix on arbitrary subsets of any algebraically closed field (including at infinity) and recovers as special cases all previously used pencils in the literature from the 1970s onward. The framework is applied to pencils arising in rational approximations for nonlinear eigenvalue problems.
Significance. If the unification holds, the work supplies a single rigorous foundation that explains the zero-pole preservation properties of disparate pencil constructions used for rational and nonlinear eigenproblems. This could standardize analysis and algorithm design in numerical linear algebra, particularly where rational approximations are employed.
major comments (2)
- [§3] §3, Definition 3.2: the statement that the new definition recovers all prior pencils as special cases requires an explicit bijection or parameter choice for each historical construction (e.g., the pencils of Gohberg et al. and the recent rational-approximation pencils); without this mapping the unification claim is not fully substantiated.
- [§4.2] §4.2, Theorem 4.5: the proof that zeros and poles are preserved on an arbitrary subset S assumes that the local linearization is regular on S; the argument does not address the case when the subset intersects the pole set of the original rational matrix, which is load-bearing for the general claim.
minor comments (2)
- [§2] Notation for the field K and its algebraic closure is introduced late; moving the definitions to §2 would improve readability.
- Figure 1 caption does not state the field or the subset S used in the example.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments. Both points identify opportunities to strengthen the explicitness of the unification and the generality of the main preservation result. We address each below and will incorporate the suggested clarifications in the revised manuscript.
read point-by-point responses
-
Referee: §3, Definition 3.2: the statement that the new definition recovers all prior pencils as special cases requires an explicit bijection or parameter choice for each historical construction (e.g., the pencils of Gohberg et al. and the recent rational-approximation pencils); without this mapping the unification claim is not fully substantiated.
Authors: We agree that the unification claim is more convincingly substantiated when explicit parameter choices are supplied. In the revision we will add a dedicated subsection (or short appendix) that lists the precise choices of the local linearization parameters recovering the classical pencils of Gohberg, Lancaster and Rodman as well as the recent pencils arising from rational approximations of nonlinear eigenvalue problems. These mappings will be stated as corollaries to Definition 3.2. revision: yes
-
Referee: §4.2, Theorem 4.5: the proof that zeros and poles are preserved on an arbitrary subset S assumes that the local linearization is regular on S; the argument does not address the case when the subset intersects the pole set of the original rational matrix, which is load-bearing for the general claim.
Authors: The referee correctly identifies that the present proof of Theorem 4.5 is written under the standing assumption that the linearization is regular on S and does not separately treat the case in which S meets the pole set of the rational matrix. We will revise the statement of Theorem 4.5 to make the regularity assumption explicit and extend the proof to cover the intersecting-pole case, either by a direct argument or by reducing it to the regular situation via a suitable deflation at the poles. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper introduces a new definition of local linearizations for rational matrices that is shown to preserve zero-pole structure on arbitrary subsets of an algebraically closed field (including infinity) while recovering prior pencils as special cases. This unification follows directly from the definition and its derived properties rather than reducing to any fitted parameter, self-citation chain, or renaming of known results. No load-bearing step in the derivation is equivalent to its inputs by construction, and the central claim remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- standard math Rational matrices are considered over algebraically closed fields
invented entities (1)
-
Local linearization
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanabsolute_floor_iff_bare_distinguishability unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Definition 4.3 (Linearization in a subset of F) and Theorem 3.7 (elementary divisors of minimal polynomial system matrix recover pole/zero divisors in Σ)
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Use of g-reversals and local equivalence at ∞ (Definition 3.9, Proposition 4.17)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
R. Alam, N. Behera, Linearizations for rational matrix functions and Rosenbroc k sys- tem polynomials, SIAM J. Matrix Anal. Appl. 37(1) (2016) 354–380
work page 2016
-
[2]
A. Amiraslani, R. M. Corless, P. Lancaster, Linearization of matrix polynomials ex- pressed in polynomial bases , IMA J. Numer. Anal. 29 (2009) 141–157
work page 2009
-
[3]
A. Amparan, S. Marcaida, I. Zaballa, On the structure invariants of proper rational matrices with prescribed finite poles, Linear and Multilinear Algebra 61(11) (2013) 1464– 1486
work page 2013
-
[4]
A. Amparan, S. Marcaida, I. Zaballa, Finite and infinite structures of rational matrices: a local approach, Electron. J. Linear Algebra 30 (2015) 196–226
work page 2015
-
[5]
A. Amparan, F. M. Dopico, S. Marcaida, I. Zaballa, Strong linearizations of rational matrices, SIAM J. Matrix Anal. Appl. 39(4) (2018) 1670–1700
work page 2018
-
[6]
D. J. Cullen, Local system equivalence, Math. Systems Theory 19 (1986) 67-78
work page 1986
-
[7]
R. Das, R. Alam, Recovery of minimal bases and minimal indices of rational ma trices from Fiedler-like pencils , Linear Algebra Appl. 566 (2019) 34–60. 48
work page 2019
-
[8]
R. Das, R. Alam, Affine spaces of strong linearizations for rational matrices and the recovery of eigenvectors and minimal bases , Linear Algebra Appl. 569 (2019) 335–368
work page 2019
-
[9]
F. De Ter´ an, F. M. Dopico, D. S. Mackey, Spectral equivalence of matrix polynomials and the index sum theorem , Linear Algebra Appl. 459 (2014) 264–333
work page 2014
-
[10]
F. M. Dopico, P. W. Lawrence, J. P´ erez, P. Van Dooren, Block Kronecker lineariza- tions of matrix polynomials and their backward errors , Numer. Math. 140 (2018) 373– 426
work page 2018
-
[11]
F. M. Dopico, S. Marcaida, M. C. Quintana, Strong linearizations of rational matrices with polynomial part expressed in an orthogonal basis , Linear Algebra Appl. 570 (2019) 1–45
work page 2019
-
[12]
F. M. Dopico, S. Marcaida, M. C. Quintana, P. Van Dooren, Block full rank lin- earizations of rational matrices with application to ratio nal approximations of nonlinear eigenvalue problems. In preparation
- [13]
-
[14]
G. D. Forney, Jr., Minimal bases of rational vector spaces, with applications to mul- tivariable linear systems, SIAM J. Control 13(3) (1975) 493–520
work page 1975
-
[15]
F. R. Gantmacher, The Theory of Matrices. Vols. 1, 2 , Chelsea Publishing Co., New York, 1959
work page 1959
-
[16]
I. Gohberg, P. Lancaster, L. Rodman, Matrix Polynomials, Academic Press, 1982
work page 1982
-
[17]
S. G¨ uttel, F. Tisseur,The nonlinear eigenvalue problem , Acta Numer. 26 (2017) 1–94
work page 2017
-
[18]
S. G¨ uttel, R. Van Beeumen, K. Meerbergen, W. Michiels, NLEIGS: A class of fully rational Krylov methods for nonlinear eigenvalue problems , SIAM J. Sci. Comput. 36(6) (2014) A2842–A2864
work page 2014
-
[19]
Kailath, Linear Systems , Prentice Hall, New Jersey, 1980
T. Kailath, Linear Systems , Prentice Hall, New Jersey, 1980
work page 1980
-
[20]
N. P. Karampetakis, S. Vologiannidis, Infinite elementary divisors structure - pre- serving transformations for polynomial matrices , International J. Appl. Math. Comput. Sci. 13 (2003) 493–503
work page 2003
-
[21]
Automatic rational approximation and linearization of nonlinear eigenvalue problems
P. Lietaert, J. P´ erez, B. Vandereycken, K. Meerbergen , Automatic rational approx- imation and linearization of nonlinear eigenvalue problem s, submitted. Available as arXiv:1801.08622v2
work page internal anchor Pith review Pith/arXiv arXiv
-
[22]
D. Lu, X. Huang, Z. Bai, Y. Su, A Pad´ e approximate linearization algorithm for solving the quadratic eigenvalue problem with low-rank dam ping, Int. J. Numer. Meth. Engng. 103 (2015) 840–858
work page 2015
-
[23]
McMillan, Introduction to formal realizability theory I, Bell System Tech
B. McMillan, Introduction to formal realizability theory I, Bell System Tech. J. 31 (1952) 217–279. 49
work page 1952
-
[24]
McMillan, Introduction to formal realizability theory II, Bell System Tech
B. McMillan, Introduction to formal realizability theory II, Bell System Tech. J. 31 (1952) 541–600
work page 1952
-
[25]
V. Mehrmann, H. Voss, Nonlinear eigenvalue problems: A challenge for modern eige n- value methods, GAMM–Mitt. 27 (2004) 121–152
work page 2004
-
[26]
C. B. Moler, G. W. Stewart, An algorithm for generalized matrix eigenvalue problems, SIAM J. Numer. Anal. 10 (1973) 241–256
work page 1973
-
[27]
Newman, Integral Matrices, Academic Press, New York and London, 1972
M. Newman, Integral Matrices, Academic Press, New York and London, 1972
work page 1972
-
[28]
H. H. Rosenbrock, State-space and Multivariable Theory, Thomas Nelson and Sons, London, 1970
work page 1970
- [29]
-
[30]
Y. Su, Z. Bai, Solving rational eigenvalue problems via linearization , SIAM J. Matrix Anal. Appl. 32 (1) (2011) 201–216
work page 2011
-
[31]
R. Van Beeumen, O. Marques, E. G. Ng, C. Yang, Z. Bai, L. Ge , O. Kononenko, Z. Li, C.-K. Ng, L. Xiao, Computing resonant modes of accelerator cavities by solving nonlinear eigenvalue problems via rational approximation , J. Comput. Phys. 374 (2018) 1031–1043
work page 2018
-
[32]
Van Dooren, The computation of Kronecker’s canonical form of a singular p encil, Linear Algebra Appl
P. Van Dooren, The computation of Kronecker’s canonical form of a singular p encil, Linear Algebra Appl. 27 (1979) 103–140
work page 1979
-
[33]
Van Dooren, The generalized eigenstructure problem in linear system the ory, IEEE Trans
P. Van Dooren, The generalized eigenstructure problem in linear system the ory, IEEE Trans. Automat. Control 26 (1981) 111–129
work page 1981
-
[34]
P. Van Dooren, P. Dewilde, J. Vandewalle, On the determination of the Smith- McMillan form of a rational matrix from its Laurent expansio n, IEEE Trans. Circuit Syst. 26(3) (1979) 180–189
work page 1979
-
[35]
A. I. G. Vardulakis, Linear Multivariable Control, John Wiley and Sons, New York, 1991
work page 1991
-
[36]
Verghese, Comments on ‘Properties of the system matrix of a generalized state- space system’, Int
G. Verghese, Comments on ‘Properties of the system matrix of a generalized state- space system’, Int. J. Control 31(5) (1980) 1007–1009
work page 1980
-
[37]
G. Verghese, P. Van Dooren, T. Kailath, Properties of the system matrix of a gener- alized state-space system , Int. J. Control 30(2) (1979) 235–243
work page 1979
-
[38]
Vidyasagar, Control System Synthesis: A Factorization Approach , The MIT Press, Boston, 1985
M. Vidyasagar, Control System Synthesis: A Factorization Approach , The MIT Press, Boston, 1985
work page 1985
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.