对于分析，我需要以 Ax = b 的形式求解大型稀疏系统。 A 的维度为 3Mx3M 并且是肯定的。
目前，我正在使用＆＃39; Matrix＆＃39;包，它提供对稀疏矩阵和工作求解器的支持。我设法在大约40秒内解决我的问题。但是，目前我服务器的32个核心中有31个在嘲笑我，因为我只能用一个核心进行计算。
我尝试制作一个foreach例程，但这是非常密集的RAM /开销。我宁愿使用类似BLAS / LAPACK求解器和集成多核支持的东西。 Microsoft R Open似乎很有希望，但它只支持密集矩阵的BLAS / LAPACK。我广泛尝试谷歌搜索无济于事。
可访问的多核求解器是否可用于求解R？
中的稀疏矩阵任何帮助都非常适合。