Hi. I am solving an under-determined system with DGELS. The matrix A has dimensions MxN, with M<N. This matrix is sparse, but is non diagonal, the elements appear at some determined places.

I wanted to know if there is a way to tell lapack that this matrix is sparse, in order to not having to initialize all the entries that are zero, saving memory space, and avoiding trivial computations like a(9,9)*zero. What I would like to do is to entry only the non zero elements, and then lapack do the thing. Is there a way to do this?

Bye there!