Fast Toeplitz eigenvalue computations, joining interpolation-extrapolation matrix-less algorithms and simple-loop conjectures: the preconditioned setting

03/21/2022
by   Manuel Bogoya, et al.
0

Under appropriate technical assumptions, the simple-loop theory allows to deduce various types of asymptotic expansions for the eigenvalues of Toeplitz matrices T_n(f) generated by a function f, unfortunately, such a theory is not available in the preconditioning setting, that is for matrices of the form T_n^-1(g)T_n(l) with l,g real-valued, g nonnnegative and not identically zero almost everywhere. Independently and under the milder hypothesis that f=l/g is even and monotonic over [0,π], matrix-less algorithms have been developed for the fast eigenvalue computation of large preconditioned matrices of the type above, within a linear complexity in the matrix order: behind the high efficiency of such algorithms there are the expansions as in the case g≡ 1, combined with the extrapolation idea, and hence we conjecture that the simple-loop theory has to be extended in such a new setting, as the numerics strongly suggest.Here we focus our attention on a change of variable, followed by the asymptotic expansion of the new variable, and we consider new matrix-less algorithms ad hoc for the current case. Numerical experiments show a much higher precision till machine precision and the same linear computation cost, when compared with the matrix-less procedures already proposed in the literature.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
01/06/2022

Fast Toeplitz eigenvalue computations, joining interpolation-extrapolation matrix-less algorithms and simple-loop theory

Under appropriate technical assumptions, the simple-loop theory allows t...
research
11/28/2022

Fast non-Hermitian Toeplitz eigenvalue computations, joining matrix-less algorithms and FDE approximation matrices

The present work is devoted to the eigenvalue asymptotic expansion of th...
research
11/07/2022

Asymptotics of the Sketched Pseudoinverse

We take a random matrix theory approach to random sketching and show an ...
research
08/05/2020

The Eight Epochs of Math as regards past and future Matrix Computation

This paper gives a personal assessment of Epoch making advances in Matri...
research
05/24/2023

On the eigenvalues of Toeplitz matrices with two off-diagonals

Consider the Toeplitz matrix T_n(f) generated by the symbol f(θ)=f̂_r e^...
research
03/14/2022

A note on the eigenvalues, singular values, and eigenvectors of Toeplitz and Hankel matrices

In a series of recent papers the spectral behavior of the matrix sequenc...

Please sign up or login with your details

Forgot password? Click here to reset