Source author record

Gilles Pagès

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.PR q-fin.CP math.ST q-fin.MF q-fin.PR q-fin.TR Statistics Theory Computational Engineering, Finance, and Science Machine Learning math.OC

Catalog footprint

Research overview

34works

10topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2026arXiv

Strong Solutions and Quantization-Based Numerical Schemes for a Class of Non-Markovian Volatility Models

We investigate a class of non-Markovian processes that hold particular relevance in the realm of mathematical finance. This family encompasses path-dependent volatility models, including those pioneered by [Platen and Rendek, 2018] and, more recently, by [Guyon and Lekeufack, 2023]. Our study unfolds in two principal phases. In the first phase, we introduce a functional quantization scheme based on an extended version of the Lamperti transformation that we propose to handle the presence of a memory term incorporated into the diffusion coefficient. In the second phase, we study the problem of existence and uniqueness of a strong solution for the SDEs related to the examples that motivate our study, in order to provide a theoretical basis to correctly apply the proposed numerical schemes.

preprint2023arXiv

Langevin algorithms for Markovian Neural Networks and Deep Stochastic control

Stochastic Gradient Descent Langevin Dynamics (SGLD) algorithms, which add noise to the classic gradient descent, are known to improve the training of neural networks in some cases where the neural network is very deep. In this paper we study the possibilities of training acceleration for the numerical resolution of stochastic control problems through gradient descent, where the control is parametrized by a neural network. If the control is applied at many discretization times then solving the stochastic control problem reduces to minimizing the loss of a very deep neural network. We numerically show that Langevin algorithms improve the training on various stochastic control problems like hedging and resource management, and for different choices of gradient descent methods.

preprint2022arXiv

Convergence of Langevin-Simulated Annealing algorithms with multiplicative noise

We study the convergence of Langevin-Simulated Annealing type algorithms with multiplicative noise, i.e. for $V : \mathbb{R}^d \to \mathbb{R}$ a potential function to minimize, we consider the stochastic equation $dY_t = - σσ^\top \nabla V(Y_t) dt + a(t)σ(Y_t)dW_t + a(t)^2Υ(Y_t)dt$ , where $(W_t)$ is a Brownian motion, where $σ: \mathbb{R}^d \to \mathcal{M}_d(\mathbb{R})$ is an adaptive (multiplicative) noise, where $a : \mathbb{R}^+ \to \mathbb{R}^+$ is a function decreasing to $0$ and where $Υ$ is a correction term. This setting can be applied to optimization problems arising in Machine Learning. The case where $σ$ is a constant matrix has been extensively studied however little attention has been paid to the general case. We prove the convergence for the $L^1$ -Wasserstein distance of $Y_t$ and of the associated Euler-scheme $\bar{Y}_t$ to some measure $ν^\star$ which is supported by $\text{argmin}(V)$ and give rates of convergence to the instantaneous Gibbs measure $ν_{a(t)}$ of density $\propto \exp(-2V(x)/a(t)^2)$ . To do so, we first consider the case where $a$ is a piecewise constant function. We find again the classical schedule $a(t) = A\log^{-1/2}(t)$ . We then prove the convergence for the general case by giving bounds for the Wasserstein distance to the stepwise constant case using ergodicity properties.

preprint2022arXiv

Convergence of Langevin-Simulated Annealing algorithms with multiplicative noise II: Total Variation

We study the convergence of Langevin-Simulated Annealing type algorithms with multiplicative noise, i.e. for $V : \mathbb{R}^d \to \mathbb{R}$ a potential function to minimize, we consider the stochastic differential equation $dY_t = - σσ^\top \nabla V(Y_t) dt + a(t)σ(Y_t)dW_t + a(t)^2Υ(Y_t)dt$ , where $(W_t)$ is a Brownian motion, where $σ: \mathbb{R}^d \to \mathcal{M}_d(\mathbb{R})$ is an adaptive (multiplicative) noise, where $a : \mathbb{R}^+ \to \mathbb{R}^+$ is a function decreasing to $0$ and where $Υ$ is a correction term. Allowing $σ$ to depend on the position brings faster convergence in comparison with the classical Langevin equation $dY_t = -\nabla V(Y_t)dt + σdW_t$ . In a previous paper we established the convergence in $L^1$ -Wasserstein distance of $Y_t$ and of its associated Euler scheme $\bar{Y}_t$ to $\text{argmin}(V)$ with the classical schedule $a(t) = A\log^{-1/2}(t)$ . In the present paper we prove the convergence in total variation distance. The total variation case appears more demanding to deal with and requires regularization lemmas.

