为了进行分析，我需要求解 Ax=b 形式的大型稀疏系统。 A 的维度为 3Mx3M 并且是正定的。
目前，我使用“矩阵”包，它提供对稀疏矩阵和工作求解器的支持。我设法在大约 40 秒内解决了我的问题。但是，目前我服务器的 32 个核心中有 31 个在 mock 我，因为我只能使用一个核心进行计算。
我尝试制作一个 foreach 例程，但这非常占用 RAM/开销。我宁愿使用具有集成多核支持的 BLAS/LAPACK 求解器之类的东西。 Microsoft R Open 看起来很有前途，但它只支持密集矩阵的 BLAS/LAPACK。我广泛地尝试了谷歌搜索，但无济于事。
是否有可访问的多核求解器可用于求解 R 中的稀疏矩阵？
非常感谢任何帮助。
最佳答案
如果你愿意用 BLAS 库重新编译你的 R，你可以看看这样的东西:Using Intel MKL with R
关于r - R 中稀疏矩阵的多核求解，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/38462045/