[SciPy-dev] sparse comments and questions
Fri Jul 13 13:36:51 CDT 2007
Thanks for putting this in, I'll give it a try.
On 7/13/07, Nathan Bell <firstname.lastname@example.org> wrote:
> On 7/9/07, Peter Skomoroch <email@example.com> wrote:
> > In Matlab and some C++ libraries, I've used the following "sparse
> > functionality :
> > When you attempt to divide one sparse matrix by another of the same
> > return the result of elementwise division of the entries in the
> > This assumes that the two input matrices have the same sparsity
> > This is useful when implementing some algorithms like NMF using sparse
> > matrices. Right now in scipy, only division of a sparse matrix by a
> > is supported.
> I've added support for elementwise division in the latest SVN. I
> designed it so you get the same result as a dense A/B on the union of
> the sparsity structures. For example, if A[i.j] is nonzero while
> B[i.j] is zero, you get +/- inf. If A[i,j] is zero while B[i,j] is
> nonzero you'll get an implicit 0 in the result (i.e. it will not
> appear in the output matrix). This holds true for explicit zeros in
> the matrices as well. I think this is the most faithful analog of the
> dense case.
> Let me know if you discover any problems.
> Nathan Bell firstname.lastname@example.org
> Scipy-dev mailing list
Peter N. Skomoroch
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Scipy-dev