[SciPy-User] Matrix Exponentials For Very Large Sparse Matrices
Fri Jan 22 18:38:08 CST 2010
On Fri, Jan 22, 2010 at 7:17 PM, Dylan Gorman <firstname.lastname@example.org> wrote:
> Hi Folks,
> I'd like to exponentiate very large sparse matrices with scipy, and I
> would be very grateful for any suggestions. Currently, I'm running
> into memory errors exponentiating random matrices of order 10^3x10^3
> with the standard linalg.expm() routine. However, I suspect that I may
> realistically be able to handle somewhat larger matrices since the
> actual matrices I will be using are quite sparse. Ideally, I'd like to
> be able to exponentiate matrices of size 10^5-10^6 x 10^5 - 10^6.
> However, there does not seem to be any linalg.sparse.expm() function--
> is this because there is in fact no advantage to exponentiating sparse
> matrices? Or would I need to implement something by hand?
Just an idea until the experts come along.
I would try for a sparse eigenvector decomposition and then only the
eigenvalues need to be exponentiated.
However, I never tried it with sparse, and I don't know how sparse the
eigenvectors would be.
> Thank you very much,
> Dylan Gorman
> SciPy-User mailing list
More information about the SciPy-User