Python module that solves projected Lyapunov and Riccati equations via saddle point systems and ADI