The Bound-Constrained Conjugate Gradient Method for Non-negative Matrices