arrow
Return

Computing the matrix square root with the doubling algorithm

delete2025-09-03
delete0
PRE
AI
P
Pei-Chang Guo *
X
Xiaoxia Guo
Y
Yi Tang
DOI:10.1007/s10915-025-03040-7delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
The matrix quare root is a very important matrix computation problem. The structure-preserving doubling algorithm (SDA) is quadratically convergent. We demonstrate how the matrix square root is computed by SDA, where the Cayley transform is deployed and the original equation is rewritten in terms of the invariant subspace. Further, we introduce a positive parameter $$\alpha $$ when using the Cayley transform and get SDA- $$\alpha $$ for $$A^{1/2}$$ . If the appropriate parameter $$\alpha $$ is given, SDA- $$\alpha $$ for $$A^{1/2}$$ can outperform SDA for $$A^{1/2}$$ in some examples. The convergence properties of the two algorithms are given. The algorithms in this paper are not numerically sensitive when the original matrix is ill-conditioned and the algorithms work when the original matrix is a singular matrix but has a square root.
Keywords:
Structure-preserving doubling algorithm
Cayley transform
Matrix square root
Quadratic convergence

Journal

Journal of Scientific Computing cover
Journal of Scientific Computing
IF:
3.3
Papers:
723
Citations:
9.6K

Organization

S
School of Mathematical Sciences
Scholars:
583
Papers: 335
Citations: 1
Cited Papers

Cited Papers

err1997-01-01
err0
PREAI
errDario Andrea Bini; Beatrice Meini
errShare
errSave
Approximating the Logarithm of a Matrix to Specified Accuracy
err2001-01-01
err0
errOAAI
errSheung Hun Cheng; Nicholas J. Higham; Charles S. Kenney; Alan J. Laub
errShare
errSave
err1997-01-01
err0
PREAI
errNicholas J. Higham
errShare
errSave
errShare
errSave
Convergence Analysis of the Doubling Algorithm for Several Nonlinear Matrix Equations in the Critical Case
err2009-01-01
err0
PREAI
errChun-Yueh Chiang; Eric King-Wah Chu; Chun-Hua Guo; Tsung-Ming Huang; Wen-Wei Lin; Shu-Fang Xu
errShare
errSave
researcher View more