https://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.linalg.svd... 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.