In this paper, we present a variable step size implementation of exponential Rosenbrock-type methods of orders 2, 3 and 4. These integrators require the evaluation of exponential and related functions of the Jacobian matrix. To this aim, the Real Leja Points Method is used. It is shown that the properties of this method combine well with the particular requirements of Rosenbrock-type integrators. We verify our implementation with some numerical experiments in MATLAB, where we solve semilinear parabolic PDEs in one and two space dimensions. We further present some numerical experiments in FORTRAN. where we compare out-method with other methods from literature. We find a great potential Of Our method for non-normal matrices. Such matrices typically arise in parabolic problems with large advection in combination with moderate diffusion and mildly stiff reactions.
Implementation of exponential Rosenbrock-type methods / Caliari, Marco; A., Ostermann. - In: APPLIED NUMERICAL MATHEMATICS. - ISSN 0168-9274. - STAMPA. - 59:3-4(2009), pp. 568-581.
Titolo: | Implementation of exponential Rosenbrock-type methods |
Autori: | |
Data di pubblicazione: | 2009 |
Rivista: | |
Citazione: | Implementation of exponential Rosenbrock-type methods / Caliari, Marco; A., Ostermann. - In: APPLIED NUMERICAL MATHEMATICS. - ISSN 0168-9274. - STAMPA. - 59:3-4(2009), pp. 568-581. |
Handle: | http://hdl.handle.net/11562/320279 |
Appare nelle tipologie: | 01.01 Articolo in Rivista |
File in questo prodotto:
File | Descrizione | Tipologia | Licenza | |
---|---|---|---|---|
preCO09.pdf | Documento in Post-print | ![]() | Open Access Visualizza/Apri | |
CO09.pdf | Versione dell'editore | Accesso ristretto | Administrator Richiedi una copia |