A fast algorithm for the two dimensional HJB equation of stochastic control

Frédéric Bonnans, Elisabeth Ottenwaelter and Hasnaa Zidani
2004
Publication type:
Paper in peer-reviewed journals
Journal:
ESAIM: Mathematical Modelling and Numerical Analysis, vol. 38(4), pp. 723-735
Abstract:
This paper analyses the implementation of the generalized finite differences method for the HJB equation of stochastic control, introduced by two of the authors in [Bonnans and Zidani, SIAM J. Numer. Anal. 41 (2003) 1008–1021]. The computation of coefficients needs to solve at each point of the grid (and for each control) a linear programming problem. We show here that, for two dimensional problems, this linear programming problem can be solved in O(p max) operations, where p max is the size of the stencil. The method is based on a walk on the Stern-Brocot tree, and on the related filling of the set of positive semidefinite matrices of size two.
BibTeX:
@article{Bon-Ott-Zid-2004,
    author={Frédéric Bonnans and Elisabeth Ottenwaelter and Hasnaa 
           Zidani },
    title={A fast algorithm for the two dimensional HJB equation of 
           stochastic control },
    doi={10.1051/m2an:2004034 },
    journal={ESAIM: Mathematical Modelling and Numerical Analysis },
    year={2004 },
    volume={38(4) },
    pages={723--735},
}