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

Question

有人能给我看一个c++代码示例，告诉我如何调用gsl函数gsl_linalg_LU_decomp（）并获得矩阵的逆吗？非常感谢！

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项目，通常在网上都有很好的文档，所以我建议您花点时间仔细阅读一下。