preprint2022arXiv

Functional convex order for the scaled McKean-Vlasov processes

We establish the functional convex order results for two scaled McKean-Vlasov processes $X=(X_{t})_{t\in[0, T]}$ and $Y=(Y_{t})_{t\in[0, T]}$ defined on a filtered probability space $(Ω, \mathcal{F}, (\mathcal{F}_{t})_{t\geq0}, \mathbb{P})$ by $\begin{cases} dX_{t}= b(t, X_{t}, μ_{t})dt+σ(t, X_{t}, μ_{t})dB_{t}, \;\;X_{0}\in L^{p}(\mathbb{P}),\\ dY_{t}\,= b(t, \,Y_{t}\,,\, ν_{t})dt+θ(t, \,Y_{t}\,,\, ν_{t})dB_{t}, \;\;Y_{0}\in L^{p}(\mathbb{P}), \end{cases}$ where $p\geq2$ , for every $t\in[0, T]$ , $μ_t$ , $ν_t$ denote the probability distribution of $X_t$ , $Y_t$ respectively and the drift coefficient $b(t, x, μ)$ is affine in $x$ (scaled). If we make the convexity and monotony assumption (only) on $σ$ and if $σ\preceqθ$ with respect to the partial matrix order, the convex order for the initial random variable $X_0 \preceq_{\,cv} Y_0$ can be propagated to the whole path of process $X$ and $Y$ . That is, if we consider a convex functional $F$ defined on the path space with polynomial growth, we have $\mathbb{E}F(X)\leq\mathbb{E}F(Y)$ ; for a convex functional $G$ defined on the product space involving the path space and its marginal distribution space, we have $\mathbb{E}\,G\big(X, (μ_t)_{t\in[0, T]}\big)\leq \mathbb{E}\,G\big(Y, (ν_t)_{t\in[0, T]}\big)$ under appropriate conditions. The symmetric setting is also valid, that is, if $θ\preceq σ$ and $Y_0 \leq X_0$ with respect to the convex order, then $\mathbb{E}\,F(Y) \leq \mathbb{E}\,F(X)$ and $\mathbb{E}\,G\big(Y, (ν_t)_{t\in[0, T]}\big)\leq \mathbb{E}\,G(X, (μ_t)_{t\in[0, T]})$ . The proof is based on several forward and backward dynamic programming principles and the convergence of the Euler scheme of the McKean-Vlasov equation.

preprint2020arXiv

Convergence rate of optimal quantization grids and application to empirical measure

We study the convergence rate of the optimal quantization for a probability measure sequence $(μ_{n})_{n\in\mathbb{N}^{*}}$ on $\mathbb{R}^{d}$ converging in the Wasserstein distance in two aspects: the first one is the convergence rate of optimal quantizer $x^{(n)}\in(\mathbb{R}^{d})^{K}$ of $μ_{n}$ at level $K$ ; the other one is the convergence rate of the distortion function valued at $x^{(n)}$ , called the "performance" of $x^{(n)}$ . Moreover, we also study the mean performance of the optimal quantization for the empirical measure of a distribution $μ$ with finite second moment but possibly unbounded support. As an application, we show that the mean performance for the empirical measure of the multidimensional normal distribution $\mathcal{N}(m, Σ)$ and of distributions with hyper-exponential tails behave like $\mathcal{O}(\frac{\log n}{\sqrt{n}})$ . This extends the results from [BDL08] obtained for compactly supported distribution. We also derive an upper bound which is sharper in the quantization level $K$ but suboptimal in $n$ by applying results in [FG15].

preprint2020arXiv

New approach to greedy vector quantization

We extend some rate of convergence results of greedy quantization sequences already investigated in arXiv:1409.0732 [math.PR]. We show, for a more general class of distributions satisfying a certain control, that the quantization error of these sequences have an $n^{-\frac1d}$ rate of convergence and that the distortion mismatch property is satisfied. We will give some non-asymptotic Pierce type estimates. The recursive character of greedy vector quantization allows some improvements to the algorithm of computation of these sequences and the implementation of a recursive formula to quantization-based numerical integration. Furthermore, we establish further properties of sub-optimality of greedy quantization sequences.

