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