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

Question

有人可以告诉我一个关于如何调用gsl函数gsl_linalg_LU_decomp（）以及获取矩阵的逆的示例c ++代码吗？非常感谢！

Answer 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项目，通常在线有很好的文档，所以我建议你花点时间仔细阅读它。