I want to solve X in a linear system A*X=B with constraints that all elements in X are non-negative. It is a Non-negative Least Squares (NNLS) problem. Do you know which function can be used for it? I tried the function trnlspbc_solve (nonlinear least squares problem with linear constraints), but it does not work well for this linear least square problem. Any ideas?