//! IP-PMM (Pougkakiotis & Gondzio 2021, DOI 10.1007/s10589-020-00240-9)
//!
//! Augmented KKT (quasi-definite, upper-tri CSC):
//! K = [(Q + ρI), Aᵀ ]
//! [A, -D ] D = Σ + δI, Σ = diag(s/y)
//!
//! PMM update rule (Algorithm PEU §5.1.4):
//! r = |μ_k − μ_{k+1}| / μ_k (実 μ)
//! primal_improved = 0.95·prev_nr_p > nr_p → y_ref=y, δ *= (1−r), else δ *= (1−r/3)
//! dual_improved = 0.95·prev_nr_d > nr_d → x_ref=x, ρ *= (1−r), else ρ *= (1−r/3)
pub use solve_ippmm_inner;