TY - CHAP

T1 - Fast inversion of polynomial-vandermonde matrices for polynomial systems related to order one quasiseparable matrices

AU - Bella, T.

AU - Eidelman, Y.

AU - GohbergZ, I.

AU - Olshevsky, V.

AU - Tyrtyshnikov, E.

N1 - Publisher Copyright:
© 2013 Springer Basel.

PY - 2013

Y1 - 2013

N2 - While Gaussian elimination is well known to require O(n 3) operations to invert an arbitrary matrix, Vandermonde matrices may be inverted using O(n 2) operations by a method of Traub [24]. While this original version of the Traub algorithm was noticed to be unstable, it was shown in [12] that with a minor modification, the Traub algorithm can typically yield a very high accuracy. This approach has been extended from classical Vandermonde matrices to polynomial-Vandermonde matrices involving real orthogonal polynomials [3], [10], and Szegő polynomials [19]. In this paper we present an algorithm for inversion of a class of polynomial-Vandermonde matrices with special structure related to order one quasiseparable matrices, generalizing monomials, real orthogonal polynomials, and Szegő polynomials. We derive a fast O(n 2) inversion algorithm applicable in this general setting, and explore its reduction in the previous special cases. Some very preliminary numerical experiments are presented, demonstrating that, as observed by our colleagues in previous work, good forward accuracy is possible in some circumstances, which is consistent with previous work of this type.

AB - While Gaussian elimination is well known to require O(n 3) operations to invert an arbitrary matrix, Vandermonde matrices may be inverted using O(n 2) operations by a method of Traub [24]. While this original version of the Traub algorithm was noticed to be unstable, it was shown in [12] that with a minor modification, the Traub algorithm can typically yield a very high accuracy. This approach has been extended from classical Vandermonde matrices to polynomial-Vandermonde matrices involving real orthogonal polynomials [3], [10], and Szegő polynomials [19]. In this paper we present an algorithm for inversion of a class of polynomial-Vandermonde matrices with special structure related to order one quasiseparable matrices, generalizing monomials, real orthogonal polynomials, and Szegő polynomials. We derive a fast O(n 2) inversion algorithm applicable in this general setting, and explore its reduction in the previous special cases. Some very preliminary numerical experiments are presented, demonstrating that, as observed by our colleagues in previous work, good forward accuracy is possible in some circumstances, which is consistent with previous work of this type.

KW - Inversion of vandermonde matrices

KW - Polynomial vandermonde matrices

KW - Quasiseparable matrices

UR - http://www.scopus.com/inward/record.url?scp=84975686329&partnerID=8YFLogxK

U2 - 10.1007/978-3-0348-0639-8_8

DO - 10.1007/978-3-0348-0639-8_8

M3 - ???researchoutput.researchoutputtypes.contributiontobookanthology.chapter???

AN - SCOPUS:84975686329

T3 - Operator Theory: Advances and Applications

SP - 79

EP - 106

BT - Operator Theory

PB - Springer International Publishing

ER -