DeepAI
Log In Sign Up

Pseudospectral roaming contour integral methods for convection-diffusion equations

12/13/2020
by   Nicola Guglielmi, et al.
0

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 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ópez-Fernández 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.

READ FULL TEXT
01/04/2022

Contour Integral Methods for time treatment in Reduced Basis Methods

In this paper we discuss a reduced basis method for linear evolution PDE...
10/18/2022

Analyses of the contour integral method for time fractional subdiffusion-normal transport equation

In this work, we theoretically and numerically discuss the time fraction...
10/20/2020

Flexible subspace iteration with moments for an effective contour integration-based eigensolver

Contour integration schemes are a valuable tool for the solution of diff...
07/31/2019

Domain Decomposition for the Closest Point Method

The discretization of elliptic PDEs leads to large coupled systems of eq...
08/08/2022

Numerical Evaluation of Mittag-Leffler Functions

The Mittag-Leffler function is computed via a quadrature approximation o...
06/07/2022

The fundamental solution of a 1D evolution equation with a sign changing diffusion coefficient

In this work we investigate a 1D evolution equation involving a divergen...
04/05/2018

Computing Stieltjes constants using complex integration

The Stieltjes constants γ_n are the coefficients appearing in the Lauren...