[Date Index] [Thread Index] [Author Index]

Jacobi Matrix Exponential

• To: mathgroup at smc.vnet.net
• Subject: [mg39661] Jacobi Matrix Exponential
• From: Kyriakos Chourdakis <tuxedomoon at Yahoo.com>
• Date: Sat, 1 Mar 2003 02:47:06 -0500 (EST)
• Reply-to: k.chourdakis at qmul.ac.uk
• Sender: owner-wri-mathgroup at wolfram.com

Dear all,

I have to compute matrix exponentials of complex
matrices. The matrices are quite large, e.g. up to
1000x1000, but only the diagonal and the adjacent
elements are non zero. The first and last rows are
zero. The matrix is therefore of the form [x denotes a
non-zero complex element, not all x's are the same]:

0 0 0 0 0 0 0 .....
x x x 0 0 0 0 .....
0 x x x 0 0 0 .....
0 0 x x x 0 0 .....
0 0 0 x x x 0 .....
. . . . . . .
. . . . . . .

Does anyone have in mind, or has experimented with, a
quick way of computing such exponentials. The built-in
function is very slow.

My experiments with Pade approximations include matrix
inversion and are therefore slow as well. Perhaps a
quick way of inverting such matrices would also be
helpful.

Best,

Kyriakos


__________________________________________________
Do You Yahoo!?
Everything you'll ever need on one web page
from News and Sport to Email and Music Charts
http://uk.my.yahoo.com