MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Adaptive Monte Carlo integration algorithm(s) in Mathematica6

Dear Denis,

The Monte Carlo algorithms are described in the advanced

See "Global Adaptive Monte Carlo and Quasi Monte Carlo Strategies" in

and "MonteCarloRule" in

(You can find these documents in the Help Browser if you have version

NIntegrate's adaptive Monte Carlo algorithm is closer to MISER
(recursive stratified integration) than to VEGAS (importance

I have an implementation of an adaptive importance sampling algorithm,
please, email me if you would like to try it.

Still, one can use the MonteCarloRule's option "PointGenerator" to
implement an importance sampling algorithm (and use it either with
crude or adaptive Monte Carlo algorithm).

As for NIntegrate in 5.2, one can actually do stratified integration
through the range specification. See again the integration strategies
document above.

Anton Antonov,
Wolfram Research, Inc.

On May 8, 4:48 am, Denis <denis.... at> wrote:
> I'm happy with Mathematica 5.2 in most respects, except it doesn't
> offer any *adaptive* Monte Carlo methods for multi-dimensional
> integration.  Simple pseudo-random MonteCarlo and the non-adaptive
> quasi-Monte Carlo are no match for the 6-dimensional integrals I need
> to calculate.  Right now, as a way around that shortcoming, I am using
> a third-party C code for LePage's VEGAS algorithm and run it in
> Mathematica courtesy of MathLink.  It seems to work quite alright.
> Now that Mathematica6 seems to offer adaptive Monte Carlo natively
> within NIntegrate, my question is - which algorithm(s)?  I would like
> to know whether Mathematica implemented their own version of VEGAS or
> perhaps some other adaptive Monte Carlo algorithm...  Does anyone know
> any specifics about this?
> Thanks,
> Denis

  • Prev by Date: Re: can I keep old version running after update to Ver 6?
  • Next by Date: elliptic integral
  • Previous by thread: Adaptive Monte Carlo integration algorithm(s) in Mathematica6
  • Next by thread: Using FindFit with constraints ??