A block-shifted cyclic reduction algorithm for solving a class of quadratic matrix equations

The cyclic reduction (CR) algorithm is an efficient method for solving quadratic matrix equations that arise in quasi-birth–death (QBD) stochastic processes. However, its convergence is not guaranteed when the associated matrix polynomial has more than one eigenvalue on the unit circle. To address this limitation, we introduce a novel iteration method, referred to as the Block-Shifted CR algorithm, that improves the CR algorithm by utilizing singular value decomposition (SVD) and block shift-and-deflate techniques. This new approach extends the applicability of existing solvers to a broader class of quadratic matrix equations. Numerical experiments demonstrate the effectiveness and robustness of the proposed method.