preprint2020arXiv

New Weak Error bounds and expansions for Optimal Quantization

We propose new weak error bounds and expansion in dimension one for optimal quantization-based cubature formula for different classes of functions, such that piecewise affine functions, Lipschitz convex functions or differentiable function with piecewise-defined locally Lipschitz or $α$ -Hölder derivatives. This new results rest on the local behaviors of optimal quantizers, the $L^r$ - $L^s$ distribution mismatch problem and Zador's Theorem. This new expansion supports the definition of a Richardson-Romberg extrapolation yielding a better rate of convergence for the cubature formula. An extension of this expansion is then proposed in higher dimension for the first time. We then propose a novel variance reduction method for Monte Carlo estimators, based on one dimensional optimal quantizers.

preprint2020arXiv

Quantization-based Bermudan option pricing in the $FX$ world

This paper proposes two numerical solution based on Product Optimal Quantization for the pricing of Foreign Echange (FX) linked long term Bermudan options e.g. Bermudan Power Reverse Dual Currency options, where we take into account stochastic domestic and foreign interest rates on top of stochastic FX rate, hence we consider a 3-factor model. For these two numerical methods, we give an estimation of the $L^2$ -error induced by such approximations and we illustrate them with market-based examples that highlight the speed of such methods.

preprint2020arXiv

Stationary Heston model: Calibration and Pricing of exotics using Product Recursive Quantization

A major drawback of the Standard Heston model is that its implied volatility surface does not produce a steep enough smile when looking at short maturities. For that reason, we introduce the Stationary Heston model where we replace the deterministic initial condition of the volatility by its invariant measure and show, based on calibrated parameters, that this model produce a steeper smile for short maturities than the Standard Heston model. We also present numerical solution based on Product Recursive Quantization for the evaluation of exotic options (Bermudan and Barrier options).

preprint2019arXiv

Characterization of probability distribution convergence in Wasserstein distance by $L^{p}$ -quantization error function

We establish conditions to characterize probability measures by their $L^{p}$ -quantization error functions in both $\mathbb{R}^{d}$ and Hilbert settings. This characterization is two-fold: static (identity of two distributions) and dynamic (convergence for the $L^p$ -Wasserstein distance). We first propose a criterion on the quantization level $N$ , valid for any norm on $\mathbb{R}^{d}$ and any order $p$ based on a geometrical approach involving the Voronoï diagram. Then, we prove that in the $L^2$ -case on a (separable) Hilbert space, the condition on the level $N$ can be reduced to $N=2$ , which is optimal. More quantization based characterization cases on dimension 1 and a discussion of the completeness of a distance defined by the quantization error function can be found in the end of this paper.

preprint2018arXiv

Weak error for nested Multilevel Monte Carlo

This article discusses MLMC estimators with and without weights, applied to nested expectations of the form E [f (E [F (Y, Z)|Y ])]. More precisely, we are interested on the assumptions needed to comply with the MLMC framework, depending on whether the payoff function f is smooth or not. A new result to our knowledge is given when f is not smooth in the development of the weak error at an order higher than 1, which is needed for a successful use of MLMC estimators with weights.

preprint2016arXiv

Multilevel Richardson-Romberg extrapolation

We propose and analyze a Multilevel Richardson-Romberg (MLRR) estimator which combines the higher order bias cancellation of the Multistep Richardson-Romberg method introduced in [Pa07] and the variance control resulting from the stratification introduced in the Multilevel Monte Carlo (MLMC) method (see [Hei01, Gi08]). Thus, in standard frameworks like discretization schemes of diffusion processes, the root mean squared error (RMSE) $\varepsilon > 0$ can be achieved with our MLRR estimator with a global complexity of $\varepsilon^{-2} \log(1/\varepsilon)$ instead of $\varepsilon^{-2} (\log(1/\varepsilon))^2$ with the standard MLMC method, at least when the weak error $\mathbf{E}[Y_h]-\mathbf{E}[Y_0]$ of the biased implemented estimator $Y_h$ can be expanded at any order in $h$ and $\|Y_h - Y_0\|_2 = O(h^{\frac{1}{2}})$ . The MLRR estimator is then halfway between a regular MLMC and a virtual unbiased Monte Carlo. When the strong error $\|Y_h - Y_0\|_2 = O(h^{\fracβ{2}})$ , $β< 1$ , the gain of MLRR over MLMC becomes even more striking. We carry out numerical simulations to compare these estimators in two settings: vanilla and path-dependent option pricing by Monte Carlo simulation and the less classical Nested Monte Carlo simulation.

