Researcher profile

Jean-François Chassagneux

Jean-François Chassagneux contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
15works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Deep Runge-Kutta schemes for BSDEs

We propose a new probabilistic scheme which combines deep learning techniques with high order schemes for backward stochastic differential equations belonging to the class of Runge-Kutta methods to solve high-dimensional semi-linear parabolic partial differential equations. Our approach notably extends the one introduced in [Hure Pham Warin 2020] for the implicit Euler scheme to schemes which are more efficient in terms of discrete-time error. We establish some convergence results for our implemented schemes under classical regularity assumptions. We also illustrate the efficiency of our method for different schemes of order one, two and three. Our numerical results indicate that the Crank-Nicolson schemes is a good compromise in terms of precision, computational cost and numerical implementation.

preprint2021arXiv

A learning scheme by sparse grids and Picard approximations for semilinear parabolic PDEs

Relying on the classical connection between Backward Stochastic Differential Equations (BSDEs) and non-linear parabolic partial differential equations (PDEs), we propose a new probabilistic learning scheme for solving high-dimensional semi-linear parabolic PDEs. This scheme is inspired by the approach coming from machine learning and developed using deep neural networks in Han and al. [32]. Our algorithm is based on a Picard iteration scheme in which a sequence of linear-quadratic optimisation problem is solved by means of stochastic gradient descent (SGD) algorithm. In the framework of a linear specification of the approximation space, we manage to prove a convergence result for our scheme, under some smallness condition. In practice, in order to be able to treat high-dimensional examples, we employ sparse grid approximation spaces. In the case of periodic coefficients and using pre-wavelet basis functions, we obtain an upper bound on the global complexity of our method. It shows in particular that the curse of dimensionality is tamed in the sense that in order to achieve a root mean squared error of order $ε$, for a prescribed precision $ε$, the complexity of the Picard algorithm grows polynomially in $ε^{-1}$ up to some logarithmic factor $ |log(ε)| $ which grows linearly with respect to the PDE dimension. Various numerical results are presented to validate the performance of our method and to compare them with some recent machine learning schemes proposed in Han and al. [20] and Huré and al. [37].

preprint2021arXiv

Numerical approximation of singular Forward-Backward SDEs

In this work, we study the numerical approximation of a class of singular fully coupled forward backward stochastic differential equations. These equations have a degenerate forward component and non-smooth terminal condition. They are used, for example, in the modeling of carbon market[9] and are linked to scalar conservation law perturbed by a diffusion. Classical FBSDEs methods fail to capture the correct entropy solution to the associated quasi-linear PDE. We introduce a splitting approach that circumvent this difficulty by treating differently the numerical approximation of the diffusion part and the non-linear transport part. Under the structural condition guaranteeing the well-posedness of the singular FBSDEs [8], we show that the splitting method is convergent with a rate $1/2$. We implement the splitting scheme combining non-linear regression based on deep neural networks and conservative finite difference schemes. The numerical tests show very good results in possibly high dimensional framework.

preprint2021arXiv

Reflected BSDEs in non-convex domains

This paper establishes the well-posedness of reflected backward stochastic differential equations in the non-convex domains that satisfy a weaker version of the star-shaped property. The main results are established (i) in a Markovian framework with Hölder-continuous generator and terminal condition and (ii) in a general setting under a smallness assumption on the input data. We also investigate the connections between this well-posedness result and the theory of martingales on manifolds.

preprint2020arXiv

Switching problems with controlled randomisation and associated obliquely reflected BSDEs

We introduce and study a new class of optimal switching problems, namely switching problem with controlled randomisation, where some extra-randomness impacts the choice of switching modes and associated costs. We show that the optimal value of the switching problem is related to a new class of multidimensional obliquely reflected BSDEs. These BSDEs allow as well to construct an optimal strategy and thus to solve completely the initial problem. The other main contribution of our work is to prove new existence and uniqueness results for these obliquely reflected BSDEs. This is achieved by a careful study of the domain of reflection and the construction of an appropriate oblique reflection operator in order to invoke results from [7].

