[SciPy-Dev] scipy.sparse.linalg.svds

[Robert Lucente]

https://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.linalg.svds.html

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."

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.

I poked around and found the following 2 articles

http://fa.bianp.net/blog/2012/singular-value-decomposition-in-scipy/

https://simplyml.com/benchmarking-the-singular-value-decomposition/

I realize that I should be suspicious of bench marks and that you have to run your own.

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.

Please note that I am not an expert in SVD.