[SciPy-Dev] Iterative Least Squares for Sparse Systems

[Stéfan van der Walt]

Hi Jake

Lsqr will be included in the next release, but in the meantime you may
download lsqr.py and import it as a stand-alone module.

Lsqr is aimed at solving large, sparse, over-determined systems such as
yours.

Regards
Stefan

On Apr 14, 2010 6:58 PM, "Jake VanderPlas" <jakevdp at gmail.com> wrote:

Hello,
I've noticed a trunk project at

http://projects.scipy.org/scipy/browser/trunk/scipy/sparse/linalg/isolve/lsqr.py
This looks like a very useful piece of code - anybody know the status
of including it in a scipy release?
I ask because I need to find a solution of a very large
over-determined sparse linear system.  I suppose I could just use
scipy.optimize.leastsq on a function which returns the residual, but
that seems a bit hacked.  Is there a better way to do this at present?
  -Jake
_______________________________________________
SciPy-Dev mailing list
SciPy-Dev at scipy.org
http://mail.scipy.org/mailman/listinfo/scipy-dev
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20100419/b661e99a/attachment.html>