Page 1 of 1

LAPACKE_dgels : a way to force receiving positive solution ?

PostPosted: Thu Jul 19, 2018 9:47 am
by 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.