データの性質によるのかも知れませんが、次のように一様分布の乱数行列だと、結果はほぼ1,1になりました。
import numpy as np
A = np.random.rand(400, 400)
Ainv = np.linalg.inv(A)
AAi = np.dot(A, Ainv)
AiA = np.dot(Ainv, A)
for i in range(400):
print(AAi[i][i], AiA[i][i])
結果
(0.999999999999984, 0.9999999999998891)
(1.0000000000000198, 0.9999999999999999)
(0.9999999999999643, 1.0000000000002853)
(1.0000000000000118, 1.0000000000000686)
(1.0000000000001563, 0.9999999999999708)
(1.000000000000004, 1.0000000000000784)
以下略
ただ、inv()では精度が出ないと言う話は聞いたことがあり、その場合solve()を使えとのことです。
Ainv = np.linalg.solve(A, np.eye(400))
これでもほぼ1,1が得られました。なお、np.show_config()は次の通り。実行時ではなくコンパイル時の設定らしいので、参考になるかわかりません。lapack, blas, openblas, atlasなどとりあえず全部インストールした状態です。
lapack_info:
libraries = ['lapack', 'lapack']
library_dirs = ['/usr/lib']
language = f77
lapack_opt_info:
libraries = ['lapack', 'lapack', 'blas', 'blas']
library_dirs = ['/usr/lib']
language = c
define_macros = [('NO_ATLAS_INFO', 1), ('HAVE_CBLAS', None)]
openblas_lapack_info:
NOT AVAILABLE
blas_info:
libraries = ['blas', 'blas']
library_dirs = ['/usr/lib']
define_macros = [('HAVE_CBLAS', None)]
language = c
atlas_3_10_blas_threads_info:
NOT AVAILABLE
atlas_threads_info:
NOT AVAILABLE
atlas_3_10_threads_info:
NOT AVAILABLE
atlas_blas_info:
NOT AVAILABLE
atlas_3_10_blas_info:
NOT AVAILABLE
atlas_blas_threads_info:
NOT AVAILABLE
openblas_info:
NOT AVAILABLE
blas_mkl_info:
NOT AVAILABLE
blas_opt_info:
libraries = ['blas', 'blas']
library_dirs = ['/usr/lib']
language = c
define_macros = [('NO_ATLAS_INFO', 1), ('HAVE_CBLAS', None)]
atlas_info:
NOT AVAILABLE
atlas_3_10_info:
NOT AVAILABLE
lapack_mkl_info:
NOT AVAILABLE
mkl_info:
NOT AVAILABLE
np.linalg.cond(A)
の結果を教えて頂けませんでしょうか。np.linalg.cond(A)
の結果は3.16e+18
です。A = A[:10,:10]
とすると1366. 5
,A = A[:100,:100]
とすると1.03e+11
などとなっています。np.linalg.cond
の値が1e+19前後になるような行列を作成してnp.linalg.inv
をやって見ましたが、AAi
やAiA
がそこまでずれることはなかったです。gistかpastebinあたりにそのようになる行列の例を載せることは可能でしょうか。