preprint2016arXiv

Vibrato and automatic differentiation for high order derivatives and sensitivities of financial options

This paper deals with the computation of second or higher order greeks of financial securities. It combines two methods, Vibrato and automatic differentiation and compares with other methods. We show that this combined technique is faster than standard finite difference, more stable than automatic differentiation of second order derivatives and more general than Malliavin Calculus. We present a generic framework to compute any greeks and present several applications on different types of financial contracts: European and American options, multidimensional Basket Call and stochastic volatility models such as Heston's model. We give also an algorithm to compute derivatives for the Longstaff-Schwartz Monte Carlo method for American options. We also extend automatic differentiation for second order derivatives of options with non-twice differentiable payoff. 1. Introduction. Due to BASEL III regulations, banks are requested to evaluate the sensitivities of their portfolios every day (risk assessment). Some of these portfolios are huge and sensitivities are time consuming to compute accurately. Faced with the problem of building a software for this task and distrusting automatic differentiation for non-differentiable functions, we turned to an idea developed by Mike Giles called Vibrato. Vibrato at core is a differentiation of a combination of likelihood ratio method and pathwise evaluation. In Giles [12], [13], it is shown that the computing time, stability and precision are enhanced compared with numerical differentiation of the full Monte Carlo path. In many cases, double sensitivities, i.e. second derivatives with respect to parameters, are needed (e.g. gamma hedging). Finite difference approximation of sensitivities is a very simple method but its precision is hard to control because it relies on the appropriate choice of the increment. Automatic differentiation of computer programs bypass the difficulty and its computing cost is similar to finite difference, if not cheaper. But in finance the payoff is never twice differentiable and so generalized derivatives have to be used requiring approximations of Dirac functions of which the precision is also doubtful. The purpose of this paper is to investigate the feasibility of Vibrato for second and higher derivatives. We will first compare Vibrato applied twice with the analytic differentiation of Vibrato and show that it is equivalent, as the second is easier we propose the best compromise for second derivatives: Automatic Differentiation of Vibrato. In [8], Capriotti has recently investigated the coupling of different mathematical methods -- namely pathwise and likelihood ratio methods -- with an Automatic differ

preprint2016arXiv

Weighted Multilevel Langevin Simulation of Invariant Measures

We investigate a weighted Multilevel Richardson-Romberg extrapolation for the ergodic approximation of invariant distributions of diffusions adapted from the one introduced in~[Lemaire-Pagès, 2013] for regular Monte Carlo simulation. In a first result, we prove under weak confluence assumptions on the diffusion, that for any integer $R\ge2$ , the procedure allows us to attain a rate $n^{\frac{R}{2R+1}}$ whereas the original algorithm convergence is at a weak rate $n^{1/3}$ . Furthermore, this is achieved without any explosion of the asymptotic variance. In a second part, under stronger confluence assumptions and with the help of some second order expansions of the asymptotic error, we go deeper in the study by optimizing the choice of the parameters involved by the method. In particular, for a given $\varepsilon\textgreater{}0$, we exhibit some semi-explicit parameters for which the number of iterations of the Euler scheme required to attain a Mean-Squared Error lower than $\varepsilon^2$ is about $\varepsilon^{-2}\log(\varepsilon^{-1})$ . Finally, we numerically this Multilevel Langevin estimator on several examples including the simple one-dimensional Ornstein-Uhlenbeck process but also on a high dimensional diffusion motivated by a statistical problem. These examples confirm the theoretical efficiency of the method.

preprint2015arXiv

Greedy vector quantization

