README.rst

rpylib

Scope

This Python pricing library was developed as a companion code for the paper:

"A Weak MLMC Scheme for Levy-copula-driven SDEs with Applications to the Pricing of Credit, Equity and Interest Rate Derivatives".

The results consist of the numerical analysis of:

• the benchmark of the Continuous-Time Markov Chain (CTMC) scheme approximation against the series representation [1]
• the benchmark of the CTMC scheme against the closed-form formula for First-to-Default CDS [2]
• the weak and strong convergence of the multilevel CTMC scheme as well as the convergence rate of the cost w.r.t the rmse compared to the standard Monte-Carlo; these results mimic those of Giles [3] for diffusion processes

The different convergence rates considered in our case are dependent on the Blumenthal-Getoor index of the underlying Levy process.

Results

The main results are presented in the form of 4 graphs (as in Giles [3] [4]) |

• log2(vl), the log level variances in function of the level l
• log2|ml| the log level means in function of the level l
• Nl (optimal number of Monte-Carlo paths for the level l) in function of the level l
• the total costs of the multilevel Monte-Carlo and the standard Monte-Carlo in function of the rmse (root-mean square error)

MLMC applied to CGMY with beta=1.5

Scripts

For the paper:

See the slurm folder.

Other scripts:

Other scripts are available in rpylib/scripts/statistics. These scripts allow to plot the distribution of the spot underlying of the Levy process simulated by Monte-Carlo (either directly from the SDE or from the CTMC scheme).

---

Contact:

Any feedback on this project will be appreciated, please log a new Issue or email me.

[1]	Levy Copulas: Review of Recent Results, P. Tankov

[2]	A Structural Jump Threshold Framework for Credit Risk, P. Garreau, A. Kercheval

[3]	(1, 2) Multilevel Monte Carlo Path Simulation, M.B. Giles

[4]	Multilevel Monte Carlo methods, M.B. Giles