Weyl distributions, spectral properties, and circulant approximation results for quaternion block multilevel Toeplitz matrix sequences
Ayoub Lailoune, Valerio Loi, Stefano Serra-Capizzano
TL;DR
The paper extends Weyl-type eigenvalue and singular value distribution theory to single-axis quaternion block multilevel Toeplitz sequences generated by quaternion-valued symbols. It builds a robust quaternion a.c.s. framework and shows that circulant matrix sequences provide effective approximations for preconditioning, all anchored by a symplectic embedding that transfers quaternionic problems to the complex GLT setting. Central contributions include precise canonical eigenvalue formulations, Schatten-norm estimates, and distribution results for both eigenvalues and singular values under left/right and sandwich kernels, along with practical circulant-a.c.s. constructions. Numerical experiments corroborate the theory across continuous and $L^{1}$ symbols, Hermitian and non-Hermitian cases, highlighting potential for fast, structure-preserving solvers in large-scale quaternion linear systems.
Abstract
The present work contains a comprehensive treatment of Weyl eigenvalue and singular value distributions for single-axis quaternion block multilevel Toeplitz matrix sequences generated by $s\times t$ quaternion matrix-valued, $d$-variate, Lebesgue integrable generating functions. Furthermore, in view of concrete applications, we are interested in preconditioning and matrix approximation results. To this end, a crucial step is the extension of the notion of an approximating class of sequences (a.c.s.) to the case of matrix sequences with quaternion entries, since it allows us to decompose the difference between a matrix and its preconditioner into low-norm plus (relatively) low-rank terms. As a specific example, we consider classes of quaternion block multilevel circulant matrix sequences as an a.c.s. for quaternion block multilevel Toeplitz matrix sequences. These approximation results lay the foundations for fast preconditioning methods when dealing with large quaternion linear systems stemming from modern applications. We conclude our study with numerical experiments and directions for future research.
