[SciPy-User] Matrix Exponentials For Very Large Sparse Matrices

[Charles R Harris]

On Fri, Jan 22, 2010 at 5:17 PM, Dylan Gorman <dgorman at berkeley.edu> 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?
>
>
Out of curiosity, do the matrices have any special structure? For instance,
are they banded or symmetric? Also, why to you want to exponentiate them?

Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20100122/0ae14492/attachment.html>