Re: [Numpy-discussion] lstsq underdetermined behaviour

2018-11-18 Thread Romesh Abeysuriya
Thanks both! Yes, I guess it's typically 'least squares' referring to the residual vector, and 'minimum norm' referring to the solution vector. That's certainly how the documentation for `dgelsd` frames it. In my case, the minimum norm solution can be sensibly interpreted (and in particular, it gua

[Numpy-discussion] lstsq underdetermined behaviour

2018-11-18 Thread Romesh Abeysuriya
Hi all, I'm solving an underdetermined system using `numpy.linalg.lstsq` and trying to track down its behavior for underdetermined systems. In previous versions of numpy (e.g. 1.14) in `linalg.py` the definition for `lstsq` calls `dgelsd` for real inputs, which I think means that the underdetermin