问R中稀疏矩阵的多核求解

EN

为了进行分析，我需要求解一个Ax=b形式的大型稀疏系统，a的维数为3Mx3M，且为正定。

现在，我使用了“矩阵”包，它提供了对稀疏矩阵和工作求解器的支持。我设法在40秒内解决了我的问题。然而，目前我的服务器的32个核心中有31个在嘲笑我，因为我只能用一个核心来计算。

我试着做了一个前端例程，但是这非常占用内存/开销。我宁愿使用像BLAS/LAPACK这样的集成多核支持的解决程序。Microsoft似乎很有前途，但它只支持BLAS/LAPACK用于密集矩阵。我试过大量搜索，但没有结果。

一个可访问的多核求解器是否可用于求解R中的稀疏矩阵？

任何帮助都是非常有用的。