需要一個使用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]); 

如果你確實需要一個實際的倒數，然後按照我交聯以獲得LU分解，然後調用函數gsl_linalg_LU_invert的例子。 gsl圖書館是一個GNU項目，通常在網上有很好的文檔，所以我建議你花一些時間仔細閱讀一下。