[SciPy-Dev] least squares solvers

Sturla Molden sturla at molden.no
Fri Jan 27 04:47:15 EST 2012


Why is the non-linear LS solver in sp.optimize called leastsq, whereas 
the linear solver in sp.linalg is called lstsq? Wouldn't a consistent 
name be better?

sp.linalg.lstsq uses SVD, by lapack driver *gelss. If we don't need 
singular values, solving by QR is faster (lapack driver *gels).

(I actually use Fortran just to get DGELS instead of DGELSS.)


Sturla



More information about the SciPy-Dev mailing list