A Stable Algorithm for Multi-dimensional Padé Systems and the Inversion of Generalized Sylvester Matrices

• Author(s) / Creator(s)

  • Jones, Anthony
  • Cabay, Stan
  • Labahn, George

• Technical report TR94-07. For k+1 power series a0(z), ..., ak(z), we present a new iterative, look-ahead algorithm for numerically computing Pade'-Hermite systems and simultaneous Pade' systems along a diagonal of the associated Pade' tables. The algorithm computes the systems at all those points along the diagonal at which the associated striped Sylvester and mosaic Sylvester matrices are well-conditioned. It is shown that a good estimate for the condition numbers of these Sylvester matrices at a point is easily determined from the Pade'-Hermite system and simultaneous Pad'e system computed at that point. The operation and the stability of the algorithm is controlled by a single parameter tau which serves as a threshold in deciding if the Sylvester matrices at a point are sufficiently well-conditioned. We show that the algorithm is weakly stable, and provide bounds for the error in the computed solutions as a function of tau. Experimental results are given which show that the bounds reflect the actual behavior of the error. The algorithm requires O(||n||^2 + s^2 ||n||) operations, to compute Pade'-Hermite and simultaneous Pade' systems of type n=[n0, ..., nk], where ||n|| = n0 + ... + nk and s is the largest step-size taken along the diagonal. An additional application of the algorithm is the stable inversion of striped and mosaic Sylvester matrices. | TRID-ID TR94-07

• Date created

  1994

• Subjects / Keywords

  • Mosiac Sylvester inverses
  • Striped Sylvester inverses
  • Simultaneous Pade' approximants
  • Pade'-Hermite approximants
  • Numerical algorithm
  • Striped Hankel inverses
  • Numerical stability

• Type of Item

  Report

• DOI

  https://doi.org/10.7939/R3NG4GV7G

• License

  Attribution 3.0 International