A Randomized Homotopy for the Hermitian Eigenpair Problem

We describe and analyze a randomized homotopy algorithm for the Hermitian eigenvalue problem. Given an $$n\times n$$n×n Hermitian matrix $$A$$A, the algorithm returns, almost surely, a pair $$(\lambda ,v)$$(λ,v) which approximates, in a very strong sense, an eigenpair of $$A$$A. We prove that the expected cost of this algorithm, where the expectation is both over the random choices of the algorithm and a probability distribution on the input matrix $$A$$A, is $$\mathcal{{O}}(n^6)$$O(n6), that is, cubic on the input size. Our result relies on a cost assumption for some pseudorandom number generators whose rationale is argued by us.

[1]  Diego Armentano,et al.  Stochastic perturbations and smooth condition numbers , 2010, J. Complex..

[2]  S. Smale Newton’s Method Estimates from Data at One Point , 1986 .

[3]  Michael Shub,et al.  Complexity of Bezout’s Theorem VI: Geodesics in the Condition (Number) Metric , 2007, Found. Comput. Math..

[4]  Sheehan Olver,et al.  Fast inverse transform sampling in one and two dimensions , 2013, 1307.1223.

[5]  K. P. Choi On the medians of gamma distributions and an equation of Ramanujan , 1994 .

[6]  Carlos Beltrán,et al.  Fast Linear Homotopy to Find Approximate Zeros of Polynomial Systems , 2011, Found. Comput. Math..

[7]  S. Smale,et al.  Complexity of Bézout’s theorem. I. Geometric aspects , 1993 .

[8]  James Demmel,et al.  Applied Numerical Linear Algebra , 1997 .

[9]  S. Smale Mathematical problems for the next century , 1998 .

[10]  Percy Deift,et al.  Some Open Problems in Random Matrix Theory and the Theory of Integrable Systems. II , 2007, 1703.04931.

[11]  L. M. Pardo,et al.  Smale’s 17th problem: Average polynomial time to compute affine and projective solutions , 2008 .

[12]  A. Edelman,et al.  Random matrix theory , 2005, Acta Numerica.

[13]  Diego Armentano,et al.  Complexity of Path-Following Methods for the Eigenvalue Problem , 2011, Found. Comput. Math..

[14]  Felipe Cucker,et al.  Fast computation of zeros of polynomial systems with bounded degree under finite-precision , 2014, Math. Comput..

[15]  Tien-Yien Li,et al.  Homotopy method for generalized eigenvalue problems Ax= ΛBx☆ , 1987 .

[16]  Felipe Cucker,et al.  On a problem posed by Steve Smale , 2009, 0909.2114.

[17]  Xingjie Helen Li,et al.  Numerical Solution of Dyson Brownian Motion and a Sampling Scheme for Invariant Matrix Ensembles , 2013, 1306.1179.

[18]  Beresford N. Parlett,et al.  The Symmetric Eigenvalue Problem (Classics in Applied Mathematics, Number 20) , 1999 .

[19]  L. Devroye Non-Uniform Random Variate Generation , 1986 .

[20]  Wolfgang Hörmann,et al.  Automatic Nonuniform Random Variate Generation , 2011 .