preprint2019arXiv

A numerical scheme for the quantile hedging problem

We consider the numerical approximation of the quantile hedging price in a non-linear market. In a Markovian framework, we propose a numerical method based on a Piecewise Constant Policy Timestepping (PCPT) scheme coupled with a monotone finite difference approximation. We prove the convergence of our algorithm combining BSDE arguments with the Barles & Jakobsen and Barles & Souganidis approaches for non-linear equations. In a numerical section, we illustrate the efficiency of our scheme by considering a financial example in a market with imperfections.

preprint2019arXiv

Weak quantitative propagation of chaos via differential calculus on the space of measures

Consider the metric space $(\mathcal{P}_2(\mathbb{R}^d),W_2)$ of square integrable laws on $\mathbb{R}^d$ with the topology induced by the 2-Wasserstein distance $W_2$. Let $Φ: \mathcal{P}_2( \mathbb{R}^d) \to \mathbb{R}$ be a function and $μ_N$ be the empirical measure of a sample of $N$ random variables distributed as $μ$. The main result of this paper is to show that under suitable regularity conditions, we have \[ |Φ(μ) - \mathbb{E}Φ(μ_N)|= \sum_{j=1}^{k-1}\frac{C_j}{N^j} + O(\frac{1}{N^k}), \] for some positive constants $C_1, \ldots, C_{k-1}$ that do not depend on $N$, where $k$ corresponds to the degree of smoothness. We distinguish two cases: a) $μ_N$ is the empirical measure of $N$-samples from $μ$; b) $μ$ is a marginal law of McKean-Vlasov stochastic differential equation in which case $μ_N$ is an empirical law of marginal laws of the corresponding particle system. The first case is studied using functional derivatives on the space of measures. The second case relies on an Itô-type formula for the flow of probability measures and is intimately connected to PDEs on the space of measures, called the master equation in the literature of mean-field games. We state the general regularity conditions required for each case and analyse the regularity in the case of functionals of the laws of McKean-Vlasov SDEs. Ultimately, this work reveals quantitative estimates of propagation of chaos for interacting particle systems. Furthermore, we are able to provide weak propagation of chaos estimates for ensembles of interacting particles and show that these may have some remarkable properties.

preprint2018arXiv

A sparse grid approach to balance sheet risk measurement

In this work, we present a numerical method based on a sparse grid approximation to compute the loss distribution of the balance sheet of a financial or an insurance company. We first describe, in a stylised way, the assets and liabilities dynamics that are used for the numerical estimation of the balance sheet distribution. For the pricing and hedging model, we chose a classical Black & Scholes model with a stochastic interest rate following a Hull & White model. The risk management model describing the evolution of the parameters of the pricing and hedging model is a Gaussian model. The new numerical method is compared with the traditional nested simulation approach. We review the convergence of both methods to estimate the risk indicators under consideration. Finally, we provide numerical results showing that the sparse grid approach is extremely competitive for models with moderate dimension.

preprint2016arXiv

A backward dual representation for the quantile hedging of Bermudan options

Within a Markovian complete financial market, we consider the problem of hedging a Bermudan option with a given probability. Using stochastic target and duality arguments, we derive a backward numerical scheme for the Fenchel transform of the pricing function. This algorithm is similar to the usual American backward induction, except that it requires two additional Fenchel transformations at each exercise date. We provide numerical illustrations.

preprint2016arXiv

Numerical simulation of quadratic BSDEs

This article deals with the numerical approximation of Markovian backward stochastic differential equations (BSDEs) with generators of quadratic growth with respect to $z$ and bounded terminal conditions. We first study a slight modification of the classical dynamic programming equation arising from the time-discretization of BSDEs. By using a linearization argument and BMO martingales tools, we obtain a comparison theorem, a priori estimates and stability results for the solution of this scheme. Then we provide a control on the time-discretization error of order $\frac{1}{2}-\varepsilon$ for all $\varepsilon>0$. In the last part, we give a fully implementable algorithm for quadratic BSDEs based on quantization and illustrate our convergence results with numerical examples.

