Numerical considerations in computing invariant subspaces
SIAM Journal on Matrix Analysis and Applications
Matrix computations (3rd ed.)
Solution of the matrix equation AX + XB = C [F4]
Communications of the ACM
Accuracy and Stability of Numerical Algorithms
Accuracy and Stability of Numerical Algorithms
Hi-index | 0.00 |
This paper presents a new numerical algorithm for solving the Sylvester equation involved in higher-order perturbation methods developed for solving stochastic dynamic general equilibrium models. The new algorithm surpasses other methods used so far (including the very popular doubling algorithm) in terms of computational time, memory consumption, and numerical stability.