In this paper a novel contour integral method is proposed for linear convection-diffusion equations. The method is based on the inversion of the Laplace transform and makes use of a contour given by an elliptic arc joined symmetrically to two half-lines. The trapezoidal rule is the chosen integration method for the numerical inversion of the Laplace transform, due to its well-known fast convergence properties when applied to analytic functions. Error estimates are provided as well as careful indications about the choice of several involved parameters. The method selects the elliptic arc in the integration contour by an algorithmic strategy based on the computation of pseudospectral level sets of the discretized differential operator. In this sense the method is general and can be applied to any linear convection-diffusion equation without knowing any a priori information about its pseudospectral geometry. Numerical experiments performed on the Black-Scholes (1D) and Heston (2D) equations show that the method is competitive with other contour integral methods available in the literature.

Numerical inverse laplace transform for convection-diffusion equations

Guglielmi N.;
2020

Abstract

In this paper a novel contour integral method is proposed for linear convection-diffusion equations. The method is based on the inversion of the Laplace transform and makes use of a contour given by an elliptic arc joined symmetrically to two half-lines. The trapezoidal rule is the chosen integration method for the numerical inversion of the Laplace transform, due to its well-known fast convergence properties when applied to analytic functions. Error estimates are provided as well as careful indications about the choice of several involved parameters. The method selects the elliptic arc in the integration contour by an algorithmic strategy based on the computation of pseudospectral level sets of the discretized differential operator. In this sense the method is general and can be applied to any linear convection-diffusion equation without knowing any a priori information about its pseudospectral geometry. Numerical experiments performed on the Black-Scholes (1D) and Heston (2D) equations show that the method is competitive with other contour integral methods available in the literature.
Contour integral methods; Laplace transform; Numerical inversion of laplace transform; Pseudospectra; Quadrature for analytic functions; Trapezoidal rule
File in questo prodotto:
File Dimensione Formato  
2020_MathComput_89_Guglielmi.pdf

non disponibili

Tipologia: Altro materiale allegato
Licenza: Non pubblico
Dimensione 559.98 kB
Formato Adobe PDF
559.98 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/20.500.12571/710
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 4
  • ???jsp.display-item.citation.isi??? ND
social impact