BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Xavier Warin

Xavier Warin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.OCmath.PRq-fin.RMMachine Learningmath.NAq-fin.CPq-fin.PMNeural and Evolutionary Computingq-fin.PRComputationecon.THmath.APNumerical Analysisphysics.soc-phPopulations and Evolution

Catalog footprint

What is connected

21works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

21 published item(s)

preprint2022arXiv

Deep learning for efficient frontier calculation in finance

We propose deep neural network algorithms to calculate efficient frontier in some Mean-Variance and Mean-CVaR portfolio optimization problems. We show that we are able to deal with such problems when both the dimension of the state and the dimension of the control are high. Adding some additional constraints, we compare different formulations and show that a new projected feedforward network is able to deal with some global constraints on the weights of the portfolio while outperforming classical penalization methods. All developed formulations are compared in between. Depending on the problem and its dimension, some formulations may be preferred.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

DeepSets and their derivative networks for solving symmetric PDEs

Machine learning methods for solving nonlinear partial differential equations (PDEs) are hot topical issues, and different algorithms proposed in the literature show efficient numerical approximation in high dimension. In this paper, we introduce a class of PDEs that are invariant to permutations, and called symmetric PDEs. Such problems are widespread, ranging from cosmology to quantum mechanics, and option pricing/hedging in multi-asset market with exchangeable payoff. Our main application comes actually from the particles approximation of mean-field control problems. We design deep learning algorithms based on certain types of neural networks, named PointNet and DeepSet (and their associated derivative networks), for computing simultaneously an approximation of the solution and its gradient to symmetric PDEs. We illustrate the performance and accuracy of the PointNet/DeepSet networks compared to classical feedforward ones, and provide several numerical results of our algorithm for the examples of a mean-field systemic risk, mean-variance problem and a min/max linear quadratic McKean-Vlasov control problem.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Incentives, lockdown, and testing: from Thucydides's analysis to the COVID-19 pandemic

In this work, we provide a general mathematical formalism to study the optimal control of an epidemic, such as the COVID-19 pandemic, via incentives to lockdown and testing. In particular, we model the interplay between the government and the population as a principal-agent problem with moral hazard, à la Cvitanić, Possamaï, and Touzi [27], while an epidemic is spreading according to dynamics given by compartmental stochastic SIS or SIR models, as proposed respectively by Gray, Greenhalgh, Hu, Mao, and Pan [45] and Tornatore, Buccellato, and Vetro [88]. More precisely, to limit the spread of a virus, the population can decrease the transmission rate of the disease by reducing interactions between individuals. However, this effort, which cannot be perfectly monitored by the government, comes at social and monetary cost for the population. To mitigate this cost, and thus encourage the lockdown of the population, the government can put in place an incentive policy, in the form of a tax or subsidy. In addition, the government may also implement a testing policy in order to know more precisely the spread of the epidemic within the country, and to isolate infected individuals. In terms of technical results, we demonstrate the optimal form of the tax, indexed on the proportion of infected individuals, as well as the optimal effort of the population, namely the transmission rate chosen in response to this tax. The government's optimisation problem then boils down to solving an Hamilton-Jacobi-Bellman equation. Numerical results confirm that if a tax policy is implemented, the population is encouraged to significantly reduce its interactions. If the government also adjusts its testing policy, less effort is required on the population side, individuals can interact almost as usual, and the epidemic is largely contained by the targeted isolation of positively-tested individuals.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical resolution of McKean-Vlasov FBSDEs using neural networks

We propose several algorithms to solve McKean-Vlasov Forward Backward Stochastic Differential Equations. Our schemes rely on the approximating power of neural networks to estimate the solution or its gradient through minimization problems. As a consequence, we obtain methods able to tackle both mean field games and mean field control problems in moderate dimension. We analyze the numerical behavior of our algorithms on several examples including non linear quadratic models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Neural networks-based backward scheme for fully nonlinear PDEs

We propose a numerical method for solving high dimensional fully nonlinear partial differential equations (PDEs). Our algorithm estimates simultaneously by backward time induction the solution and its gradient by multi-layer neural networks, while the Hessian is approximated by automatic differentiation of the gradient at previous step. This methodology extends to the fully nonlinear case the approach recently proposed in \cite{HPW19} for semi-linear PDEs. Numerical tests illustrate the performance and accuracy of our method on several examples in high dimension with nonlinearity on the Hessian term including a linear quadratic control problem with control on the diffusion coefficient, Monge-Amp{è}re equation and Hamilton-Jacobi-Bellman equation in portfolio optimization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deep backward schemes for high-dimensional nonlinear PDEs

