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

Nils Wagner nwagner@iam.uni-stuttgart...
Sun Jan 24 03:11:19 CST 2010

On Fri, 22 Jan 2010 20:33:15 -0500
  Joshua Stults <joshua.stults@gmail.com> wrote:
> On Fri, Jan 22, 2010 at 8:24 PM, Dylan Gorman 
><dgorman@berkeley.edu> wrote:
>> Dear Chuck and Joshua,
>> It's a problem in quantum simulation. I'm trying to 
>>solve d(rho)/dt =
>> L*rho for a sparse matrix L. L should be symmetric, and 
>>in principle I
>> just need to compute (e^L*t)*(rho(0)) or something--not 
>> explicitly.
> If you end up rolling your own, it sounds like 'Method 
>20: Krylov
> space methods', in 'Nineteen Dubious Ways to Compute the 
> of a Matrix' is for you:
> http://www.cs.cornell.edu/cv/researchpdf/19ways+.pdf

You might be interested in



More information about the SciPy-User mailing list