ERA

Download the full-sized PDF of A Stable Algorithm for Multi-dimensional Padé Systems and the Inversion of Generalized Sylvester MatricesDownload the full-sized PDF

Analytics

Share

Permanent link (DOI): https://doi.org/10.7939/R3NG4GV7G

Download

Export to: EndNote  |  Zotero  |  Mendeley

Communities

This file is in the following communities:

Computing Science, Department of

Collections

This file is in the following collections:

Technical Reports (Computing Science)

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

Descriptions

Author or creator
Cabay, Stan
Jones, Anthony
Labahn, George
Additional contributors
Subject/Keyword
striped Hankel inverses
Pade'-Hermite approximants
simultaneous Pade' approximants
mosiac Sylvester inverses
numerical stability
numerical algorithm
striped Sylvester inverses
Type of item
Computing Science Technical Report
Computing science technical report ID
TR94-07
Language
English
Place
Time
Description
Technical report TR94-07. For k+1 power series a_0(z), ..., a_k(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=[n_0, ..., n_k], where ||n|| = n_0 + ... + n_k 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.
Date created
1994
DOI
doi:10.7939/R3NG4GV7G
License information
Creative Commons Attribution 3.0 Unported
Rights

Citation for previous publication

Source
Link to related item

File Details

Date Uploaded
Date Modified
2014-04-29T19:23:30.597+00:00
Audit Status
Audits have not yet been run on this file.
Characterization
File format: pdf (Portable Document Format)
Mime type: application/pdf
File size: 1241008
Last modified: 2015:10:12 17:29:46-06:00
Filename: TR94-07.pdf
Original checksum: 866b0c78ffed00c145720dd1db396adb
Well formed: true
Valid: true
Page count: 55
Activity of users you follow
User Activity Date