Improved Padé approximation using SVD

"Jake Bowhay" <[email protected]> Mon, 04 Mar 2024 11:34:08 -0000
Newsgroups gmane.comp.python.scientific.devel
Message-ID <[email protected]>
SGkgYWxsLA0KDQpJIHdhbnQgdG8gcHJvcG9zZSBpbXByb3ZlbWVudHMgdG8gdGhlIGN1cnJlbnQg
aW1wbGVtZW50YXRpb24gb2YgUGFkw6kgYXBwcm94aW1hbnRzIG9mZmVyZWQgd2l0aGluIFNjaVB5
LiBDdXJyZW50bHksIHRoZSBtZXRob2Qgd2Ugb2ZmZXIgc29sdmVzIGEgdG9lcGxpdHogc3lzdGVt
IHdoaWNoIGlzIG9mdGVuIGhpZ2hseSBpbGwtY29uZGl0aW9uZWQuIFRoaXMgbGVhZHMgdG8gc3B1
cmlvdXMgcG9sZS16ZXJvIHBhaXJzIG9yIEZyb2lzc2FydCBkb3VibGV0cyB3aGljaCBtZWFucyBp
biBwcmFjdGljZSBvZnRlbiB0aGUgaW1wbGVtZW50YXRpb24gb2ZmZXJlZCBpbiBTY2lQeSBpcyB1
bnVzYWJsZS4gRm9yIGV4YW1wbGUsIEkgYW0gaW50ZXJlc3RlZCBpbiBQYWTDqSBhcHByb3hpbWF0
aW9ucyB0byB0aGUgc29sdXRpb24gb2YgdGhlIFZvbHRlcnJhIGludGVncm8tZGlmZmVyZW50aWFs
IGVxdWF0aW9ucy4gSGVyZSB0aGUgbWV0aG9kIGluIFNjaVB5IHF1aWNrbHkgYnJlYWtzIGRvd24u
DQoNCkkgd291bGQgbGlrZSB0byBwcm9wb3NlIHdlIHJlcGxhY2Ugb3VyIGN1cnJlbnQgbWV0aG9k
IHdpdGggYSBtZXRob2QgYmFzZWQgb24gdGhlIHBhcGVyIGJ5IFRyZWZldGhlbiBldCBhbC4gWzFd
LiBUaGlzIHBhcGVyIGhhcyBxdWl0ZSBhIGxvdCBvZiBjaXRhdGlvbnMgYW5kIGlzIGFscmVhZHkg
aW1wbGVtZW50ZWQgaW4gdGhlIHBvcHVsYXIgQ2hlYmZ1biBwYWNrYWdlLiBIb3BlZnVsbHksIHRo
aXMgcGFydCBpcyBmYWlybHkgdW5jb250cm92ZXJzaWFsIGFzIGl0IGlzIGEgY2xlYXIgaW1wcm92
ZW1lbnQgdG8gc29tZXRoaW5nIHdlIGFscmVhZHkgb2ZmZXIuIEhvd2V2ZXIsIHRoZXJlIGFyZSBh
IGNvdXBsZSBvZiBxdWVzdGlvbnMgYWJvdXQgYWRkaW5nIHRoaXM6DQoNCjEuIExvY2F0aW9uLiAg
Q3VycmVudGx5LCBvdXIgUGFkw6kgYXBwcm94aW1hdGlvbiBsaXZlcyB3aXRoaW4gc2NpcHkuaW50
ZXJwb2xhdGUuIENsZWFybHksIHRoaXMgaXMgbm90IGFuIGludGVycG9sYXRpb24gbWV0aG9kIHNv
IHRoaXMgd291bGQgYmUgYSBnb29kIG9wcG9ydHVuaXR5IHRvIHRoaW5rIGFib3V0IGlmIHRoZXJl
IGlzIGEgYmV0dGVyIHBsYWNlIGZvciB0aGlzIHRvIGxpdmUuDQoNCjIuIEltcGxlbWVudCB3aXRo
aW4gdGhlIG9sZCBzaWduYXR1cmUgb3Igc3RhcnQgYWdhaW4uIEl0IHdhcyBzdWdnZXN0ZWQgdGhh
dCBpdCB3b3VsZCBiZSBlYXNpZXIgdG8gc3RhcnQgYWZyZXNoIGFuZCBkZXByZWNhdGUgdGhlIG9s
ZCBmdW5jdGlvbi4gSSB0aGluayB0aGlzIHdvdWxkIGJlIHByZWZlcmFibGUgYXMgd2Ugd291bGRu
J3QgYmUgc2hhY2tsZWQgYnkgb2xkIGRlc2lnbiBkZWNpc2lvbnMuDQoNCkFzIGV2ZXIsIHBsZWFz
ZSBmZWVsIGZyZWUgdG8gam9pbiBpbiB0aGUgZGlzY3Vzc2lvbiBlaXRoZXIgaGVyZSBvciBvbiB0
aGUgb3JpZ2luYWwgaXNzdWUgWzJdLg0KDQpSZWdhcmRzLA0KSmFrZQ0KDQpbMV0gaHR0cHM6Ly9w
ZW9wbGUubWF0aHMub3guYWMudWsvdHJlZmV0aGVuL3B1YmxpY2F0aW9uL1BERi8yMDExXzE0NC5w
ZGYNClsyXSBodHRwczovL2dpdGh1Yi5jb20vc2NpcHkvc2NpcHkvaXNzdWVzLzIwMDY0DQpfX19f
X19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fXwpTY2lQeS1EZXYgbWFp
bGluZyBsaXN0IC0tIHNjaXB5LWRldkBweXRob24ub3JnClRvIHVuc3Vic2NyaWJlIHNlbmQgYW4g
ZW1haWwgdG8gc2NpcHktZGV2LWxlYXZlQHB5dGhvbi5vcmcKaHR0cHM6Ly9tYWlsLnB5dGhvbi5v
cmcvbWFpbG1hbjMvbGlzdHMvc2NpcHktZGV2LnB5dGhvbi5vcmcvCk1lbWJlciBhZGRyZXNzOiBn
Y3BzZC1zY2lweS1kZXZAbS5nbWFuZS1teC5vcmcK