Abstract
An efficient and numerically stable implementation of a known algorithm is suggested for finding the extremal positive definite solutions of the matrix equation X + A*X-1A = I, if such solutions exist. The convergence rate is analyzed. A new algorithm that avoids matrix inversion is presented. Numerical examples are given to illutrate the effectiveness of the algorithms.
| Original language | English |
|---|---|
| Pages (from-to) | 1167-1174 |
| Number of pages | 8 |
| Journal | SIAM Journal on Scientific Computing |
| Volume | 17 |
| Issue number | 5 |
| DOIs | |
| State | Published - Sep 1996 |
| Externally published | Yes |
Keywords
- Iteration
- Matrix equation
- Positive definite solution