[Quantlib-dev] Implementation of Milstein discretization scheme

[<fho...@gm...>]

Is there any suggestion how to implement a Milstein scheme? The main difference between Euler and Milstein discretization is that for, e.g., a geometric Brownian motion it incorporates a drift correction which involves the realization of the standard normal variate which also applies to the corresponding diffusion part. One would have to pass a dw parameter to the discretization itself. In the multidimensional case one would have to pass off-diagonal terms as well (see Peter Jäckel's book on 'Monte Carlo Methods in Finance'). Predictor-corrector mechanisms seem to be covered by the already existing setup for the Euler discretization. Any suggestions?

Regards
Frank
-- 
Ist Ihr Browser Vista-kompatibel? Jetzt die neuesten 
Browser-Versionen downloaden: http://www.gmx.net/de/go/browser