We investigate the greedy version of the $L^p$ -optimal vector quantization problem for an $\mathbb{R}^d$ -valued random vector $X\!\in L^p$ . We show the existence of a sequence $(a_N)_{N\ge 1}$ such that $a_N$ minimizes $a\mapsto\big \|\min_{1\le i\le N-1}|X-a_i|\wedge |X-a|\big\|_{L^p}$ ( $L^p$ -mean quantization error at level $N$ induced by $(a_1,\ldots,a_{N-1},a)$ ). We show that this sequence produces $L^p$ -rate optimal $N$ -tuples $a^{(N)}=(a_1,\ldots,a_{_N})$ ( $i.e.$ the $L^p$ -mean quantization error at level $N$ induced by $a^{(N)}$ goes to $0$ at rate $N^{-\frac 1d}$ ). Greedy optimal sequences also satisfy, under natural additional assumptions, the distortion mismatch property: the $N$ -tuples $a^{(N)}$ remain rate optimal with respect to the $L^q$ -norms, $p\le q <p+d$ . Finally, we propose optimization methods to compute greedy sequences, adapted from usual Lloyd's I and Competitive Learning Vector Quantization procedures, either in their deterministic (implementable when $d=1$ ) or stochastic versions.

preprint2015arXiv

Recursive marginal quantization of the Euler scheme of a diffusion process

We propose a new approach to quantize the marginals of the discrete Euler diffusion process. The method is built recursively and involves the conditional distribution of the marginals of the discrete Euler process. Analytically, the method raises several questions like the analysis of the induced quadratic quantization error between the marginals of the Euler process and the proposed quantizations. We show in particular that at every discretization step $t\_k$ of the Euler scheme, this error is bounded by the cumulative quantization errors induced by the Euler operator, from times $t\_0=0$ to time $t\_k$ . For numerics, we restrict our analysis to the one dimensional setting and show how to compute the optimal grids using a Newton-Raphson algorithm. We then propose a closed formula for the companion weights and the transition probabilities associated to the proposed quantizations. This allows us to quantize in particular diffusion processes in local volatility models by reducing dramatically the computational complexity of the search of optimal quantizers while increasing their computational precision with respect to the algorithms commonly proposed in this framework. Numerical tests are carried out for the Brownian motion and for the pricing of European options in a local volatility model. A comparison with the Monte Carlo simulations shows that the proposed method may sometimes be more efficient (w.r.t. both computational precision and time complexity) than the Monte Carlo method.

preprint2015arXiv

Sharp rate for the dual quantization problem

In this paper we establish the sharp rate of the optimal dual quantization problem. The notion of dual quantization was recently introduced in the paper [8], where it was shown that, at least in an Euclidean setting, dual quantizers are based on a Delaunay triangulation, the dual counterpart of the Voronoi tessellation on which "regular" quantization relies. Moreover, this new approach shares an intrinsic stationarity property, which makes it very valuable for numerical applications. We establish in this paper the counterpart for dual quantization of the celebrated Zador theorem, which describes the sharp asymptotics for the quantization error when the quantizer size tends to infinity. The proof of this theorem relies among others on an extension of the so-called Pierce Lemma by means of a random quantization argument.

preprint2014arXiv

Convex order for path-dependent derivatives: a dynamic programming approach

We investigate the (functional) convex order of for various continuous martingale processes, either with respect to their diffusions coefficients for Lévy-driven SDEs or their integrands for stochastic integrals. Main results are bordered by counterexamples. Various upper and lower bounds can be derived for path wise European option prices in local volatility models. In view of numerical applications, we adopt a systematic (and symmetric) methodology: (a) propagate the convexity in a {\em simulatable} dominating/dominated discrete time model through a backward induction (or linear dynamical principle); (b) Apply functional weak convergence results to numerical schemes/time discretizations of the continuous time martingale satisfying (a) in order to transfer the convex order properties. Various bounds are derived for European options written on convex pathwise dependent payoffs. We retrieve and extend former results obtains by several authors since the seminal 1985 paper by Hajek . In a second part, we extend this approach to Optimal Stopping problems using a that the Snell envelope satisfies (a') a Backward Dynamical Programming Principle to propagate convexity in discrete time; (b') satisfies abstract convergence results under non-degeneracy assumption on filtrations. Applications to the comparison of American option prices on convex pathwise payoff processes are given obtained by a purely probabilistic arguments.

preprint2014arXiv

Functional quantization-based stratified sampling methods

