需要一个示例代码来使用 gsl LU 分解获得方阵的逆

【问题标题】：Need an example code for getting the inverse of a square matrix using gsl LU decomposition需要一个示例代码来使用 gsl LU 分解获得方阵的逆
【发布时间】：2012-03-16 03:56:30
【问题描述】：

有人可以向我展示一个关于如何调用 gsl 函数 gsl_linalg_LU_decomp() 并与获取矩阵逆相关的示例 C++ 代码吗？非常感谢！

【问题讨论】：

标签： c++ linear-algebra gsl

【解决方案1】：

我假设您不需要实际的逆，但您需要解决 Ax=b 类型的问题。如果是这样，那么有一个很好的例子here。如果您使用 STL 容器存储数据，例如std::vector，那么你需要传递一个指向第一个数据条目的指针，比如

std::vector<double> vec(length,val);
gsl_needs_ptr_to_double(&vec[0]);

如果您确实需要 A 的实际逆，请按照我链接的示例获取 LU 分解，然后调用函数gsl_linalg_LU_invert。 gsl 库是一个 GNU 项目，通常在网上都有很好的文档记录，所以我建议您花点时间阅读一下。

【讨论】：