The block decomposition of a vandermonde matrix and its applications

Springer Science and Business Media LLC - Tập 21 - Trang 505-517 - 1981
W. P. Tang1, G. H. Golub
1Computer Science Department, Stanford University, Stanford, USA

Tóm tắt

This paper presents a new algorithm for the solution of linear equations with a Vandermonde coefficient matrix. The algorithm can also be used to solve the dual problem. Since the algorithm uses a block decomposition of the matrix, it is especially suitable for parallel computation. A variation of the block decomposition leads to the efficient solution of the interpolation problem with complex-conjugate interpolation points where the coefficients of the interpolating polynomial are real. In addition the algorithm can be used to solve some kinds of confluent Vandermonde systems.

Tài liệu tham khảo

Å. Björck and V. Pereyra,Solution of Vandermonde systems of equations, Math. Comp. 24 (1970), 893–904. Å. Björck and T. Elfving,Algorithms for confluent Vandermonde systems of equations, Numer. Math. 21 (1973), 130–137. G. Galimberti and V. Pereyra,Solving confluent Vandermonde systems of Hermite type, Numer. Math. 18 (1971), 44–60. W. Gautschi,On the inverse of Vandermonde and confluent Vandermonde matrix I, II, Numer. Math., v.4, 1962, pp. 117–123; ibid., v. 5, 1963, pp. 425–430. S. Å. Gustafson,Rapid computation of general interpolation formulae and mechanical quadrature rules, Comm. ACM, 14 (1971), 797–801. R. D. Skeel,Iterative refinement implies numerical stability for Gaussian elimination, Math. Comp. 34 (1980), 817–832.