我有一个Lanscos实现的降秩SVD，它可以分布并处理具有数亿行和列的非常稀疏的矩阵。可以想象，有时需要几天的时间。
我最近在https://research.fb.com/fast-randomized-svd/遇到了Facebook快速随机SVD，代码在http://tygert.com/software.html，我想尝试一下，但是fortran代码仅支持密集矩阵输入，在我看来我无法使用它。
那么，有人知道我可以使用的此代码或类似代码的版本（最好是分布式fortran）吗？
或者，由于我发现的代码包括源代码，有没有一种简单的方法可以从密集矩阵代码创建稀疏矩阵版本？
FWIW，该代码基于BLAS和LAPACK构建。
感谢您提供的任何帮助！