まあ,6 でちょっとやってみたけれど。
from numpy import * from scipy import io, sparse A = sparse.lil_matrix((3, 3)) A[0,1] = 3 A[1,0] = 2 A[2,2] = 5 |
こんな感じ?
solve=scipy.sparse.linalg.factorized(a) x1=solve(b1) x2=solve(b2) ... xn=solve(bn) |
固有値問題もとりあえず大丈夫そう。 例の問題をどれくらいの効率で解けるかが気になる。
桂田 祐史