Pythonによる行列演算について

Question

Pythonの行列演算について質問です。

Pythonでは、Numpy, Scipy のサブパッケージであるlinalgを使えば行列の逆行列を求めることができますが(linalg.inv) 、実際に動いてるプログラムは、行列の形に対してある程度最適されたアルゴリズムで動くのでしょうか。それとも、掃き出し法などの計算時間はかかるが一般に解ける方法で解いているのでしょうか。その場合、自分でLU分解などの前処理をするのが一般的なのでしょうか。
ご回答お待ちしております。

Answer 1

数値解析は苦手なのですが…

scipyは線形代数演算ライブラリとして BLAS/LAPACK を使用してります。

linalg.inv() 関数の実装ですが、ソースコードを見る限り LAPACKのgetrf()とgetri()がコールされております。

getrf は LU分解を行う関数
getri は LU分解された行列から逆行列を求める関数

のようですので、LU分解により逆行列を求めているようですね。

Comments

[atnet]

ありがとうございます。大変参考になりました。