LAPACKE_dgels : a way to force receiving positive solution ?

[jasonchampion]

Hi everyone,

I am currently using LAPACKE to take advantages of LAPACK within my C program. I am using the function LAPACKE_dgels in order to solve an underestimated (or sometimes overestimated ^^) system of linear equation. It works! But... if the solution than I receive among the infinity available, it contains a few negative values which is inconsistent and unphysical for my case. Thus my question is the following:
Is there a way to force that routine, or to find another one, allowing me to receive a solution containing only zero or positive values ?

Thank you in advance for your help!

Jason.