Quasiseparable Hessenberg reduction of real diagonal plus low rank matrices and applications

We present a novel algorithm to perform the Hessenberg reduction of an n×nn×n matrix A   of the form A=D+UV⁎A=D+UV⁎ where D is diagonal with real entries and U and V   are n×kn×k matrices with k≤nk≤n. The algorithm has a cost of O(n2k)O(n2k) arithmetic operations and is based on the quasiseparable matrix technology. Applications are shown to solving polynomial eigenvalue problems and some numerical experiments are reported in order to analyze the stability of the approach.