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

Question

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

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

Answer 1

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

如果您需要代码来执行此操作，您可以查看数字配方作为示例。但是，解决稠密线性系统是一个非常丰富的主题，所以它可能超出了这里的范围来解决它的所有方面。