Symplectic Möbius integrators for LQ optimal control problems

The paper presents symplectic Möbius integrators for Riccati equations. All proposed methods preserve symmetry, positivity and quadratic invariants for the Riccati equations, and non-stationary Lyapunov functions. In addition, an efficient and numerically stable discretization procedure based on reinitialization for the associated linear Hamiltonian system is proposed.