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 プロジェクトであり、一般的にオンラインで十分に文書化されているため、時間をかけて少し読んでみることをお勧めします。