Two Iterative algorithms for the matrix sign function based on the adaptive filtering technology

In this paper, two new efficient algorithms for calculating the sign function of the large-scale sparse matrix are proposed by combining filtering algorithm with Newton method and Newton Schultz method respectively. Through the theoretical analysis of the error diffusion in the iterative process, we designed an adaptive filtering threshold, which can ensure that the filtering has little impact on the iterative process and the calculation result. Numerical experiments are consistent with our theoretical analysis, which shows that the computational efficiency of our method is much better than that of Newton method and Newton Schultz method, and the computational error is of the same order of magnitude as that of the two methods.

[1]  Feng Wu,et al.  High-performance computation of the exponential of a large sparse matrix , 2021, SIAM J. Matrix Anal. Appl..

[2]  C. Chun,et al.  Construction of stable and globally convergent schemes for the matrix sign function , 2019, Linear Algebra and its Applications.

[3]  Ashim Kumar,et al.  A fourth-order method for computing the sign function of a matrix with application in the Yang–Baxter-like matrix equation , 2019, Comput. Appl. Math..

[4]  H. B. Jebreen.  Constructing a High-Order Globally Convergent Iterative Method for Calculating the Matrix Sign Function , 2018 .

[5]  Fazlollah Soleymani,et al.  A Novel Iterative Method for Polar Decomposition and Matrix Sign Function , 2015 .

[6]  J. A. Ezquerro,et al.  Iterative methods for computing the matrix square root , 2015 .

[7]  Fazlollah Soleymani,et al.  A fast convergent numerical method for matrix sign function with application in SDEs , 2015, J. Comput. Appl. Math..

[8]  F. Khaksar Haghani,et al.  A generalized Steffensen's method for matrix sign function , 2015, Appl. Math. Comput..

[9]  F. W. Williams,et al.  A fast precise integration method for structural dynamics problems , 2012 .

[10]  Nicholas J. Higham,et al.  Functions of matrices - theory and computation , 2008 .

[11]  M. Sh. Misrikhanov,et al.  Matrix sign function in the problems of analysis and design of the linear systems , 2008 .

[12]  Enrique S. Quintana-Ortí,et al.  Efficient algorithms for generalized algebraic Bernoulli equations based on the matrix sign function , 2007, Numerical Algorithms.

[13]  Enrique S. Quintana-Ortí,et al.  Parallel solution of Riccati matrix equations with the matrix sign function , 1998, Autom..

[14]  C. Pearce,et al.  On the matrix-sign-function method for solving algebraic Riccati equations , 1997 .

[15]  JUDITH D. GARDINER,et al.  A Stabilized Matrix Sign Function Algorithm for Solving Algebraic Riccati Equations , 1997, SIAM J. Sci. Comput..

[16]  J. D. Roberts,et al.  Linear model reduction and solution of the algebraic Riccati equation by use of the sign function , 1980 .

[17]  高强,et al.  Fast precise integration method for hyperbolic heat conduction problems , 2013 .

[18]  M. Benzi,et al.  DECAY BOUNDS AND ( ) ALGORITHMS FOR APPROXIMATING FUNCTIONS OF SPARSE MATRICES , 2007 .

[19]  E. S. Quintana‐Ortí,et al.  Solving stable generalized Lyapunov equations with the matrix sign function , 1999 .

[20]  Ronald R. Coifman,et al.  Wavelet-Like Bases for the Fast Solution of Second-Kind Integral Equations , 1993, SIAM J. Sci. Comput..

[21]  R. Byers.  Solving the algebraic Riccati equation with the matrix sign function , 1987 .