使用R中的多個核心計算SVD

Question

我想在R大型稀疏矩陣（17k x 2m）上運行svd()，並且我可以訪問羣集。有沒有一種直接的方法可以用多核來計算R中的SVD？

RScaLAPACK軟件包（http://www.inside-r.org/packages/cran/RScaLAPACK）似乎使這成爲可能，但它似乎不再支持（http://cran.r-project.org/web/packages/RScaLAPACK/），我認爲這是有原因的。

Answer 1

rARPACK是您需要的軟件包。像魅力一樣工作（即使矩陣比您的規格大得多）。超快速，因爲它通過C和C++並行化。

Answer 2

rARPACK是一種選擇，但要確保你有一個優化的多核BLAS庫，因爲所有的並行計算部分本身並不在rARPACK之內，而是在BLAS中。

另外，要小心，rARPACK只計算局部 SVD，這意味着它僅計算最大k奇異值和相關的奇異向量。如果您確實需要完整的SVD，您仍然可以使用svd()。