In this article, we propose several quantization-based stratified sampling methods to reduce the variance of a Monte Carlo simulation. Theoretical aspects of stratification lead to a strong link between optimal quadratic quantization and the variance reduction that can be achieved with stratified sampling. We first put the emphasis on the consistency of quantization for partitioning the state space in stratified sampling methods in both finite and infinite dimensional cases. We show that the proposed quantization-based strata design has uniform efficiency among the class of Lipschitz continuous functionals. Then a stratified sampling algorithm based on product functional quantization is proposed for path-dependent functionals of multi-factor diffusions. The method is also available for other Gaussian processes such as Brownian bridge or Ornstein-Uhlenbeck processes. We derive in detail the case of Ornstein-Uhlenbeck processes. We also study the balance between the algorithmic complexity of the simulation and the variance reduction factor

preprint2013arXiv

Pointwise convergence of the Lloyd algorithm in higher dimension

We establish the pointwise convergence of the iterative Lloyd algorithm, also known as $k$ -means algorithm, when the quadratic quantization error of the starting grid (with size $N\ge 2$ ) is lower than the minimal quantization error with respect to the input distribution is lower at level $N-1$ . Such a protocol is known as the splitting method and allows for convergence even when the input distribution has an unbounded support. We also show under very light assumption that the resulting limiting grid still has full size $N$ . These results are obtained without continuity assumption on the input distribution. A variant of the procedure taking advantage of the asymptotic of the optimal quantizer radius is proposed which always guarantees the boundedness of the iterated grids.

preprint2012arXiv

Asymptotics of the maximal radius of an $L^r$ -optimal sequence of quantizers

Let $P$ be a probability distribution on $\mathbb{R}^d$ (equipped with an Euclidean norm $|\cdot|$ ). Let $r> 0$ and let $(α_n)_{n \geq1}$ be an (asymptotically) $L^r(P)$ -optimal sequence of $n$ -quantizers. We investigate the asymptotic behavior of the maximal radius sequence induced by the sequence $(α_n)_{n \geq1}$ defined for every $n \geq1$ by $ρ(α_n) = \max{|a|, a \inα_n}$ . When $\card(\supp(P))$ is infinite, the maximal radius sequence goes to $\sup{|x|, x \in\operatorname{supp}(P)}$ as $n$ goes to infinity. We then give the exact rate of convergence for two classes of distributions with unbounded support: distributions with hyper-exponential tails and distributions with polynomial tails. In the one-dimensional setting, a sharp rate and constant are provided for distributions with hyper-exponential tails.

preprint2012arXiv

Functional co-monotony of processes with applications to peacocks and barrier options

We show that several general classes of stochastic processes satisfy a functional co-monotony principle, including processes with independent increments, Brownian diffusions, Liouville processes. As a first application, we recover some recent results about peacock processes obtained by Hirsch et al. which were themselves motivated by a former work of Carr et al. about the sensitivity of Asian Call options with respect to their volatility and residual maturity (seniority). We also derive semi-universal bounds for various barrier options.

preprint2012arXiv

Optimal posting price of limit orders: learning by trading

Considering that a trader or a trading algorithm interacting with markets during continuous auctions can be modeled by an iterating procedure adjusting the price at which he posts orders at a given rhythm, this paper proposes a procedure minimizing his costs. We prove the a.s. convergence of the algorithm under assumptions on the cost function and give some practical criteria on model parameters to ensure that the conditions to use the algorithm are fulfilled (using notably the co-monotony principle). We illustrate our results with numerical experiments on both simulated data and using a financial market dataset.

preprint2012arXiv

Stochastic Approximation with Averaging Innovation Applied to Finance

The aim of the paper is to establish a convergence theorem for multi-dimensional stochastic approximation when the "innovations" satisfy some "light" averaging properties in the presence of a pathwise Lyapunov function. These averaging assumptions allow us to unify apparently remote frameworks where the innovations are simulated (possibly deterministic like in Quasi-Monte Carlo simulation) or exogenous (like market data) with ergodic properties. We propose several fields of applications and illustrate our results on five examples mainly motivated by Finance.

preprint2011arXiv

GPGPUs in computational finance: Massive parallel computing for American style options

The pricing of American style and multiple exercise options is a very challenging problem in mathematical finance. One usually employs a Least-Square Monte Carlo approach (Longstaff-Schwartz method) for the evaluation of conditional expectations which arise in the Backward Dynamic Programming principle for such optimal stopping or stochastic control problems in a Markovian framework. Unfortunately, these Least-Square Monte Carlo approaches are rather slow and allow, due to the dependency structure in the Backward Dynamic Programming principle, no parallel implementation; whether on the Monte Carlo levelnor on the time layer level of this problem. We therefore present in this paper a quantization method for the computation of the conditional expectations, that allows a straightforward parallelization on the Monte Carlo level. Moreover, we are able to develop for AR(1)-processes a further parallelization in the time domain, which makes use of faster memory structures and therefore maximizes parallel execution. Finally, we present numerical results for a CUDA implementation of this methods. It will turn out that such an implementation leads to an impressive speed-up compared to a serial CPU implementation.

