Muszę odwrócić dużą macierz rzadką. Nie mogę uciec od inwersji matrycy, jedynym skrótem byłoby po prostu zorientowanie się w głównych przekątnych elementach i zignorowanie elementów poza diagonalnych (wolałbym nie, ale jako rozwiązanie byłoby do przyjęcia).Odwracanie dużych rzadkich macierzy za pomocą scipy
Matryce muszę odwracające są zazwyczaj duże (40000 * 40000) i które mają tylko kilka non-niezerowych przekątnych. Moje obecne podejście jest budowanie wszystko skąpe, a następnie
posterior_covar = np.linalg.inv (hessian.todense())
to wyraźnie zajmuje dużo czasu i dużo pamięci.
żadnych wskazówek, czy jest to tylko kwestia cierpliwości lub podejmowania problemu mniejsze?
Opcja 'scipy' doc mówi, że nie trzeba do zagęszczenia matrycę, więc jestem nieco zdezorientowany. http://docs.scipy.org/doc/scipy-0.8.x/reference/sparse.html – BenDundee
Wersja 0.12 scipy (aktualnie w fazie testów beta) ma funkcję scipy.sparse.linalg.inv. –