gsl LU分解を使用して正方行列の逆行列を得るためのサンプルコードが必要

Question

gsl関数gsl_linalg_LU_decomp（）を呼び出す方法と行列の逆行列を取得する方法に関するサンプルC++コードを親切に教えてもらえますか？とても感謝しています！

Answer 1

実際の逆関数は必要ないと仮定していますが、Ax = bという型の問題を解決する必要があります。もしそうなら、かなり良い例があるでしょうhere。データにSTLコンテナを使用している場合（例： std::vector、その後、あなたはAの実際の逆を必要とした場合は、その後、私はLU分解を取得し、機能gsl_linalg_LU_invertを呼び出すためにリンクされている例に従う

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

のような最初のデータエントリへのポインタを渡す必要があります。 gslライブラリはGNUプロジェクトであり、一般的にオンラインでよく文書化されていますので、ちょっと読んでみることをお勧めします。