求解大型稠密线性方程组的R函数？

对不起，我可能是瞎子，但我找不到解决一个相当常见问题的具体方法：

我想实施

solve(A,b)

与 A. 作为一个大的方阵，上面的命令使用了我所有的内存并发出一个错误（b是一个具有相应长度的向量）。我得到的矩阵不是稀疏的，因为会有大量的零块等。 一定有一些函数实现了逐步迭代方案，这样即使可用内存有限也能找到解决方案

我在稀疏矩阵上找到了好几篇文章，当然还有矩阵包，但无法找到一个满足我需要的函数。我也看到了，但是 比格姆 生成完整的线性模型拟合。我只需要一个简单的解决方案。我将不得不重复这一步几次，所以它将是伟大的保持尽可能苗条

---

我已经担心“老问题的重复”和“看这里”的评论，但我真的非常感谢您的帮助。

如果您能找出您想要的算法的名称，将更容易看到它是否在较低级别的库中实现，如LAPACK、LINPACK、Armadillo。。。然后查找R包装（例如在

矩阵

或

RCppArmadillo

包中）。“核心外”似乎很有用：@flodel:假设它足够大和足够密集，正如我上面所描述的：它使用了所有的内存，稀疏矩阵方法没有帮助。更精确地说：取一个实际上没有零元素的10kx10k矩阵。应用

solve（）！我查看了一下，但还没有在R中找到实现。矩阵包的文档显然没有提到“内存”这个词，RccpArmadillo需要C++编程，而我还没有能力。我是否监督过任何直接在plain R中实现这种方法的东西？查看High Performance Task视图会得到
pbdDMAT
包，它可能会满足您的需要。我还没有深入研究。