preprint2015arXiv

A Probabilistic approach to classical solutions of the master equation for large population equilibria

We analyze a class of nonlinear partial differential equations (PDEs) defined on $\mathbb{R}^d \times \mathcal{P}_2(\mathbb{R}^d),$ where $\mathcal{P}_2(\mathbb{R}^d)$ is the Wasserstein space of probability measures on $\mathbb{R}^d$ with a finite second-order moment. We show that such equations admit a classical solutions for sufficiently small time intervals. Under additional constraints, we prove that their solution can be extended to arbitrary large intervals. These nonlinear PDEs arise in the recent developments in the theory of large population stochastic control. More precisely they are the so-called master equations corresponding to asymptotic equilibria for a large population of controlled players with mean-field interaction and subject to minimization constraints. The results in the paper are deduced by exploiting this connection. In particular, we study the differentiability with respect to the initial condition of the flow generated by a forward-backward stochastic system of McKean-Vlasov type. As a byproduct, we prove that the decoupling field generated by the forward-backward system is a classical solution of the corresponding master equation. Finally, we give several applications to mean-field games and to the control of McKean-Vlasov diffusion processes.

preprint2014arXiv

Numerical stability analysis of the Euler scheme for BSDEs

In this paper, we study the qualitative behaviour of approximation schemes for Backward Stochastic Differential Equations (BSDEs) by introducing a new notion of numerical stability. For the Euler scheme, we provide sufficient conditions in the one-dimensional and multidimensional case to guarantee the numerical stability. We then perform a classical Von Neumann stability analysis in the case of a linear driver $f$ and exhibit necessary conditions to get stability in this case. Finally, we illustrate our results with numerical applications.

preprint2014arXiv

Runge-Kutta schemes for backward stochastic differential equations

We study the convergence of a class of Runge-Kutta type schemes for backward stochastic differential equations (BSDEs) in a Markovian framework. The schemes belonging to the class under consideration benefit from a certain stability property. As a consequence, the overall rate of the convergence of these schemes is controlled by their local truncation error. The schemes are categorized by the number of intermediate stages implemented between consecutive partition time instances. We show that the order of the schemes matches the number $p$ of intermediate stages for $p\le3$. Moreover, we show that the so-called order barrier occurs at $p=3$, that is, that it is not possible to construct schemes of order $p$ with $p$ stages, when $p>3$. The analysis is done under sufficient regularity on the final condition and on the coefficients of the BSDE.

preprint2013arXiv

Linear multi-step schemes for BSDEs

We study the convergence rate of a class of linear multi-step methods for BSDEs. We show that, under a sufficient condition on the coefficients, the schemes enjoy a fundamental stability property. Coupling this result to an analysis of the truncation error allows us to design approximation with arbitrary order of convergence. Contrary to the analysis performed in \cite{zhazha10}, we consider general diffusion model and BSDEs with driver depending on $z$. The class of methods we consider contains well known methods from the ODE framework as Nystrom, Milne or Adams methods. We also study a class of Predictor-Correctot methods based on Adams methods. Finally, we provide a numerical illustration of the convergence of some methods.

preprint2013arXiv

When terminal facelift enforces Delta constraints

This paper deals with the super-replication of non path-dependent European claims under additional convex constraints on the number of shares held in the portfolio. The corresponding super-replication price of a given claim has been widely studied in the literature and its terminal value, which dominates the claim of interest, is the so-called facelift transform of the claim. We investigate under which conditions the super-replication price and strategy of a large class of claims coincide with the exact replication price and strategy of the facelift transform of this claim. In one dimension, we observe that this property is satisfied for any local volatility model. In any dimension, we exhibit an analytical necessary and sufficient condition for this property, which combines the dynamics of the stock together with the characteristics of the closed convex set of constraints. To obtain this condition, we introduce the notion of first order viability property for linear parabolic PDEs. We investigate in details several practical cases of interest: multidimensional Black Scholes model, non-tradable assets or short selling restrictions.