Which function could solve linear least square with inequality constraints?

[Guo__Qiandong]

I want to solve X in a linear system A*X=B with a constraint that all elements in X are non-negative. It is a Non-negative Least Squares (NNLS) problem. Do you know which function in the MKL library for this problem? I tried to use trnlspbc_solve(), the Nonlinear Least Squares Problem with Linear Constraints, but it doesn't work well. Any ideas? 

[mecej4]

You may choose a suitable NNLS solver from the list provided at http://plato.asu.edu/sub/nonlsq.html , under the heading "THE CONSTRAINED LSQ-PROBLEM".

[Guo__Qiandong]

Thank you for your reply, mecej4! I know "Lawson, C.L. and R.J. Hanson, Solving Least Squares Problems, Prentice-Hall, 1974" could solve it. But I prefer to use MKL to do it. Does MKL have any solutions? When I use trnlspbc_solve(), it is time-consuming and the result is a little bit different than NNLS. Any comments about using trnlspbc_solve() to do linear least square problem with constraints? 

[mecej4]

At this point it would be useful to have a working example source code that exhibits slow performance with trnlspbc_solve(). Even a mathematical problem definition would help. I have used trnlsp() a number of times, and it has always performed quite well. Your problem may use specific functions and bounds that trigger a slowdown, which is why a full description is needed.