preprint2010arXiv

Optimal split of orders across liquidity pools: a stochastic algorithm approach

Evolutions of the trading landscape lead to the capability to exchange the same financial instrument on different venues. Because of liquidity issues, the trading firms split large orders across several trading destinations to optimize their execution. To solve this problem we devised two stochastic recursive learning procedures which adjust the proportions of the order to be sent to the different venues, one based on an optimization principle, the other on some reinforcement ideas. Both procedures are investigated from a theoretical point of view: we prove a.s. convergence of the optimization algorithm under some light ergodic (or "averaging") assumption on the input data process. No Markov property is needed. When the inputs are i.i.d. we show that the convergence rate is ruled by a Central Limit Theorem. Finally, the mutual performances of both algorithms are compared on simulated and real data with respect to an "oracle" strategy devised by an "insider" who knows a priori the executed quantities by every venues.

preprint2009arXiv

Approximation of the distribution of a stationary Markov process with application to option pricing

We build a sequence of empirical measures on the space D(R_+,R^d) of R^d-valued càdlàg functions on R_+ in order to approximate the law of a stationary R^d-valued Markov and Feller process (X_t). We obtain some general results of convergence of this sequence. Then, we apply them to Brownian diffusions and solutions to Lévy driven SDE's under some Lyapunov-type stability assumptions. As a numerical application of this work, we show that this procedure gives an efficient way of option pricing in stochastic volatility models.

preprint2006arXiv

Distortion mismatch in the quantization of probability measures

We elucidate the asymptotics of the L^s-quantization error induced by a sequence of L^r-optimal n-quantizers of a probability distribution P on R^d when s>r. In particular we show that under natural assumptions, the optimal rate is preserved as long as s<r+d (and for every s in the case of a compactly supported distribution). We derive some applications of these results to the error bounds for quantization based quadrature formulae in numerical integration on R^d and on the Wiener space.

preprint2006arXiv

Moment estimates for Lévy Processes

For real Lévy processes $(X\_t)\_{t \geq 0}$ having no Brownian component with Blumenthal-Getoor index $β$ , the estimate $\E \sup\_{s \leq t} | X\_s - a\_p s |^p \leq C\_p t$ for every $t \in [0,1]$ and suitable $a\_p \in \R$ has been established by Millar \cite{MILL} for $β< p \leq 2$ provided $X\_1 \in L^p$ . We derive extensions of these estimates to the cases $p > 2$ and $p \leqβ$ .

preprint2005arXiv

A penalized bandit algorithm

We study a two armed-bandit algorithm with penalty. We show the convergence of the algorithm and establish the rate of convergence. For some choices of the parameters, we obtain a central limit theorem in which the limit distribution is characterized as the unique stationary distribution of a discontinuous Markov process.

preprint2005arXiv

A two armed bandit type problem revisited

In a recent paper, M. Benaïm and G. Ben Arous solve a multi-armed bandit problem arising in the theory of learning in games. We propose an short elementary proof of this result based on a variant of the Kronecker Lemma.

preprint2005arXiv

Functional quantization and metric entropy for Riemann-Liouville processes

We derive a high-resolution formula for the $L^2$ -quantization errors of Riemann-Liouville processes and the sharp Kolmogorov entropy asymptotics for related Sobolev balls. We describe a quantization procedure which leads to asymptotically optimal functional quantizers. Regular variation of the eigenvalues of the covariance operator plays a crucial role.

preprint2005arXiv

How fast is the bandit?

In this paper we investigate the rate of convergence of the so-called two-armed bandit algorithm in a financial context of asset allocation. The behaviour of the algorithm turns out to be highly non-standard: no CLT whatever the time scale, possible existence of two rate regimes.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Gilles Pagès

Research overview

Is this your work?

34 publications

Strong Solutions and Quantization-Based Numerical Schemes for a Class of Non-Markovian Volatility Models

Langevin algorithms for Markovian Neural Networks and Deep Stochastic control

