Abstract
Arbitrage-free prices u of European contracts on risky assets whose log-returns are modelled by Levy processes satisfy a parabolic partial integro-differential equation (PIDE) . This PIDE is localized to bounded domains and the error due to this localization is estimated. The localized PIDE is discretized by the θ -scheme in time and a wavelet Galerkin method with N degrees of freedom in log-price space. The dense matrix for can be replaced by a sparse matrix in the wavelet basis, and the linear systems in each implicit time step are solved approximatively with GMRES in linear complexity. The total work of the algorithm for M time steps is bounded by O(MN( log(N))2 ) operations and O(N log(N)) memory. The deterministic algorithm gives optimal convergence rates (up to logarithmic terms) for the computed solution in the same complexity as finite difference approximations of the standard Black–Scholes equation. Computational examples for various Levy price processes are presented.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callMore From: ESAIM: Mathematical Modelling and Numerical Analysis
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
