RADI: A low-rank ADI-type algorithm for large scale algebraic Riccati equations

Abstract

This paper introduces a new algorithm for solving large-scale continuous-time algebraic Riccati equations (CARE). The advantage of the new algorithm is in its immediate and efficient low-rank formulation, which is a generalization of the Cholesky-factored variant of the Lyapunov ADI method. We discuss important implementation aspects of the algorithm, such as reducing the use of complex arithmetic and shift selection strategies. We show that there is a very tight relation between the new algorithm and three other algorithms for CARE previously known in the literature -- all of these seemingly different methods in fact produce exactly the same iterates when used with the same parameters: they are algorithmically different descriptions of the same approximation sequence to the Riccati solution.

Cited by in corpus (9)

• A numerical comparison of solvers for large-scale, continuous-time algebraic Riccati equations and LQR problems
• Robust output-feedback stabilization for incompressible flows using low-dimensional ${\mathcal{H}}_{\infty }$-controllers
• Multiscale differential Riccati equations for linear quadratic regulator problems
• A low-rank solution method for Riccati equations with indefinite quadratic terms
• Parametric model order reduction using pyMOR
• On a family of low-rank algorithms for large-scale algebraic Riccati equations
• Residual-based iterations for the generalized Lyapunov equation
• The intrinsic Toeplitz structure and its applications in algebraic Riccati equations
• Deep polytopic autoencoders for low-dimensional linear parameter-varying approximations and nonlinear feedback design