History of confluent Vandermonde matrices and inverting them algorithms
Pith reviewed 2026-05-23 22:45 UTC · model grok-4.3
The pith
A numerical algorithm inverts confluent Vandermonde matrices in quadratic time for any allowed parameters including large root multiplicities.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper establishes the historical record of confluent Vandermonde matrices since 1891 and supplies a ready-to-implement numerical procedure that inverts any such matrix in quadratic time for arbitrary parameter values allowed by the definition, including high multiplicities, while avoiding symbolic computations entirely.
What carries the argument
The numerical inversion algorithm that operates directly on matrix entries without symbolic operations or special-case handling for multiplicities.
If this is right
- The algorithm can be coded directly in any general-purpose programming language or mathematical package without needing computer algebra support.
- Linear systems involving these matrices become solvable in quadratic time even when root multiplicities are large.
- Historical context for the matrix class is now collected in one place for researchers needing background.
Where Pith is reading between the lines
- The same numerical approach might be adapted to invert other structured matrices that arise from repeated roots in polynomials.
- Applications in polynomial interpolation with multiple nodes could gain a direct quadratic-time inversion step.
- The history survey may prompt similar reviews for related classes of confluent or generalized Vandermonde matrices.
Load-bearing premise
The surveyed historical facts are accurate and the presented numerical algorithm produces correct inverses for every permitted set of parameters including high multiplicities.
What would settle it
Run the algorithm on a confluent Vandermonde matrix whose characteristic polynomial has a root of multiplicity 10 or higher, then compare the computed inverse against an independently verified result while measuring that runtime grows quadratically with matrix size.
read the original abstract
The author was encouraged to write this review by numerous enquiries from researchers all over the world, who needed a ready-to-use algorithm for the inversion of confluent Vandermonde matrices which works in quadratic time for any values of the parameters allowed by the definition, including the case of large root multiplicities of the characteristic polynomial. Article gives the history of the title special matrix since 1891 and surveys algorithms for solving linear systems with the title class matrix and inverting it. In particular, it presents, also by example, a numerical algorithm which does not use symbolic computations and is ready to be implemented in a general-purpose programming language or in a specific mathematical package.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript surveys the history of confluent Vandermonde matrices since 1891 and reviews algorithms for solving linear systems and inverting such matrices. It presents a numerical algorithm for inversion claimed to run in quadratic time for any allowed parameter values (including large root multiplicities) without symbolic computation, illustrated by example.
Significance. If the algorithm's claimed quadratic-time performance and correctness for arbitrary multiplicities hold under verification, the work would supply a practical, ready-to-implement tool addressing documented researcher needs in numerical linear algebra.
major comments (2)
- [Algorithm presentation] Algorithm presentation (the section introducing the new numerical method): the quadratic-time claim for arbitrary multiplicities rests only on an illustrative example; no general complexity analysis, inductive argument, or matrix-algebraic derivation establishing O(n²) runtime or correctness for high multiplicities is supplied.
- [Abstract and algorithm section] Abstract and algorithm section: no verification data, test cases with explicit multiplicities, or implementation pseudocode are provided, preventing assessment of functionality beyond the single example.
Simulated Author's Rebuttal
We thank the referee for the thoughtful comments on our historical survey and the presented algorithm. We address the major comments point by point below, agreeing where revisions are needed to strengthen the work.
read point-by-point responses
-
Referee: [Algorithm presentation] Algorithm presentation (the section introducing the new numerical method): the quadratic-time claim for arbitrary multiplicities rests only on an illustrative example; no general complexity analysis, inductive argument, or matrix-algebraic derivation establishing O(n²) runtime or correctness for high multiplicities is supplied.
Authors: We acknowledge that the manuscript supports the quadratic-time claim primarily through the structure and operations shown in the detailed illustrative example rather than a standalone formal analysis. The algorithm proceeds via block-wise constructions tied to the multiplicities, with the total arithmetic operations scaling quadratically with matrix dimension n as the sum of multiplicity-related terms remains O(n²). To address the concern directly, the revised manuscript will include an explicit general complexity subsection deriving the O(n²) bound from the operation counts in the recursive inverse construction, along with a brief argument for correctness across arbitrary multiplicities. revision: yes
-
Referee: [Abstract and algorithm section] Abstract and algorithm section: no verification data, test cases with explicit multiplicities, or implementation pseudocode are provided, preventing assessment of functionality beyond the single example.
Authors: The manuscript presents the algorithm via one concrete example with chosen multiplicities to show its numerical, non-symbolic nature and readiness for implementation. We agree this limits independent verification. In revision we will add explicit pseudocode for the full procedure, plus additional test cases covering a range of multiplicity configurations (including high multiplicities), with reported numerical results to confirm functionality. revision: yes
Circularity Check
No circularity; historical survey with example-based algorithm presentation
full rationale
The manuscript is a review of the history of confluent Vandermonde matrices since 1891 together with a survey of existing algorithms for their inversion. It presents one numerical algorithm 'also by example' without any derivation chain, first-principles result, fitted parameter, or prediction that reduces to its own inputs. No load-bearing self-citation, uniqueness theorem, or ansatz is invoked to justify a new claim. The quadratic-time statement for arbitrary multiplicities is asserted on the basis of the exhibited procedure rather than derived from prior equations within the paper, so no circular reduction exists. This is the normal non-finding for a survey paper whose central content is archival and illustrative rather than deductive.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
- [1]
-
[2]
A.C. Ahlin, A bivariate generalization of Hermite’s interpolation formula, Mathematics of Computation 18 (1964) 264-273
work page 1964
-
[3]
Aitken, XII—Further Numerical Studies in Algebraic Equations and Matrices, Proc
A.C. Aitken, XII—Further Numerical Studies in Algebraic Equations and Matrices, Proc. Royal Society of Edinburgh 51 (1932) 80-90
work page 1932
-
[4]
Aitken, XX—Studies in Practical Mathematics
A.C. Aitken, XX—Studies in Practical Mathematics. II. The Evaluation of the Latent Roots and Latent Vectors of a Matrix, Proc. Royal Society of Edinburgh 57 (1938) 269-304
work page 1938
-
[5]
Aitken, Determinants and Matrices, Third Edition, Interscience Publishers, New York, USA, 1944
A.C. Aitken, Determinants and Matrices, Third Edition, Interscience Publishers, New York, USA, 1944
work page 1944
-
[6]
Bellman, Introduction to Matrix Analysis, 2nd
R. Bellman, Introduction to Matrix Analysis, 2nd. ed., SIAM, Philadelphia, USA, 1997
work page 1997
- [7]
- [8]
-
[9]
A. Cayley, A Memoir on the Theory of Matrices, Philosophical Transactions of the Royal Society of London 148 (1857) 17-37
-
[10]
F-Ch. Chang, The inverse of the generalized Vandermonde matrix through the partial fraction expansion, IEEE Transactions on Automatic Control 19 (2) (1974) 151-152
work page 1974
- [11]
-
[12]
D. Coppersmith, S. Winograd, Matrix multiplication via Arithmetic Progressions, J. Symbolic Computation 9 (1990) 251-280
work page 1990
-
[13]
F. Csaki, Some notes on the inversion of confluent Vandermonde matrices, IEEE Transactions on Automatic Control 20 (1) (1975) 154-157
work page 1975
-
[14]
Cullis, Matrices and Determinoids, Vol
C.E. Cullis, Matrices and Determinoids, Vol. 1-3, Cambridge University Press, Cambridge, UK, 1913, 1918, 1925
work page 1913
-
[15]
Mac Duffee, The Theory of Matrices, Springer Verlag, Berlin, Germany, 1933
C.C. Mac Duffee, The Theory of Matrices, Springer Verlag, Berlin, Germany, 1933
work page 1933
-
[16]
R.A. Frazer, W.J. Duncan, A.R. Collar, Elementary matrices and some applications to dynamics and differential equations, Cambridge University Press, Cambridge, UK, 1938
work page 1938
-
[17]
Fuchs, The explicit inverse of the stiffness matrix, Int
M.B. Fuchs, The explicit inverse of the stiffness matrix, Int. J. Solids Struct. 29 (1992) 2101-2113
work page 1992
-
[18]
G. Galimberti, V. Pereyra, Solving confluent Vandermonde systems of Hermite type, Numerische Mathematik 18 (1971) 44–60. J.S. Respondek History of confluent Vandermonde matrices … Page 14 of 20
work page 1971
-
[19]
F.R. Gantmacher, Applications of the Theory of Matrices, Interscience Publishers, Inc., New York, USA, 1959
work page 1959
-
[20]
Gantmacher, The Theory of Matrices, Vol
F.R. Gantmacher, The Theory of Matrices, Vol. I-II, Chelsea Publishing Company, New York, USA, 1960
work page 1960
-
[21]
Gautschi, On inverses of Vandermonde and confluent Vandermonde matrices
W. Gautschi, On inverses of Vandermonde and confluent Vandermonde matrices. II, Numerische Mathematik 5 (1963) 425-430
work page 1963
-
[22]
I. Goknar, Obtaining the inverse of the generalized Vandermonde matrix of the most general type, IEEE Transactions on Automatic Control 18 (5) (1973) 530-532
work page 1973
-
[23]
G.H. Golub, Ch.F. Van Loan, Matrix Computations, Fourth Edition, The Johns Hopkins University Press, Baltimore, USA, 2013
work page 2013
-
[24]
L.A. González-Serrano, E.A. Maximenko, Bialternant formula for Schur polynomials with repeating variables, arXiv:2312.15680, 2023 (submitted for publication)
-
[25]
H. Gorecki, On switching instants in minimum-time control problem, One-dimensional case n-tuple eigenvalue, Bull. de L’Acad. Pol. Des. Sci. 16 (1968) 23–30
work page 1968
-
[26]
T.T. Ha, J.A. Gibson, A note on the determinant of a functional confluent Vandermonde matrix and controllability, Linear Algebra and its Applications 30 (1980) 69-75
work page 1980
-
[27]
N. Halidias, Computing the Minimum Polynomial, the Function and the Drazin Inverse of a Matrix with Matlab, Asian Journal of Research in Computer Science 17 (5) (2024) 1-9
work page 2024
-
[28]
W.A. Harris, J.P. Fillmore, D.R. Smith, Matrix Exponentials-Another Approach, SIAM Review 43 (4) (2001) 694-706
work page 2001
-
[29]
Hawkins, The Theory of Matrices in the 19th Century, Proc
T. Hawkins, The Theory of Matrices in the 19th Century, Proc. Int. Congress of Mathematicians, Vancouver, Vol. 2 (1974) 561-570
work page 1974
-
[30]
M.Ch. Hermite, M. Borchardt, Sur la formule d'interpolation de Lagrange, Journal für die reine und angewandte Mathematik 84 (1878) 70-79
-
[31]
Higham, Functions of Matrices, Theory and Computation, SIAM, Philadelphia, USA, 2008
N.J. Higham, Functions of Matrices, Theory and Computation, SIAM, Philadelphia, USA, 2008
work page 2008
-
[32]
R.A. Horn, Ch.R. Johnson, Matrix Analysis, Second Edition, Cambridge University Press, New York, USA, 2013
work page 2013
-
[33]
S-H. Hou, E.S-H. Hou, A Recursive Algorithm for Triangular Factorization of Inverse of Confluent Vandermonde Matrices, AIP Conf. Proc. Vol. 1089 (1) (2009) 277–288
work page 2009
- [34]
-
[35]
A.S. Householder, The Theory of Matrices in Numerical Analysis, First edition, Blaisdell Publishing Company, New York, USA, 1964
work page 1964
-
[36]
T. Kaczorek, Vectors and matrices in automation and electrical engineering, 2nd ed., WNT, Warsaw, Poland, 1998 (in Polish). J.S. Respondek History of confluent Vandermonde matrices … Page 15 of 20
work page 1998
-
[37]
Kalman, The Generalized Vandermonde Matrix, Mathematics Magazine 57 (1) (1984) 15-21
D. Kalman, The Generalized Vandermonde Matrix, Mathematics Magazine 57 (1) (1984) 15-21
work page 1984
-
[38]
I. Kaufman, The inversion of the Vandermonde matrix and transformation to the Jordan canonical form, IEEE Transactions on Automatic Control 14 (6) (1969) 774-777
work page 1969
- [39]
-
[40]
G. Kalogeropoulos, P. Psarrakos, A note on the controllability of higher-order linear systems, Applied Mathematics Letters 17 (12) (2004) 1375-1380
work page 2004
- [41]
-
[42]
R-C. Li, Lower Bounds for the Condition Number of a Real Confluent Vandermonde Matrix, Mathematics of Computation 75 (256) (2006) 1987–1995
work page 2006
-
[43]
K.A. Linsay, C.E. Rooney, A note on compound matrices, J. Computat. Phys. 103 (1992) 472-477
work page 1992
-
[44]
G.G. Lorentz, K. Jetter, S.D. Riemenschneider, Birkhoff Interpolation, 1st ed., Encyclopedia of Mathematics and its Applications Book 19, Cambridge University Press, UK, 1984
work page 1984
-
[45]
Lu, Fast Solution of Confluent Vandermonde Linear Systems, SIAM J
H. Lu, Fast Solution of Confluent Vandermonde Linear Systems, SIAM J. Matrix Anal. Appl. 15 (4) (1994) 1277-1289
work page 1994
-
[46]
H. Lu, Fast Algorithms For Confluent Vandermonde Linear Systems and Generalized Trummers Problem, SIAM J. Matrix Anal. Appl. 16 (2) (1995) 655-674
work page 1995
-
[47]
L. Lupas, On the computation of the generalized Vandermonde matrix inverse, IEEE Transactions on Automatic Control 20 (4) (1975) 559-561
work page 1975
- [48]
-
[49]
Malik, Compound matrices to the tree-generating problem, IEEE Trans
R.N. Malik, Compound matrices to the tree-generating problem, IEEE Trans. Circuit Theory 17 (1970) 149-151
work page 1970
-
[50]
M.E.A. El-Mikkawy, Explicit inverse of a generalized Vandermonde matrix, Applied Mathematics and Computation 146 (2003) 643–651
work page 2003
- [51]
-
[52]
Muir, A treatise on the theory of determinants, MacMillan and co., London, UK, 1882
T. Muir, A treatise on the theory of determinants, MacMillan and co., London, UK, 1882
-
[53]
Muir, The Theory of Determinants in the Historical Order of Its Development
T. Muir, The Theory of Determinants in the Historical Order of Its Development. Part I. Determinants in General. Leibnitz (1693) to Cayley (1841), MacMillan and co., London, UK, 1890
-
[54]
Muir, The Theory of Determinants in the Historical Order of Development
T. Muir, The Theory of Determinants in the Historical Order of Development. Volumes Two to Four. The Periods 1841 to 1860, 1861 to 1880, 1880 to 1900, Dover Publications, Inc., New York, USA, 1911, 1920, 1923. J.S. Respondek History of confluent Vandermonde matrices … Page 16 of 20
work page 1900
-
[55]
Muir, Contributions to the History of Determinants
T. Muir, Contributions to the History of Determinants. 1900-1920, Blackie & Son Limited, London and Glasgow, UK, 1930
work page 1900
- [56]
-
[57]
Nambiar, J.D Keating, Application of compound matrices to linear systems, IEEE Trans
K.K. Nambiar, J.D Keating, Application of compound matrices to linear systems, IEEE Trans. Circuit Theory 17 (1970) 626-628
work page 1970
-
[58]
Nambiar, Hall's theorem and compound matrices, Math
K.K. Nambiar, Hall's theorem and compound matrices, Math. Comput. Modelling 25 (1997) 23-24
work page 1997
-
[59]
V. Pan, Structured Matrices and Polynomials Unified Superfast Algorithms, Springer Science+Business Media, LLC, New York, USA, 2001
work page 2001
-
[60]
Pan, How Bad Are Vandermonde Matrices? SIAM J
V. Pan, How Bad Are Vandermonde Matrices? SIAM J. Matrix Anal. Appl. 37 (2) (2016) 676-694
work page 2016
-
[61]
U. Prells, M.I. Friswell, S.D. Garvey, Use of geometric algebra: compound matrices and the determinant of the sum of two matrices, Proc. Royal Society 459 (2003) 273-285
work page 2003
-
[62]
J.S. Respondek, Approximate controllability of the n-th order infinite dimensional systems with controls delayed by the control devices, Int. J. Syst. Sci. 39 (8) (2008) 765–782
work page 2008
-
[63]
Respondek, On the confluent Vandermonde matrix calculation algorithm, Appl
J.S. Respondek, On the confluent Vandermonde matrix calculation algorithm, Appl. Math. Lett. 24 (2011) 103–106
work page 2011
-
[64]
J.S. Respondek, Numerical recipes for the high efficient inverse of the confluent Vandermonde matrices, Appl. Math. Comput. 218 (2011) 2044–2054
work page 2011
-
[65]
R. Schappelle, The inverse of the confluent Vandermonde matrix, IEEE Transactions on Automatic Control 17 (5) (1972) 724-725
work page 1972
-
[66]
Schendel, Das alternirende Exponentialdifferenzenproduct, Zeitschrift Math
L. Schendel, Das alternirende Exponentialdifferenzenproduct, Zeitschrift Math. Phys. (1891) 84-94
-
[67]
Scott, Theory of Determinants, Cambridge University Press, Cambridge, UK, 1880
R.F. Scott, Theory of Determinants, Cambridge University Press, Cambridge, UK, 1880
-
[68]
B. Shen, H. Tan, Z. Wang, T. Huang, Quantized/Saturated Control for Sampled-Data Systems Under Noisy Sampling Intervals: A Confluent Vandermonde Matrix Approach, IEEE Transactions on Automatic Control 62 (9) (2017) 4753-4759
work page 2017
-
[69]
A. Spitzbart, A Generalization of Hermite's Interpolation Formula, The American Mathematical Monthly 67 (1) (1960) 42-46
work page 1960
-
[70]
J.J. Sylvester, Additions to the articles in the September number of this journal, "On a new class of theorems," and on Pascal's theorem, Phil. Mag. 3 (37) (1850) 363-370
-
[71]
J.J. Sylvester, On the Relation between the Minor Determinants of Linearly Equivalent Quadratic Functions, Phil. Mag. 4 (1) (1851) 295-305
-
[72]
T-Y Tam, X. Liu, Matrix Inequalities and Their Extensions to Lie Groups, 1st Edition, Chapman & Hall/CRC Monographs and Research Notes in Mathematics, CRC Press, Taylor & Francis Group, New York, USA, 2018. J.S. Respondek History of confluent Vandermonde matrices … Page 17 of 20
work page 2018
- [73]
-
[74]
H.W. Turnbull, The Theory of Determinants, Matrices and Invariants, Blackie & Son, London & Glasgow, UK, 1928
work page 1928
-
[75]
H.W. Turnbull, A.C. Aitken, An Introduction to the Theory of Canonical Matrices, London, Glasgow and Bombay: Blackie and Son, 1932
work page 1932
-
[76]
Vogt, Sur l’apolarité des formes binaires
M. Vogt, Sur l’apolarité des formes binaires. Nouvelles annales de mathématiques 4 (1) (1901) 337-365
work page 1901
-
[77]
Wedderburn, Lectures on Matrices
J.H.M. Wedderburn, Lectures on Matrices. American Mathematical Society, Colloquium Publications, Providence, Rhode Island, USA, 1934
work page 1934
- [78]
-
[79]
X. Zhong, Y. Zhaoyong, A Fast Algorithm for Inversion of Confluent Vandermonde-Like Matrices Involving Polynomials That Satisfy a Three-Term Recurrence Relation, SIAM J. Matrix Anal. Appl. 19 (3) (1998) 797-806. J.S. Respondek History of confluent Vandermonde matrices … Page 18 of 20 Appendix – example of execution of the algorithm Let us invert the matri...
work page 1998
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.