A low-rank Krylov squared Smith method for large-scale discrete-time Lyapunov equations

The squared Smith method is adapted to solve large-scale discrete-time Lyapunov matrix equations. The adaptation uses a Krylov subspace to generate the squared Smith iteration in a low-rank form. A restarting mechanism is employed to cope with the increase of memory storage of the Krylov basis. Theoretical aspects of the algorithm are presented. Several numerical illustrations are reported.