Institution: | a Dipartimento di Matematica e Informatica, Università di Udine, I-33100 Udine, Italy b Dipartimento di Matematica, Università degli Studi di Pisa, Via Buonarroti 2, I-56127 Pisa, Italy |
Abstract: | A classical result of structured numerical linear algebra states that the inverse of a nonsingular semiseparable matrix is a tridiagonal matrix. Such a property of a semiseparable matrix has been proved to be useful for devising linear complexity solvers, for establishing recurrence relations among its columns or rows and, moreover, for efficiently evaluating its characteristic polynomial. In this paper, we provide sparse structured representations of a semiseparable matrix A which hold independently of the fact that A is singular or not. These relations are found by pointing out the band structure of the inverse of the sum of A plus a certain sparse perturbation of minimal rank. Further, they can be used to determine in a computationally efficient way both a reflexive generalized inverse of A and its characteristic polynomial. |