Re: Random points in NIntegrate
- To: mathgroup at smc.vnet.net
- Subject: [mg56501] Re: Random points in NIntegrate
- From: Ismail Turan <isturan at gmail.com>
- Date: Tue, 26 Apr 2005 21:53:43 -0400 (EDT)
- References: <d4hvel$1g7$1@smc.vnet.net> <426D56CE.9020908@arcor.de>
- Reply-to: Ismail Turan <isturan at gmail.com>
- Sender: owner-wri-mathgroup at wolfram.com
Thanks for your quick replies. I should have mentioned that I had
Mathematica 4.0. So, Sow and Reap are not defined within this version.
Furthermore, I can't use the option "'EvaluationMonitor" inside
NIntegrate.
I realized that Peter wrote a package ("GoodReapSow.nb"") which
implements Reap and Sow commands in Mathematica 4. I can use that
package even though it is not efficient according to Peter. However, I
still need to figure out how to define an option like
"EvaluationMonitor" for NIntegrate. Does anybody have an idea?
Thanks a lot,
ismail
PS: Peter I tried to run your "NUnitSphere" notebook but it seems
that it is incompatible with mathematica 4.
On 4/25/05, Peter Pein <petsie at arcor.de> wrote:
> Ismail T. wrote:
> > Hi,
> >
> > I need to do the following. I am using mathematica to take a
> > multi-dimensional integral with NIntegrate and QuasiMonteCarlo method
> > is default for such cases. There are two things that I couldn't figure
> > out. First, I want to print the random points that the mathematica
> > chooses to take each data points. This is particularly important to
> > see whether mathematica is stuck only some parts of the region of
> > integration where it assumes that the biggest cntribution comes from.
> >
> > Secondly, how can I make mathematica to take sample (random) points
> > only in a specific part of the region of integration without applying
> > any particular cut to the limits of integration?
> >
> > Thanks a lot in advance,
> >
> > ismail
> >
> I created a small notebook to illustrate tue use of Sow/Reap to get a
> list of the evaluation points. Please have a look at
> http://people.freenet.de/Peter_Berlin/Mathe/NUnitSphere/index.html
> or download an unevaluated notebook file from
> http://people.freenet.de/Peter_Berlin/Mathe/NUnitSphere/NUnitSphere_unevaluated.nb
>
> --
> Peter Pein
> Berlin
>