Non-Negative Conjugate Gradients

Thomas Schmelzer, Martin Stoll

The conjugate gradient method (CG) solves a symmetric positive definite (SPD) system $Ax=b$, but does not respect $x\ge 0$. We develop non-negative conjugate gradients, a solver fo…