Abstract:We consider the quaternion linear system for the unknown matrix X , where A , B are given , matrices with quaternion entries, motivated by applications that arise from fields such as quantum mechanics and signal processing. Our primary concern is the large-scale setting when n is large so that direct solutions are not feasible. We describe a block Krylov subspace method for the iterative solution of these quaternion linear systems. One difference compared to usual block Krylov subspace methods over complex Euclidean spaces is that the multiplication of quaternion scalars is not commutative. We describe a block quaternion Arnoldi process, taking noncommutativity features of quaternions into account, to generate an orthonormal basis for the quaternion Krylov space , where and is an initial guess for the solution. Then the best solution of in the least-squares sense is sought in the generated Krylov space. We explain how these least-squares problems over quaternion Krylov spaces can be solved efficiently by means of Householder reflectors. Most notably, we analyze rigorously the convergence of the proposed block quaternion GMRES approach when A is diagonalizable, and in the more general setting when A is not necessarily diagonalizable by making use of the Jordan form of A . Finally, we report numerical results that confirm the validity of the deduced theoretical convergence results, in particular illustrate that the proposed block quaternion Krylov subspace method converges quickly when A has clustered eigenvalues.

A block quaternion GMRES method and its convergence analysis

Gl-QFOM and Gl-QGMRES: two efficient algorithms for quaternion linear systems with multiple right-hand sides

Generalized Block Lanczos Methods for Large Unsymmetric Eigenproblems

A spectrally preconditioned and initially deflated variant of the restarted block GMRES method for solving multiple right-hand sides linear systems

On the backward stability of s-step GMRES

Convergence of block iterative methods for linear systems with generalized H-matrices

Generalized Jacobi Method for Computing Eigenvalues of Dual Quaternion Hermitian Matrices

A Variation on the Block Arnoldi Method for Large Unsymmetric Matrix Eigenproblems

The Convergence of Krylov Subspace Methods for Large Unsymmetric Linear Systems

Block GMRES Method with Inexact Breakdowns and Deflated Restarting

QQMR: A Structure-Preserving Quaternion Quasi-Minimal Residual Method for Non-Hermitian Quaternion Linear Systems

A block Lanczos method for large-scale quadratic minimization problems with orthogonality constraints

Structure Preserving Quaternion Biconjugate Gradient Method

The Convergence of Generalized Lanczos Methods for Large Unsymmetric Eigenproblems

Global FOM and GMRES algorithms for a class of complex matrix equations.

Convergence analysis of inner-iteration preconditioned GMRES

A quantile-based block Kaczmarz algorithm for solving large consistent linear systems

On block Gaussian sketching for the Kaczmarz method

On IGMRES: an Incomplete Generalized Minimal Residual Method for Large Unsymmetric Linear Systems

On Rayleigh Quotient Iteration for Dual Quaternion Hermitian Eigenvalue Problem

Randomized Large-Scale Quaternion Matrix Approximation: Practical Rangefinders and One-Pass Algorithm