We generalize ideas in the recent literature and develop new ones in order to propose a general class of contour integral methods for linear convection-diffusion PDEs and in particular for those arising in finance. These methods aim to provide a numerical approximation of the solution by computing its inverse Laplace transform. The choice of the integration contour is determined by the computation of a few suitably weighted pseudo-spectral level sets of the leading operator of the equation. Parabolic and hyperbolic profiles proposed in the literature are investigated and compared to the elliptic contour originally proposed by Guglielmi, L\'opez-Fern\'andez and Nino. In summary, the article (i) provides a comparison among three different integration profiles; (ii) proposes a new fast pseudospectral roaming method; (iii) optimizes the selection of time windows on which one may arbitrarily approximate the solution by no extra computational cost with respect to the case of a fixed time instant; (iv) focuses extensively on computational aspects and it is the reference of the MATLAB code https://github.com/MattiaManucci/Contour_Integral_Methods.git, where all algorithms described here are implemented.
翻译:我们在最近的文献中推广了各种想法,并制定了新的想法,以便提出一个用于线性对流扩散PDE特别是金融中产生的PDE的等离整体法一般类别。这些方法的目的是通过计算其反Laplace变换,提供解决方案的数值近似值。集成轮廓的选择是通过计算该等式主要操作者几组适当加权的伪光谱级来决定的。对文献中提议的抛线和双曲剖面图进行了调查,并与古格利埃尔米、L\'opez-Fern\'andez和Nino最初提议的椭圆形轮廓作了比较。简而言之,文章(一)提供了三种不同的整合剖面的比较;(二)提出了一个新的快速假光谱漫游方法;(三)优化了时间窗口的选择,在固定时间瞬间的情况中,人们可以不增加计算成本,任意估计解决方案;(四)广泛侧重于计算方面,这是MATLAB代码 https://github.com/mattigusical_manucicius/tourcis