Need a Lapack routine for AX=B (A not nxn)

[ctcsback]

I need something like DGESV for matrix data of doubles. However DGESV doesn't work because I need a solver that can solve AX=B where A can be any n*m matrix.

So basically I have:

A matrix : nxm
B matrix : nx1

and I need to find X : mx1

I feel that this should have a simple solution. Any suggestions?

[Julien Langou]

Please check out dgels. J.