Numerical stability and instability in matrix sign function based algorithms

This paper uses a forward and backward error analysis to try to identify some classes of matrices for . P which the matrix sign function is a numerically stable algorithm for extracting invariant subspaces roper scaling is essential to numerical stability as well as to rapid convergence. g a Roberts [21] and Beavers and Denman [7] introduced the matrix sign function as a means of solvin lgebraic Riccati equations and Lyapunov equations. The matrix sign function has since attracted the ] a attention of control engineers and some applied mathematicians ([1] to [21]). Balzer [3], Barraud [5 nd Byers [9] have suggested strategies for accelerating convergence. Denman and Beavers [11] - - 0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000002 extended matrix sign function algorithms to a list of invariant subspace related calculations. Howland e a [16] used the matrix sign function to count eigenvalues in boxes in the complex plane. Some of th lgorithms have been refined and extended by Attarzadeh [2], Bierman [8], Byers [9]. Gardiner and d d Laub [14] have extended the use of the matrix sign function to generalized Riccati equations an iscrete Riccati equations. Higham [15] has used matrix sign function techniques to calculate polar decomposition