Abstract
A numerical method for computing the square root of a symmetric positive definite matrix is developed in this paper. It is based on the Padé approximation of √1 + x in the prime fraction form. A precise analysis allows us to determine the minimum number of terms required in the Padé approximation for a given error tolerance. Theoretical studies and numerical experiments indicate that the method is more efficient than the standard method based on the spectral decomposition, unless the condition number is very large.
| Original language | English |
|---|---|
| Pages (from-to) | 833-845 |
| Journal | SIAM Journal on Matrix Analysis and Applications |
| Volume | 19 |
| Issue number | 3 |
| DOIs | |
| Publication status | Published - Jul 1998 |
Research Keywords
- Matrix square root
- Padé approximation
- Prime fraction form
Publisher's Copyright Statement
- COPYRIGHT TERMS OF DEPOSITED FINAL PUBLISHED VERSION FILE: © 1998 Society for Industrial and Applied Mathematics