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

Question

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

Answer 1

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

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

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