Estoy interesado en calcular la solución de un sistema lage de EDO utilizando un método krylov como en [1]. Dicho método involucra funciones relacionadas con la exponencial (las llamadas funciones ). Básicamente consiste en calcular la acción de la función de matriz mediante la construcción de un subespacio de Krylov utilizando la iteración de Arnoldi y la proyección de la función en este subespacio. Esto reduce el problema para calcular el exponencial de una matriz de Hessenberg mucho más pequeña.
Soy consciente de que hay varios algoritmos para calcular el exponencial (ver [2] [3] y referencias en él). Me pregunto si hay un algoritmo especial para calcular el exponencial que puede aprovechar el hecho de que la matriz es Hessenberg.
[1] Sidje, RB (1998). Expokit: un paquete de software para calcular exponenciales de matriz. Transacciones ACM en software matemático (TOMS), 24 (1), 130-156.
[2] Moler, C. y Van Loan, C. (1978). Diecinueve formas dudosas de calcular el exponencial de una matriz. Revisión de SIAM, 20 (4), 801-836.
[3] Moler, C. y Van Loan, C. (2003). Diecinueve formas dudosas de calcular el exponencial de una matriz, veinticinco años después. Revisión de SIAM, 45 (1), 3-49.