c++ - 使用 boost 求解(密集)线性系统 Ax=b

Question

我有一个密集的 Ax=b 类型的方程组要在我的 C++ 程序中求解，我希望在 boost 中使用 UBLAS 来实现该解决方案。在其他一些问题中，我发现人们正在使用扩展 LAPACK，但不幸的是，它似乎不是我标准 boost 安装的一部分(至少在 Debian 中)，而且我不允许添加更多依赖项。

有人可以粘贴不使用 LAPACK 或任何其他库的解决方案吗？

Answer 1

不幸的是，您正在求解一个需要 LAPACK 或编写您自己的代码的线性系统。如果你不想要 LAPACK，你唯一的选择就是编写你自己的求解器。这样的求解器当然可以使用 uBLAS。

如果需要代码可以看numerical recipes举个例子。但是，求解稠密线性系统是一个非常丰富的主题，因此可能超出了此处的范围来解决它的所有方面。