A note on the O(n)-storage implementation of the GKO algorithm and its adaptation to Trummer-like matrices

We propose a new O(n)-space implementation of the GKO-Cauchy algorithm for the solution of linear systems where the coefficient matrix is Cauchy-like. Moreover, this new algorithm makes a more efficient use of the processor cache memory; for matrices of size larger than n ≈ 500 − 1000, it outperforms the customary GKO algorithm. We present an applicative case of Cauchy-like matrices with nonreconstructible main diagonal. In this special instance, the O(n) space algorithms can be adapted nicely to provide an efficient implementation of basic linear algebra operations in terms of the low displacement-rank generators.

• Type: Journal
• Year: 2010
• Where: NA
• Authors: Federico Poloni
• Keywords: Algorithms, Customary Gko Algorithm, NA 2010, Processor Cache Memory