An Algorithm for Finding Best Matches in Logarithmic Expected Time
ACM Transactions on Mathematical Software (TOMS)
Foundations of Quantization for Probability Distributions
Foundations of Quantization for Probability Distributions
An Object-Oriented Random-Number Package with Many Long Streams and Substreams
Operations Research
Regression methods for pricing complex American-style options
IEEE Transactions on Neural Networks
Special Issue for the Workshop on High Performance Computational Finance
Concurrency and Computation: Practice & Experience
Pricing American options with least squares Monte Carlo on GPUs
WHPCF '13 Proceedings of the 6th Workshop on High Performance Computational Finance
Recent progress and challenges in exploiting graphics processors in computational fluid dynamics
The Journal of Supercomputing
Hi-index | 0.00 |
The pricing of American style and multiple exercise options is a very challenging problem in mathematical finance. One usually employs a least squares 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 squares MC approaches are rather slow and allow, because of the dependency structure in the backward dynamic programming principle, no parallel implementation neither on the MC level nor 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 MC level. Moreover, we are able to develop for first-order autoregressive processes a further parallelization in the time domain, which makes use of faster memory structures and therefore maximizes parallel execution. Furthermore, we discuss the generation of random numbers in parallel on a GPGPU architecture, which is the crucial tool for the application of massive parallel computing architectures in mathematical finance. Finally, we present numerical results for a CUDA implementation of these methods. It will turn out that such an implementation leads to an impressive speed-up compared with a serial CPU implementation. Copyright © 2011 John Wiley & Sons, Ltd.