We propose new machine learning schemes for solving high dimensional nonlinear partial differential equations (PDEs). Relying on the classical backward stochastic differential equation (BSDE) representation of PDEs, our algorithms estimate simultaneously the solution and its gradient by deep neural networks. These approximations are performed at each time step from the minimization of loss functions defined recursively by backward induction. The methodology is extended to variational inequalities arising in optimal stopping problems. We analyze the convergence of the deep learning schemes and provide error estimates in terms of the universal approximation of neural networks. Numerical results show that our algorithms give very good results till dimension 50 (and certainly above), for both PDEs and variational inequalities problems. For the PDEs resolution, our results are very similar to those obtained by the recent method in \cite{weinan2017deep} when the latter converges to the right solution or does not diverge. Numerical tests indicate that the proposed methods are not stuck in poor local minimaas it can be the case with the algorithm designed in \cite{weinan2017deep}, and no divergence is experienced. The only limitation seems to be due to the inability of the considered deep neural networks to represent a solution with a too complex structure in high dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Discretization and Machine Learning Approximation of BSDEs with a Constraint on the Gains-Process

We study the approximation of backward stochastic differential equations (BSDEs for short) with a constraint on the gains process. We first discretize the constraint by applying a so-called facelift operator at times of a grid. We show that this discretely constrained BSDE converges to the continuously constrained one as the mesh grid converges to zero. We then focus on the approximation of the discretely constrained BSDE. For that we adopt a machine learning approach. We show that the facelift can be approximated by an optimization problem over a class of neural networks under constraints on the neural network and its derivative. We then derive an algorithm converging to the discretely constrained BSDE as the number of neurons goes to infinity. We end by numerical experiments. Mathematics Subject Classification (2010): 65C30, 65M75, 60H35, 93E20, 49L25.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Risk management with machine-learning-based algorithms

We propose some machine-learning-based algorithms to solve hedging problems in incomplete markets. Sources of incompleteness cover illiquidity, untradable risk factors, discrete hedging dates and transaction costs. The proposed algorithms resulting strategies are compared to classical stochastic control techniques on several payoffs using a variance criterion. One of the proposed algorithm is flexible enough to be used with several existing risk criteria. We furthermore propose a new moment-based risk criteria.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Fast and stable multivariate kernel density estimation by fast sum updating

Kernel density estimation and kernel regression are powerful but computationally expensive techniques: a direct evaluation of kernel density estimates at $M$ evaluation points given $N$ input sample points requires a quadratic $\mathcal{O}(MN)$ operations, which is prohibitive for large scale problems. For this reason, approximate methods such as binning with Fast Fourier Transform or the Fast Gauss Transform have been proposed to speed up kernel density estimation. Among these fast methods, the Fast Sum Updating approach is an attractive alternative, as it is an exact method and its speed is independent of the input sample and the bandwidth. Unfortunately, this method, based on data sorting, has for the most part been limited to the univariate case. In this paper, we revisit the fast sum updating approach and extend it in several ways. Our main contribution is to extend it to the general multivariate case for general input data and rectilinear evaluation grid. Other contributions include its extension to a wider class of kernels, including the triangular, cosine and Silverman kernels, its combination with parsimonious additive multivariate kernels, and its combination with a fast approximate k-nearest-neighbors bandwidth for multivariate datasets. Our numerical tests of multivariate regression and density estimation confirm the speed, accuracy and stability of the method. We hope this paper will renew interest for the fast sum updating approach and help solve large-scale practical density estimation and regression problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Branching diffusion representation of semilinear PDEs and Monte Carlo approximation

We provide a representation result of parabolic semi-linear PD-Es, with polynomial nonlinearity, by branching diffusion processes. We extend the classical representation for KPP equations, introduced by Skorokhod (1964), Watanabe (1965) and McKean (1975), by allowing for polynomial nonlinearity in the pair $(u, Du)$, where $u$ is the solution of the PDE with space gradient $Du$. Similar to the previous literature, our result requires a non-explosion condition which restrict to "small maturity" or "small nonlinearity" of the PDE. Our main ingredient is the automatic differentiation technique as in Henry Labordere, Tan and Touzi (2015), based on the Malliavin integration by parts, which allows to account for the nonlinearities in the gradient. As a consequence, the particles of our branching diffusion are marked by the nature of the nonlinearity. This new representation has very important numerical implications as it is suitable for Monte Carlo simulation. Indeed, this provides the first numerical method for high dimensional nonlinear PDEs with error estimate induced by the dimension-free Central limit theorem. The complexity is also easily seen to be of the order of the squared dimension. The final section of this paper illustrates the efficiency of the algorithm by some high dimensional numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical approximation of a cash-constrained firm value with investment opportunities

We consider a singular control problem with regime switching that arises in problems of optimal investment decisions of cash-constrained firms. The value function is proved to be the unique viscosity solution of the associated Hamilton-Jacobi-Bellman equation. Moreover, we give regularity properties of the value function as well as a description of the shape of the control regions. Based on these theoretical results, a numerical deterministic approximation of the related HJB variational inequality is provided. We finally show that this numerical approximation converges to the value function. This allows us to describe the investment and dividend optimal policies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Asset Liability Management problem of a nuclear operator : a numerical stochastic optimization approach

