“…Direct methods are widely used for solving the KKT systems in well-established optimization codes based on IP methods, such as MOSEK [61], LOQO [77,78], OOQP [37], KNITRO-Interior/Direct [12], IPOPT [79]. In this case, modified Cholesky factorizations [19,40,72], that are capable of handling the case the matrix is not positive definite, are usually applied to the condensed system, while LBL T factorizations, where L is unit lower triangular and B is symmetric block diagonal with 1 × 1 or 2 × 2 blocks, are applied to the KKT system.…”