Positive definiteness in the numerical solution of Riccati differential equations