×

A weakly stable algorithm for Padé approximants and the inversion of Hankel matrices. (English) Zbl 0779.65009

Authors’ summary: A new algorithm, algorithm NPADE, is presented for numerically computing Padé approximants in a weakly stable fashion. By this, it is meant that if the problem is well conditioned, then algorithm NPADE produces a good solution. No restrictions are imposed on the problem being solved. Except in certain pathological cases, the cost of the algorithm is \(O(n^ 2)\), where \(n\) is the maximum degree of the polynomials comprising the Padé approximant.
The operation of algorithm NPADE is controlled by a single parameter. Bounds are obtained for the computed solution and it is seen that they are a function of this parameter. Experimental results show that the bounds, while crude, reflect the actual behavior of the error. In addition, it is shown how better bounds can easily be obtained a posteriori. As another application of algorithm NPADE, it is shown that it can be used to compute stably, in a weak sense, the inverse of a Hankel or Toeplitz matrix.
Reviewer: D.Braess (Bochum)

MSC:

65D15 Algorithms for approximation of functions
65F05 Direct numerical methods for linear systems and matrix inversion
41A21 Padé approximation
PDFBibTeX XMLCite
Full Text: DOI