Convergence of Langevin-Simulated Annealing algorithms with multiplicative noise

Convergence of Langevin-Simulated Annealing algorithms with multiplicative noise II: Total Variation

Functional convex order for the scaled McKean-Vlasov processes

Convergence rate of optimal quantization grids and application to empirical measure

New approach to greedy vector quantization

New Weak Error bounds and expansions for Optimal Quantization

Quantization-based Bermudan option pricing in the $FX$ world

Stationary Heston model: Calibration and Pricing of exotics using Product Recursive Quantization

Characterization of probability distribution convergence in Wasserstein distance by $L^{p}$ -quantization error function

Weak error for nested Multilevel Monte Carlo

Multilevel Richardson-Romberg extrapolation

Vibrato and automatic differentiation for high order derivatives and sensitivities of financial options

Weighted Multilevel Langevin Simulation of Invariant Measures

Greedy vector quantization

Recursive marginal quantization of the Euler scheme of a diffusion process

Sharp rate for the dual quantization problem

Convex order for path-dependent derivatives: a dynamic programming approach

Functional quantization-based stratified sampling methods

Pointwise convergence of the Lloyd algorithm in higher dimension

Asymptotics of the maximal radius of an $L^r$ -optimal sequence of quantizers

Functional co-monotony of processes with applications to peacocks and barrier options

Optimal posting price of limit orders: learning by trading

Stochastic Approximation with Averaging Innovation Applied to Finance

GPGPUs in computational finance: Massive parallel computing for American style options

Optimal split of orders across liquidity pools: a stochastic algorithm approach

Approximation of the distribution of a stationary Markov process with application to option pricing

Distortion mismatch in the quantization of probability measures

Moment estimates for Lévy Processes

A penalized bandit algorithm

A two armed bandit type problem revisited

Functional quantization and metric entropy for Riemann-Liouville processes

How fast is the bandit?

See the researcher in context

Building this map preview

Gilles Pagès

Research overview

Is this your work?

34 publications

Strong Solutions and Quantization-Based Numerical Schemes for a Class of Non-Markovian Volatility Models

Langevin algorithms for Markovian Neural Networks and Deep Stochastic control

Convergence of Langevin-Simulated Annealing algorithms with multiplicative noise

Convergence of Langevin-Simulated Annealing algorithms with multiplicative noise II: Total Variation

Functional convex order for the scaled McKean-Vlasov processes

Convergence rate of optimal quantization grids and application to empirical measure

New approach to greedy vector quantization

New Weak Error bounds and expansions for Optimal Quantization

Quantization-based Bermudan option pricing in the FXFXFX world

Stationary Heston model: Calibration and Pricing of exotics using Product Recursive Quantization

Characterization of probability distribution convergence in Wasserstein distance by LpL^{p}Lp-quantization error function

Weak error for nested Multilevel Monte Carlo

Multilevel Richardson-Romberg extrapolation

Vibrato and automatic differentiation for high order derivatives and sensitivities of financial options

Weighted Multilevel Langevin Simulation of Invariant Measures

Greedy vector quantization

Recursive marginal quantization of the Euler scheme of a diffusion process

Sharp rate for the dual quantization problem

Convex order for path-dependent derivatives: a dynamic programming approach

Functional quantization-based stratified sampling methods

Pointwise convergence of the Lloyd algorithm in higher dimension

Asymptotics of the maximal radius of an LrL^rLr-optimal sequence of quantizers

Functional co-monotony of processes with applications to peacocks and barrier options

Optimal posting price of limit orders: learning by trading

Stochastic Approximation with Averaging Innovation Applied to Finance

GPGPUs in computational finance: Massive parallel computing for American style options

Optimal split of orders across liquidity pools: a stochastic algorithm approach

Approximation of the distribution of a stationary Markov process with application to option pricing

Distortion mismatch in the quantization of probability measures

Moment estimates for Lévy Processes

A penalized bandit algorithm

A two armed bandit type problem revisited

Functional quantization and metric entropy for Riemann-Liouville processes

How fast is the bandit?

See the researcher in context

Building this map preview

Quantization-based Bermudan option pricing in the $FX$ world

Characterization of probability distribution convergence in Wasserstein distance by $L^{p}$ -quantization error function

Asymptotics of the maximal radius of an $L^r$ -optimal sequence of quantizers