Numerically stable LDLT factorizations in interior point methods for convex quadratic programmingD. GoldfarbK. Scheinberg2008IMA Journal of Numerical Analysis