解決（密集）線性系統Ax = b使用boost

Question

我有一個密集的Ax = b類型方程組來解決我的C++程序，並且我希望在boost中使用UBLAS實現解決方案。在其他一些問題中，我發現人們使用擴展LAPACK，但不幸的是，它似乎不是我的標準增強安裝的一部分（至少在Debian中），我不允許添加更多的依賴關係。

有人可以粘貼不使用LAPACK或任何其他庫的解決方案嗎？

Answer 1

不幸的是，你正在解決一個線性系統，要麼需要LAPACK，要麼編寫自己的代碼。如果你不想LAPACK，你唯一的選擇是編寫自己的求解器。這樣的求解器當然可以使用uBLAS。

如果您需要執行此操作的代碼，您可以查看numerical recipes的示例。但是，解決密集的線性系統是一個非常豐富的主題，所以它可能超出瞭解決它的所有方面的範圍。