Abstract
New theoretical results are presented about the principal matrix pth root. In particular, we show that the pth root is related to the matrix sign function and to the Wiener-Hopf factorization, and that it can be expressed as an integral over the unit circle. These results are used in the design and analysis of several new algorithms for the numerical computation of the pth root. We also analyze the convergence and numerical stability properties of Newton's method for the inverse pth root. Preliminary computational experiments are presented to compare the methods. © Springer 2005.
Original language | English |
---|---|
Pages (from-to) | 349-378 |
Number of pages | 29 |
Journal | Numerical Algorithms |
Volume | 39 |
Issue number | 4 |
DOIs | |
Publication status | Published - Aug 2005 |
Keywords
- Cyclic reduction
- Graeffe iteration
- Laurent polynomial
- Matrix pth root
- Matrix sign function
- Newton's method
- Wiener-Hopf factorization