Abstract
A quasi-Toeplitz M-matrix A is an infinite M-matrix that can be written as the sum of a semi-infinite Toeplitz matrix and a correction matrix. This paper is concerned with computing the square root of invertible quasi-Toeplitz M-matrices which preserves the quasi-Toeplitz structure. We show that the Toeplitz part of the square root can be easily computed through evaluation/interpolation. This advantage allows us to propose algorithms solely for the computation of correction part, whence we propose a fixed-point iteration and a structure-preserving doubling algorithm. Additionally, we show that the correction part can be approximated by solving a nonlinear matrix equation with coefficients of finite size followed by extending the solution to infinity. Numerical experiments showing the efficiency of the proposed algorithms are performed.



Similar content being viewed by others
Data Availability
All data used in the manuscript is numerically generated using MATLAB. The MATLAB source code used to generate the numerical results is available at https://github.com/JieMeng00/structured_sqrtm_square_root_m-matrices.
References
Alefeld, G., Schneider, N.: On square roots of \(M\)-matrices. Linear Algebra Appl. 42, 119–132 (1982)
B\(\ddot{\rm o}\)ttcher, A. and Grudsky, S.M.: Spectral Properties of Banded Toeplitz Matrices. SIAM, Philadelphia, PA (2005)
Bini, D.A., Iannazzo, B., and Meng, J.: Algorithms for approximating means of semi-definite quasi-Toeplitz matrices, in: International Conference on Geometric Science of Information, GSI 2021: Geometric Science of Information, 2021, pp. 405–414
Bini, D.A., Iannazzo, B. and Meini, B.: Numerical Solution of Algebraic Riccati Equations. Volume 9 of Fundamentals of Algorithms. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA (2012)
Bini, D.A. and Meini, B.: A defect-correction algorithm for quadratic matrix equations, with applications to quasi-Toeplitz matrices. arXiv preprint (2022)
Bini, D.A., Massei, S., Meini, B.: On functions of quasi Toeplitz matrices. Sb. Math. 208, 56–74 (2017)
Bini, D.A., Massei, S., Meini, B.: Semi-infinite quasi-Toeplitz matrices with applications to QBD stochastic processes. Math. Comp. 87, 2811–2830 (2018)
Bini, D.A., Massei, S., Meini, B., Robol, L.: On quadratic matrix equations with infinite size coefficients encountered in QBD stochastic processes. Numer. Linear Algebra Appl. 25, e2128 (2018)
Bini, D.A., Massei, S., Robol, L.: Quasi-Toeplitz matrix arithmetic: a MATLAB toolbox. Numer. Algorithms 81, 741–769 (2019)
Bini, D.A., Meini, B., Meng, J.: Solving quadratic matrix equations arising in random walks in the quarter plane. SIAM J. Matrix Anal. Appl. 41, 691–714 (2020)
Bini, D.A., Massei, S., Meini, B., Robol, L.: A computational framework for two-dimensional random walks with restarts. SIAM J. Sci. Comput. 42(4), A2108–A2133 (2020)
Bini, D.A., Iannazzo, B., Meng, J.: Geometric mean of quasi-Toeplitz matrices. BIT 63, 20 (2023)
Bini, D.A., Iannazzo, B., Meini, B., Meng, J., Robol, L.: Computing eigenvalues of semi-infinite quasi-Toeplitz matrices. Numer. Algorithms 92, 89–118 (2023)
Higham, N.J.: Functions of Matrices: Theory and Computation. Society for Industrial and Applied Mathematics, Philadelphia, PA, USA (2008)
Huang, T.-M., Li, R.-C. and Lin, W.-W.: Structure-preserving Doubling Algorithms for Nonlinear Matrix Equations. Volume 14 of Fundamentals of Algorithms. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA (2018)
Kadison, R.V. and Ringrose, J.R.: Fundamentals of the Theory of Operator Algebras. Vol. I, volume 100 of Pure and Applied Mathematics. Academic Press, Inc. [Harcourt Brace Jovanovich, Publishers], New York, 1983. Elementary theory
Kannan, M.R., Sivakumar, K.C.: On certain positivity classes of operators. Numer. Funct. Anal. Optim. 37, 206–224 (2017)
Kim, H.-M., Meng, J.: Structured perturbation analysis for an infinite size quasi-Toeplitz matrix equation with applications. BIT. 61, 859–879 (2021)
Marek, I.: Frobenius theory of positive operators: comparison theorems and applications. SIAM J. Appl. Math. 19, 607–628 (1970)
Marek, I.: On square roots of M-operators. Linear Algebra Appl. 223(224), 501–520 (1995)
Marek, I., Szyld, D.B.: Splittings of \(M\)-operators: irreducibility and the index of the iteration operator. Numer. Funct. Anal. Optim. 11, 529–553 (1990)
Meng, J.: Theoretical and computational properties of semi-infinite quasi-Toeplitz \(M\)-matrices. Linear Algebra Appl. 653, 66–85 (2022)
Motyer, A.J., Taylor, P.G.: Decay rates for quasi-birth-and-death processes with countably many phases and tridiagonal block generators. Adv. Appl. Prob. 38, 522–544 (2006)
Robol, L.: Rational Krylov and ADI iteration for infinite size quasi-Toeplitz matrix equations. Linear Algebra Appl. 604, 210–235 (2020)
Shivakumar, P.N., Sivakumar, K.C. and Zhang,Y.: Infinite Matrices and Their Recent Applications. Springer International Publishing Switzerland (2016)
Acknowledgements
The authors wish to thank the anonymous referees for providing useful comments that helped to improve this presentation.
Funding
The work of the first author was partly supported by the National Natural Science Foundation of China under grant Nos. 12001262, 12371379, 12361080 and by Jiangxi Provincial Natural Science Foundation under grant No. 20224BAB211006. The work of the third author was partly supported by the National Natural Science Foundation of China under grant No. 12201591 and Qingdao Natural Science Foundation under project 23-2-1-158-zyyd-jch.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Chen, H., Kim, HM. & Meng, J. Algorithms for Square Root of Semi-Infinite Quasi-Toeplitz M-Matrices. J Sci Comput 99, 66 (2024). https://doi.org/10.1007/s10915-024-02491-8
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-024-02491-8