(I think John probably needs to be CC'd, since I don't know if he's subscribed to this list; so I'm CCing him this message. John, Robert Kern's query is below.) Robert Kern wrote:
On Fri, Feb 27, 2009 at 15:52, <jason-sage@creativetrax.com> wrote:
John Cremona posted the following message to the sage development list about matrix exponentials. I'm copying it to here since it asks about the scipy matrix exponential method (we say numpy below, but we really mean scipy...)
John Cremona wrote:
I have just been to a colloquium talk by numerical analyst Nick Higham (Manchester) called "How to compute and not to compute a matrix exponential". He has new methods which are now in mathematica, matlab and NAG but (apparantly) nowhere else.
Are the good methods in this paper?
A New Scaling and Squaring Algorithm for the Matrix Exponential (with Awad Al-Mohy), MIMS EPrint 2009.9, January 2009. [new]
http://eprints.ma.man.ac.uk/1217/01/covered/MIMS_ep2009_9.pdf