[Numpy-discussion] Linear least squares
Tue Jan 8 16:35:12 CST 2013
On Tue, Jan 8, 2013 at 1:17 PM, Till Stensitz <firstname.lastname@example.org> wrote:
> i did some profiling and testing of my data-fitting code.
> One of its core parts is doing some linear least squares,
> until now i used np.linalg.lstsq. Most of time the size
> a is (250, 7) and of b is (250, 800).
My guess is that this depends a lot on the shape
try a is (10000, 7) and b is (10000, 1)
> Today i compared it to using pinv manually,
> to my surprise, it is much faster. I taught,
> both are svd based? Too check another computer
> i also run my test on wakari:
> Also using scipy.linalg instead of np.linalg is
> slower for both cases. My numpy and scipy
> are both from C. Gohlkes website. If my result
> is valid in general, maybe the lstsq function
> should be changed or a hint should be added
> to the documentation.
> NumPy-Discussion mailing list
More information about the NumPy-Discussion