GPU acceleration of splitting schemes applied to differential matrix equations

We consider differential Lyapunov and Riccati equations, and generalize