Re: Adaptive Monte Carlo integration algorithm(s) in Mathematica6
- To: mathgroup at smc.vnet.net
- Subject: [mg75730] Re: Adaptive Monte Carlo integration algorithm(s) in Mathematica6
- From: antononcube <antononcube at gmail.com>
- Date: Wed, 9 May 2007 04:36:23 -0400 (EDT)
- References: <email@example.com>
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
Wolfram Research, Inc.
On May 8, 4:48 am, Denis <denis.... at gmail.com> 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?
Prev by Date:
Re: can I keep old version running after update to Ver 6?
Next by Date:
Previous by thread:
Adaptive Monte Carlo integration algorithm(s) in Mathematica6
Next by thread:
Using FindFit with constraints ??