<div dir="ltr"><div>Jake Vanderplas also followed up on Fabian's post around those times. For reference, <a href="https://jakevdp.github.io/blog/2012/12/19/sparse-svds-in-python/">https://jakevdp.github.io/blog/2012/12/19/sparse-svds-in-python/</a></div><div><br></div><div>My guess is that PROPACK's license was ambiguous for a long time while sparse library was being implemented and getting mature. There are more experienced devs that can comment on that. <br></div><div><br></div><div>But maybe during the numpy.matrix revamp (<a href="https://github.com/scipy/scipy/issues/8162">https://github.com/scipy/scipy/issues/8162</a>) we can consider involving PROPACK alternatives. <br></div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jun 15, 2018 at 1:55 AM, Robert Lucente <span dir="ltr"><<a href="mailto:rlucente@pipeline.com" target="_blank">rlucente@pipeline.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><a href="https://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.linalg.svds.html" rel="noreferrer" target="_blank">https://docs.scipy.org/doc/<wbr>scipy/reference/generated/<wbr>scipy.sparse.linalg.svds.html</a><br>
<br>
states "This is a naive implementation using ARPACK as an eigensolver on A.H * A or A * A.H, depending on which one is more efficient."<br>
<br>
One of my co-workers is struggling to find an implementation for SVD for a sparse matrix. I was thinking of suggesting that he look at SciPy. The phrase "naive implementation" has made me think twice about making the suggestion. I just don't want the first intro of Python/SciPy into the organization to be a negative one.<br>
<br>
I poked around and found the following 2 articles<br>
<br>
<a href="http://fa.bianp.net/blog/2012/singular-value-decomposition-in-scipy/" rel="noreferrer" target="_blank">http://fa.bianp.net/blog/2012/<wbr>singular-value-decomposition-<wbr>in-scipy/</a><br>
<br>
<a href="https://simplyml.com/benchmarking-the-singular-value-decomposition/" rel="noreferrer" target="_blank">https://simplyml.com/<wbr>benchmarking-the-singular-<wbr>value-decomposition/</a><br>
<br>
I realize that I should be suspicious of bench marks and that you have to run your own.<br>
<br>
I would like to provide my co-worker the link to SciPy and then to 2 articles which can help guide him in creating his own bench mark. Any suggestions would be appreciated.<br>
<br>
Please note that I am not an expert in SVD.<br>
______________________________<wbr>_________________<br>
SciPy-Dev mailing list<br>
<a href="mailto:SciPy-Dev@python.org">SciPy-Dev@python.org</a><br>
<a href="https://mail.python.org/mailman/listinfo/scipy-dev" rel="noreferrer" target="_blank">https://mail.python.org/<wbr>mailman/listinfo/scipy-dev</a><br>
</blockquote></div><br></div>