LSQR algorithm with structured preconditioner for the least squares problem in quaternionic quantum theory

Abstract The solution of a linear quaternionic least squares (QLS) problem can be transformed into that of a linear least squares (LS) problem with JRS-symmetric real coefficient matrix, which is suitable to be solved by developing structured iterative methods when the coefficient matrix is large and sparse. The main aim of this work is to construct a structured preconditioner to accelerate the LSQR convergence. The preconditioner is based on structure-preserving tridiagonalization to the real counterpart of the coefficient matrix of the normal equation, and the incomplete inverse upper–lower factorization related to only one symmetric positive definite tridiagonal matrix rather than four, so it is reliable and has low storage requirements. The performances of the LSQR algorithm with structured preconditioner are demonstrated by numerical experiments.

[1]  Xiaoye S. Li,et al.  Solution of a three-body problem in quantum mechanics using sparse linear algebra on parallel computers , 2001, SC.

[2]  David Finkelstein,et al.  Quaternionic Representations of Compact Groups , 1963 .

[3]  Davod Khojasteh Salkuyeh,et al.  A PRECONDITIONER FOR THE LSQR ALGORITHM , 2008 .

[4]  J. Navarro-Pedreño.  Numerical Methods for Least Squares Problems , 1996 .

[5]  M. Benzi.  Preconditioning techniques for large linear systems: a survey , 2002 .

[6]  Timothy A. Davis,et al.  The university of Florida sparse matrix collection , 2011, TOMS.

[7]  Angelika Bunse-Gerstner,et al.  An improved preconditioned LSQR for discrete ill-posed problems , 2006, Math. Comput. Simul..

[8]  Musheng Wei,et al.  Equality constrained least squares problem over quaternion field , 2003, Appl. Math. Lett..

[9]  Stefano De Leo,et al.  Quaternionic Electron Theory: Dirac's Equation , 1998, hep-th/9806057.

[10]  Stefano De Leo.  Quaternions and special relativity , 1996 .

[11]  Yong Duan,et al.  A comparative study of iterative solutions to linear systems arising in quantum mechanics , 2010, J. Comput. Phys..

[12]  Minghui Wang.  Algorithm Q-LSQR for the least squares problem in quaternionic quantum theory , 2010, Comput. Phys. Commun..

[13]  M. Benzi,et al.  Block preconditioning of real-valued iterative algorithms for complex linear systems , 2007 .

[14]  P. B. Visscher,et al.  Quaternion-based algorithm for micromagnetics , 2002 .

[15]  Yan Feng,et al.  An iterative algorithm for least squares problem in quaternionic quantum theory , 2008, Comput. Phys. Commun..

[16]  Valeria Simoncini,et al.  Recent computational developments in Krylov subspace methods for linear systems , 2007, Numer. Linear Algebra Appl..

[17]  Jianli Zhao,et al.  A new technique of quaternion equality constrained least squares problem , 2008 .

[18]  C. D. Lin,et al.  Linear-least-squares fitting method for the solution of the time-dependent Schrodinger equation: Applications to atoms in intense laser fields , 2000 .

[19]  Å. Björck,et al.  Stability of Conjugate Gradient and Lanczos Methods for Linear Least Squares Problems , 1998, SIAM J. Matrix Anal. Appl..

[20]  Davies,et al.  Observability of quaternionic quantum mechanics. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[21]  Musheng Wei,et al.  A new structure-preserving method for quaternion Hermitian eigenvalue problems , 2013, J. Comput. Appl. Math..

[22]  L. Kolotilina,et al.  Factorized Sparse Approximate Inverse Preconditionings I. Theory , 1993, SIAM J. Matrix Anal. Appl..

[23]  Nicolas Le Bihan,et al.  Singular value decomposition of quaternion matrices: a new tool for vector-sensor signal processing , 2004, Signal Process..

[24]  Ling Xia,et al.  The application of subspace preconditioned LSQR algorithm for solving the electrocardiography inverse problem. , 2009, Medical engineering & physics.

[25]  P. Hansen,et al.  Subspace Preconditioned LSQR for Discrete Ill-Posed Problems , 2003 .

[26]  Tongsong Jiang,et al.  Two Novel Algebraic Techniques for Quaternion Least Squares Problems in Quaternionic Quantum Mechanics , 2016 .

[27]  Michele Benzi,et al.  A Robust Preconditioner with Low Memory Requirements for Large Sparse Least Squares Problems , 2003, SIAM J. Sci. Comput..

[28]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[29]  Guillaume Alléon,et al.  Sparse approximate inverse preconditioning for dense linear systems arising in computational electromagnetics , 1997, Numerical Algorithms.

[30]  Minghui Wang,et al.  Hermitian tridiagonal solution with the least norm to quaternionic least squares problem , 2010, Comput. Phys. Commun..

[31]  Zhong-zhi,et al.  A SHIFT-SPLITTING PRECONDITIONER FOR NON-HERMITIAN POSITIVE DEFINITE MATRICES , 2006 .

[32]  S. Adler,et al.  Quaternionic quantum mechanics and quantum fields , 1995 .

[33]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[34]  G. Golub,et al.  Block Preconditioning for the Conjugate Gradient Method , 1985 .

[35]  L. Xia,et al.  Combination of the LSQR method and a genetic algorithm for solving the electrocardiography inverse problem. , 2007, Physics in medicine and biology.

[36]  Yanfeng Sun,et al.  Color face recognition based on quaternion matrix representation , 2011, Pattern Recognit. Lett..

[37]  Samuel Schiminovich,et al.  Foundations of Quaternion Quantum Mechanics , 1962 .

[38]  Li Chen,et al.  Algebraic algorithms for least squares problem in quaternionic quantum theory , 2007, Comput. Phys. Commun..

[39]  Li Chen,et al.  An algebraic method for Schrödinger equations in quaternionic quantum mechanics , 2008, Comput. Phys. Commun..

[40]  Tongsong Jiang,et al.  An algebraic technique for total least squares problem in quaternionic quantum theory , 2016, Appl. Math. Lett..

[41]  A. Yu. Yeremin,et al.  Factorized Sparse Approximate Inverse Preconditioning II: Solution of 3D FE Systems on Massively Parallel Computers , 1995, Int. J. High Speed Comput..

[42]  Adler.  Scattering and decay theory for quaternionic quantum mechanics, and the structure of induced T nonconservation. , 1988, Physical review. D, Particles and fields.

[43]  Tongsong Jiang,et al.  An algebraic method for quaternion and complex Least Squares coneigen-problem in quantum mechanics , 2014, Appl. Math. Comput..

[44]  Yunzhong Shen,et al.  A Quaternion-based Geodetic Datum Transformation Algorithm , 2006 .

[45]  Jun-Feng Yin,et al.  Preconditioned GMRES methods with incomplete Givens orthogonalization method for large sparse least-squares problems , 2009 .