I ran into this problem as well a few months back.
The reason for the empty residual array when M==N is that the LAPACK
routine for Ax = b puts the solution for x in b. When M>N, the
norm-squared is parceled out into the unused (M-N) points in the b
array. When M==N, there's no room for the resids
Hi all,
The summary of ticket 937 is incomplete.
It should be "Complex matrices and lstsq".
http://projects.scipy.org/scipy/numpy/ticket/937
Nils
___
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/nu
