좋은 행렬 역변환 루틴 C

Question

수치 계산을위한 파이썬 코드의 일부로, 다소 큰 (희소) 행렬 (~ 100x100)을 여러 번 뒤집어 야합니다. 프로그램을 빠르게하고 싶습니다. 그리고 나에게 제안 된 방법 중 하나는 행렬 반전 단계에서 C로 서브 루틴을 호출하는 것입니다.

이 작업을 위해 효율적이고 잘 테스트 된 C 루틴이 있습니까?

감사합니다.

Answer 1

>>> from numpy import * 
>>> from numpy.linalg import inv 
>>> from scipy.sparse import csr_matrix 
>>> m = matrix([[3,1,5],[1,0,8],[2,1,4]]) 
>>> s = csr_matrix(m) 
>>> invs = inv(a) # Inverse sparse matrix 
>>> dot(a,inva) # Check the result, should be eye(3) within machine precision 
csr_matrix([[ 1.00000000e-00, 2.77555756e-17, 3.60822483e-16], 
      [ 0.00000000e+00, 1.00000000e+00, 0.00000000e+00], 
      [ -1.11022302e-16, 0.00000000e+00, 1.00000000e+00]]) 

정말로 필요한 역수입니까? 당신은 반전없이 목표를 달성 할 수 있습니다 : 당신이 정말이 역을 필요로

사례는 드물다. 또한 희소 행렬의 역함수가 반드시 희소하지는 않습니다. 일반적으로 반전은 LU 인수 분해보다 비용이 높으며 반올림 오류가 발생하기 쉽습니다.

-http://mail.scipy.org/pipermail/scipy-user/2007-October/013936.html

->http://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.linalg.factorized.html#scipy.sparse.linalg.factorized