Cálculo sin matriz de combinación lineal de función phi por vector en integradores exponenciales
DOI:
https://doi.org/10.5281/zenodo.13915619Palabras clave:
funciones Phi, método numérico iterativo, método sin evaluación de matriz, integrador exponencial, problemas de valor inicial de grandes dimensiones, MSC 65F60, MSC 65F10, MSC 65L05Resumen
El cálculo numérico de combinaciones lineales de productos de funciones phi por vectores surge como un problema importante en la implementación de integradores exponenciales para Problemas de Valor Inicial (PVI) de grandes dimensiones. Esta comunicación corta trata con dicho cálculo cuando es inviable evaluar y almacenar las matrices Jacobianas involucradas. Para ello, se introduce una aproximación de Krylov-Padé apropiada y se determina su error. Se realizan simulaciones numéricas con matrices Jacobianas y campos vectoriales de PVI resultantes de la discretización espacial de ecuaciones de derivadas parciales conocidas.
Descargas
Citas
Arnoldi, W.E.: The principle of minimized iterations in the solution of the matrix eigenvalue problem. Quarterly of Applied Mathematics, 9:17–29, 1951. https://www.ams.org/journals/qam/1951-09-01/S0033-569X-1951-42792-9/S0033-569X-1951-42792-9.pdf.
Brown, P.N.: A local convergence theory for combined inexact-newton/finite-difference projection methods. SIAM Journal on Numerical Analysis, 24:407–434, 1987. https://epubs.siam.org/doi/abs/10.1137/0724031.
Carbonell, F., J.C. Jimenez, and L. Pedroso: Computing multiple integrals involving matrix exponentials. Journal of Computational and Applied Mathematics, 213:300–305, 2008. https://www.sciencedirect.com/science/article/pii/S0377042707000283.
Cruz, H. de la, R.J. Biscay, F. Carbonell, T. Ozaki, and J.C. Jimenez: A higher order local linearization method for solving ordinary differential equations. Applied Mathematics and Computation, 185:197–212, 2007. https://www.sciencedirect.com/science/article/abs/pii/S0096300306008514.
Hochbruck, M. and A. Ostermann: Exponential multistep methods of Adams-type. BIT Numerical Mathematics, 51:889–908, 2011.
https://link.springer.com/article/10.1007/s10543-011-0332-6.
Jimenez, J.C., H. de la Cruz, and P. De Maio: Efficient computation of phi-functions in exponential integrators. Journal of Computational and Applied Mathematics, 374:112758, 2020. https://www.sciencedirect.com/science/article/abs/pii/S0377042720300492.
Jimenez, J.C., L. Pedroso, F. Carbonell, and V. Hernadez: Local linearization method for numerical integration of
delay differential equations. SIAM Journal of Numerical Analysis, 44:2584–2609, 2006. https://epubs.siam.org/doi/abs/10.1137/040607356.
Knoll, D.A. and Keyes D.E.: Jacobian-free Newton-Krylov methods: a survey of approaches and applications. Journal of Computational
Physics, 193:357–397, 2004. https://www.sciencedirect.com/science/article/abs/pii/S0021999103004340.
Naranjo-Noda, F.S. and J.C. Jimenez: Locally linearized runge-kutta method of dormand and prince for large systems of initial value problems. Journal of Computational Physics, 426:109946, 2021. https://www.sciencedirect.com/science/article/abs/pii/S0021999120307208.
Naranjo-Noda, F.S. and J.C. Jimenez: Computing high dimensional multiple integrals involving matrix exponentials. Journal of Computational and Applied Mathematics, 421:114844, 2023. https://www.sciencedirect.com/science/article/abs/pii/S0377042722004423.
Naranjo-Noda, F.S. and J.C. Jimenez: Jacobian-free high order local linearization methods for large systems of initial value problems. Applied Numerical Mathematics, 187:158–175, 2023. https://www.sciencedirect.com/science/article/abs/pii/S0168927423000405.
Niesen, J. and Wright WM. Algorithm 919: Algorithm 919: A Krylov subspace algorithm for evaluating the ϕ-functions appearing in exponential integrators. ACM Transactions on Mathematical Software (TOMS), 38:1–9, 2012. https://dl.acm.org/doi/abs/10.1145/2168773.2168781.
Sidje, R.B.: Expokit: A software package for computing matrix exponentials. ACM Transactions on Mathematical Software (TOMS), 24:130–156, 1998. https://dl.acm.org/doi/abs/10.1145/285861.285868.
Tokman, M.: Efficient integration of large stiff systems of ODEs with exponential propagation iterative (EPI) methods. Journal of Computational Physics, 213(2):748–776, 2006. https://www.sciencedirect.com/science/article/abs/pii/S0021999105004158.
Publicado
Versiones
Cómo citar
Número
Sección
Licencia
Derechos de autor 2024 Frank Sadan Naranjo-Noda, Juan Carlos Jiménez Sobrino
Esta obra está bajo una licencia internacional Creative Commons Atribución 4.0.
Esta licencia permite copiar y redistribuir el material en cualquier medio o formato bajo los siguientes términos: se debe dar crédito de manera adecuada, no se puede hacer uso del material con propósitos comerciales, y si remezcla, transforma o crea a partir del material, no podrá distribuir el material modificado. Bajo la licencia mencionada, los autores mantienen los derechos de autor de su trabajo.
