A stabilized superfast solver for indefinite Hankel systems
This paper deals with the indefinite Hankel systems whose size is a power of 2. Based on an inversion formula for Loewner matrices and a divide and conquer approach, a stabilized super fast solver is derived. The algorithm is stabilized via pivoting, iterative improvement, and by giving the so-called ``difficult'' interpolation points an adequate treatment. Numerical examples are presented to illustrate the effectiveness of this algorithm.
scientific article; zbMATH DE number 1371019
scientific article