使用增强解决（密集）线性系统Ax = b

Question

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

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

Answer 1

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

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