On the NPHSS-KPIK iteration method for low-rank complex Sylvester equations arising from time-periodic fractional diffusion equations

Based on the Hermitian and skew-Hermitian (HS) splitting for non-Hermitian matrices, a nonalternating preconditioned Hermitian and skew-Hermitian splitting-Krylov plus inverted Krylov subspace (NPHSS-KPIK) iteration method for solving a class of large and low-rank complex Sylvester equations arising from the two-dimensional timeperiodic fractional diffusion problem is established. The local convergence condition is proposed and the optimal parameter is given. Numerical experiments are used to show the efficiency of the NPHSS-KPIK iteration method for solving the Sylvester equations arising from the time-periodic fractional diffusion equations.

PDF

___

[1] Axelsson O, Bai ZZ, Qiu SX. A class of nested iteration schemes for linear systems with a coefficient matrix with a dominant positive definite symmetric part. Numer Algorithms 2004; 35: 351-372.
[2] Bai J, Feng XC. Fractional-order anisotropic diffusion for image denoising. IEEE Trans. Image Process 2007; 16: 24922502.
[3] Bai ZZ. On Hermitian and skew-Hermitian splitting iteration methods for continuous Sylvester equations. J Comput Math 2011; 29: 185-198.
[4] Bai ZZ, Benzi M, Chen F. Modified HSS iteration methods for a class of complex symmetric linear systems. Computing 2010; 87: 93-111.
[5] Bai ZZ, Golub GH, Ng MK. Hermitian and skew-Hermitian splitting methods for non-Hermitian positive definite linear systems. SIAM J Matrix Anal Appl 2003; 24: 603-626.
[6] Bai ZZ, Qiu SX. Splitting-MINRES methods for linear systems with the coefficient matrix with a dominant indefinite symmetric part. Math Numer Sinica 2002; 24: 113-128 (in Chinese).
[7] Bartels RH, Stewart GW. Algorithm 432: Solution of the matrix equation AX+XB=C. Comm ACM 1972; 15: 820-826.
[8] Benson DA, Wheatcraft SW, Meerschaert MM. Application of a fractional advection-dispersion equation. Water Res Research 2000; 36: 1403-1412.
[9] Golub G, Nash S, Loan CV. A Hessenberg-Schur method for the problem AX+XB=C. IEEE Trans Automat Control 1979; 24: 909-913.
[10] Gu CQ, Xue HY. A shift-splitting hierarchical identification method for solving Lyapunov matrix equations. Linear Algebra Appl 2009; 430: 1517-1530.
[11] Harman HH, Jones WH. Factor analysis by minimizing residuals (minres). Psychometrika 1966; 31: 351-368.
[12] Jbilou K. Low rank approximate solutions to large Sylvester matrix equations. Appl Math Comput 2006; 177: 365-376.
[13] Li CX, Wu SL. A single-step HSS method for non-Hermitian positive definite linear systems. Appl Math Lett 2015; 44: 26-29.
[14] Magin RL. Fractional Calculus in Bioengineering. New York, NY, USA: Begell House Publishers, 2006.
[15] Meerschaert MM, Tadjeran C. Finite difference approximations for two-sided space-fractional partial differential equations. Appl Numer Math 2006; 56: 80-90.
[16] Oldham KB, Spanier J. The Fractional Calculus. New York, NY, USA: Academic Press, 1974.
[17] Podlubny I, Chechkin A, Skovranek T, Chen YQ, Jara BMV. Matrix approach to discrete fractional calculus II: Partial fractional differential equations. J Comput Phys 2009; 228: 3137-3153.
[18] Raberto M, Scalas E, Mainardi F. Waiting-times and returns in high-frequency financial data: an empirical study. Phys A: Statistical Mechanics and its Applications 2002; 314: 749-755.
[19] Saad Y, Schultz MH. GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J Sci Statist Comput 1986; 7: 856-869.
[20] Shlesinger MF, West BJ, Klafter J. L´evy dynamics of enhanced diffusion: Application to turbulence. Phys Rev Lett 1987; 58: 1100.
[21] Simoncini V. A new iterative method for solving large-scale Lyapunov matrix equations. SIAM J Sci Comput 2007; 29: 1268-1288.
[22] Smith RA. Matrix equation XA+BX=C. SIAM J Appl Math 1968; 16: 198-201.
[23] Sokolov IM, Klafter J, Blumen A. Fractional kinetics. Phys Today 2002; 55: 48-54.
[24] Starke G, Niethammer W. SOR for AX-XB=C. Linear Algebra Appl 1991; 154: 355375.
[25] Wachspress EL. Iterative solution of the Lyapunov matrix equation. Appl Math Lett 1988; 1: 87-90.
[26] Wang X, Li WW, Mao LZ. On positive-definite and skew-Hermitian splitting iteration methods for continuous Sylvester equation AX+ XB= C. Comput Math Appl 2013; 66: 2352-2361.
[27] Zaslavsky GM, Stevens D, Weitzner H. Self-similar transport in incomplete chaos. Phys Rev E(3) 1993; 48: 1683.
[28] Zeng ML, Zhang GF. Incomplete circulant and skew-circulant splitting iteration method for time-dependent space fractional diffusion equations. Jpn J Ind Appl Math 2016; 33: 251-268.
[29] Zheng QQ, Ma CF. On normal and skew-Hermitian splitting iteration methods for large sparse continuous Sylvester equations. J Comput Appl Math 2014; 268: 145-154.