{"id":"arxiv/2607.22121","url":"https://arxiv.org/abs/2607.22121","tags":["math.OC","math.NA"],"score":0,"title":"Non-Negative Conjugate Gradients","domain":"math.OC","source":"arxiv","authors":["Thomas Schmelzer","Martin Stoll"],"reviews":0,"subject":"math.OC","summary":"The conjugate gradient method (CG) solves a symmetric positive definite (SPD) system $Ax=b$, but does not respect $x\\geq0$. We develop non-negative conjugate gradients, a solver for the bound-constrained quadratic $\\min_x \\tfrac12 x^\\top A x - b^\\top x$ subject to $x\\geq0$, $Bx=c$, $A\\succ0$, with $Bx=c$ present only in the equality-augmented variant. The solver wraps CG in a primal-dual active-set loop: each step solves the unconstrained system over the free variables, drops variables that turn negative, re-admits any bound variable with negative reduced gradient, and restarts CG. These…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-07-24","reactions":{}}