We numerically study an Asset Liability Management problem linked to the decommissioning of French nuclear power plants. We link the risk aversion of practitioners to an optimization problem. Using different price models we show that the optimal solution is linked to a de-risking management strategy similar to a concave strategy and we propose an effective heuristic to simulate the underlying optimal strategy. Besides we show that the strategy is stable with respect to the main parameters involved in the liability problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Liquidity Management with Decreasing-returns-to-scale and Secured Credit Line

This paper examines the dividend and investment policies of a cash constrained firm that has access to costly external funding. We depart from the literature by allowing the firm to issue collateralized debt to increase its investment in productive assets resulting in a performance sensitive interest rate on debt. We formulate this problem as a bi-dimensional singular control problem and use both a viscosity solution approch and a verification technique to get qualitative properties of the value function. We further solve quasi-explicitly the control problem in two special cases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Adaptive sparse grids for time dependent Hamilton-Jacobi-Bellman equations in stochastic control

We introduce some sparse grids interpolations used in Semi-Lagrangian schemes for linear and fully non-linear diffusion Hamilton Jacobi Bellman equations arising in stochastic control. We prove that the method introduced converges toward the viscosity solution of the problem and we show that some potentially high order schemes can be efficiently implemented. Numerical test in dimension $2$ to $5$ are achieved and show that deterministic methods can be used efficiently in stochastic control in moderate dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical methods for the quadratic hedging problem in Markov models with jumps

We develop algorithms for the numerical computation of the quadratic hedging strategy in incomplete markets modeled by pure jump Markov process. Using the Hamilton-Jacobi-Bellman approach, the value function of the quadratic hedging problem can be related to a triangular system of parabolic partial integro-differential equations (PIDE), which can be shown to possess unique smooth solutions in our setting. The first equation is non-linear, but does not depend on the pay-off of the option to hedge (the pure investment problem), while the other two equations are linear. We propose convergent finite difference schemes for the numerical solution of these PIDEs and illustrate our results with an application to electricity markets, where time-inhomogeneous pure jump Markov processes appear in a natural manner.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hedging Swing contract on gas markets

Swing options on the gas market are american style option where daily quantities exercices are constrained and global quantities exerciced each year constrained too. The option holder has to decide each day how much he consumes of the quantities satisfying the constraints and tries to use a strategy in order to maximize its expected profit. The pay off fonction is a spread between the spot gas market and the value of an index composed of the past average of some commodities spot or future prices. We study the valorization and the effectiveness of the dynamic hedging of such a contract.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Swing Options Valuation: a BSDE with Constrained Jumps Approach

We introduce a new probabilistic method for solving a class of impulse control problems based on their representations as Backward Stochastic Differential Equations (BSDEs for short) with constrained jumps. As an example, our method is used for pricing Swing options. We deal with the jump constraint by a penalization procedure and apply a discrete-time backward scheme to the resulting penalized BSDE with jumps. We study the convergence of this numerical method, with respect to the main approximation parameters: the jump intensity $λ$, the penalization parameter $p > 0$ and the time step. In particular, we obtain a convergence rate of the error due to penalization of order $(λp)^{α- \frac{1}{2}}, \forall α\in \left(0, \frac{1}{2}\right)$. Combining this approach with Monte Carlo techniques, we then work out the valuation problem of (normalized) Swing options in the Black and Scholes framework. We present numerical tests and compare our results with a classical iteration method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A finite dimensional approximation for pricing moving average options

We propose a method for pricing American options whose pay-off depends on the moving average of the underlying asset price. The method uses a finite dimensional approximation of the infinite-dimensional dynamics of the moving average process based on a truncated Laguerre series expansion. The resulting problem is a finite-dimensional optimal stopping problem, which we propose to solve with a least squares Monte Carlo approach. We analyze the theoretical convergence rate of our method and present numerical results in the Black-Scholes framework.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Probabilistic Numerical Method for Fully Nonlinear Parabolic PDEs

We consider the probabilistic numerical scheme for fully nonlinear PDEs suggested in \cite{cstv}, and show that it can be introduced naturally as a combination of Monte Carlo and finite differences scheme without appealing to the theory of backward stochastic differential equations. Our first main result provides the convergence of the discrete-time approximation and derives a bound on the discretization error in terms of the time step. An explicit implementable scheme requires to approximate the conditional expectation operators involved in the discretization. This induces a further Monte Carlo error. Our second main result is to prove the convergence of the latter approximation scheme, and to derive an upper bound on the approximation error. Numerical experiments are performed for the approximation of the solution of the mean curvature flow equation in dimensions two and three, and for two and five-dimensional (plus time) fully-nonlinear Hamilton-Jacobi-Bellman equations arising in the theory of portfolio optimization in financial mathematics.

0 citations0 reviewsNo code linkedOpen access