Numerical inverse Laplace transform for convection-diffusion equations

Guglielmi, Nicola; López-Fernández, María; Nino, Giancarlo

doi:10.1090/mcom/3497

Numerical inverse Laplace transform for convection-diffusion equations
HTML articles powered by AMS MathViewer

by Nicola Guglielmi, María López-Fernández and Giancarlo Nino HTML | PDF

Math. Comp. 89 (2020), 1161-1191 Request permission

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.

References

L. Banjai, M. López-Fernández, and A. Schädle, Fast and oblivious algorithms for dissipative and two-dimensional wave equations, SIAM J. Numer. Anal. 55 (2017), no. 2, 621–639. MR 3623208, DOI 10.1137/16M1070657
Fischer Black and Myron Scholes, The pricing of options and corporate liabilities, J. Polit. Econ. 81 (1973), no. 3, 637–654. MR 3363443, DOI 10.1086/260062
Benedict Dingfelder and J. A. C. Weideman, An improved Talbot method for numerical Laplace transform inversion, Numer. Algorithms 68 (2015), no. 1, 167–183. MR 3296705, DOI 10.1007/s11075-014-9895-z
Ivan P. Gavrilyuk and Vladimir L. Makarov, Exponentially convergent parallel discretization methods for the first order evolution equations, Comput. Methods Appl. Math. 1 (2001), no. 4, 333–355. MR 1892950, DOI 10.2478/cmam-2001-0022
Steven L. Heston, A closed-form solution for options with stochastic volatility with applications to bond and currency options, Rev. Financ. Stud. 6 (1993), no. 2, 327–343. MR 3929676, DOI 10.1093/rfs/6.2.327
J. C. Hull, Options, Futures and Other Derivatives, 6th ed., Prentice Hall, New Jersey, 2006.
K. J. In ’t Hout and S. Foulon, ADI finite difference schemes for option pricing in the Heston model with correlation, Int. J. Numer. Anal. Model. 7 (2010), no. 2, 303–320. MR 2587424
K. J. in ’t Hout, ADI Schemes in the Numerical Solution of the Heston PDE, Numerical Analysis and Applied Mathematics, eds. T. E. Simos et.al., AIP Conf. Proc: 936, 2007.
K. J. in ’t Hout and B. D. Welfert, Stability of ADI schemes applied to convection-diffusion equations with mixed derivative terms, Appl. Numer. Math. 57 (2007), no. 1, 19–35. MR 2279504, DOI 10.1016/j.apnum.2005.11.011
K. J. in ’t Hout and B. D. Welfert, Unconditional stability of second-order ADI schemes applied to multi-dimensional diffusion equations with mixed derivative terms, Appl. Numer. Math. 59 (2009), no. 3-4, 677–692. MR 2492284, DOI 10.1016/j.apnum.2008.03.016
K. J. in ’t Hout and J. A. C. Weideman, A contour integral method for the Black-Scholes and Heston equations, SIAM J. Sci. Comput. 33 (2011), no. 2, 763–785. MR 2785971, DOI 10.1137/090776081
Mohsin Javed and Lloyd N. Trefethen, A trapezoidal rule error bound unifying the Euler-Maclaurin formula and geometric convergence for periodic functions, Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 470 (2014), no. 2161, 20130571, 9. MR 3145072, DOI 10.1098/rspa.2013.0571
María López-Fernández, Christian Lubich, Cesar Palencia, and Achim Schädle, Fast Runge-Kutta approximation of inhomogeneous parabolic equations, Numer. Math. 102 (2005), no. 2, 277–291. MR 2206466, DOI 10.1007/s00211-005-0624-3
M. López-Fernández and C. Palencia, On the numerical inversion of the Laplace transform of certain holomorphic mappings, Appl. Numer. Math. 51 (2004), no. 2-3, 289–303. MR 2091405, DOI 10.1016/j.apnum.2004.06.015
María López-Fernández, César Palencia, and Achim Schädle, A spectral order method for inverting sectorial Laplace transforms, SIAM J. Numer. Anal. 44 (2006), no. 3, 1332–1350. MR 2231867, DOI 10.1137/050629653
G. Nino, Laplace transform methods based on pseudospectral roaming for convection-diffusion equations, PhD thesis, Gran Sasso Science Institute, L’Aquila, 2017.
Satish C. Reddy and Lloyd N. Trefethen, Pseudospectra of the convection-diffusion operator, SIAM J. Appl. Math. 54 (1994), no. 6, 1634–1649. MR 1301275, DOI 10.1137/S0036139993246982
Dongwoo Sheen, Ian H. Sloan, and Vidar Thomée, A parallel method for time discretization of parabolic equations based on Laplace transformation and quadrature, IMA J. Numer. Anal. 23 (2003), no. 2, 269–299. MR 1975267, DOI 10.1093/imanum/23.2.269
A. Talbot, The accurate numerical inversion of Laplace transforms, J. Inst. Math. Appl. 23 (1979), no. 1, 97–120. MR 526286
Lloyd N. Trefethen and Mark Embree, Spectra and pseudospectra, Princeton University Press, Princeton, NJ, 2005. The behavior of nonnormal matrices and operators. MR 2155029
Lloyd N. Trefethen and J. A. C. Weideman, The exponentially convergent trapezoidal rule, SIAM Rev. 56 (2014), no. 3, 385–458. MR 3245858, DOI 10.1137/130932132
J. A. C. Weideman, Improved contour integral methods for parabolic PDEs, IMA J. Numer. Anal. 30 (2010), no. 1, 334–350. MR 2580562, DOI 10.1093/imanum/drn074
Thomas G. Wright, EigTool, http://www.comlab.ox.ac.uk/pseudospectra/eigtool/, 2002.