[SciPy-User] Sparse Matricies and NNLS

[Calvin Morrison]

Paul,

It seems nobody wants to touch the nnls algorithm because the only
implementation that is floating around is the one from the original
publication or automatic conversions of it.

Calvin
On Mar 28, 2013 4:48 PM, "Pauli Virtanen" <pav at iki.fi> wrote:

> 28.03.2013 21:46, Calvin Morrison kirjoitti:
> [clip]
> > I was thinking about using a sparse matrix since it is relatively
> > sparse, but the problem is that even if I was to use the sparse
> > matricies, the nnls function only accepts a ndarray, not a sparse
> > matrix. (when I try and throw a sparse matrix at it I get an error.)
>
> The nnls algorithm in Scipy relies on dense matrix algebra, and is
> moreover written in Fortran. So, there is no way to tell it to use
> sparse matrices.
>
> You'll need to find an implementation of NNLS algorithm that either is
> matrix-free or works for sparse problems. If you find such code, be sure
> to reply to this list --- it might be useful to include it in Scipy,
> provided the license is compatible.
>
> --
> Pauli Virtanen
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20130328/9202d8e0/attachment.html>