public inbox for gsl-discuss@sourceware.org
 help / color / mirror / Atom feed
From: Toan T Nguyen <toan.nguyen@physics.gatech.edu>
To: gsl-discuss@sourceware.org
Subject: Monte-Carlo integration good ?
Date: Tue, 17 Aug 2010 09:27:00 -0000	[thread overview]
Message-ID: <4C69A6F5.4050001@physics.gatech.edu> (raw)

Hello all,

I'm seeking advice as on how reliable the result of my numerical
integration is. I've been using the GSL Monte-Carlo integration
routines to perform a 3D integration of a function with
singularity at the origin. The integrand contains something similar to

exp( - k x) / x

The integration range is for 0 < x < 1. (However, since k ~ 50, only the 
range
0< x < 0.02 contributes most to the integration.)

Both GSL VEGAS and MISER algorithm gives more or less the same result.
But I'm a tiny bit suspect of the number that GSL gives.

My question is how reliable is this result. Are VEGAS/MISER algorithms
able to automatically focus on the 0<x<0.02 range instead of the full
integration range 0 < x < 1? Do I need to manually adjust the parameters
of the integration subroutine to handle this situation ? So far, I simply
follow the Monte-Carlo integration example to do the integration.

Thanks

Toan


                 reply	other threads:[~2010-08-17  9:27 UTC|newest]

Thread overview: [no followups] expand[flat|nested]  mbox.gz  Atom feed

Reply instructions:

You may reply publicly to this message via plain-text email
using any one of the following methods:

* Save the following mbox file, import it into your mail client,
  and reply-to-all from there: mbox

  Avoid top-posting and favor interleaved quoting:
  https://en.wikipedia.org/wiki/Posting_style#Interleaved_style

* Reply using the --to, --cc, and --in-reply-to
  switches of git-send-email(1):

  git send-email \
    --in-reply-to=4C69A6F5.4050001@physics.gatech.edu \
    --to=toan.nguyen@physics.gatech.edu \
    --cc=gsl-discuss@sourceware.org \
    /path/to/YOUR_REPLY

  https://kernel.org/pub/software/scm/git/docs/git-send-email.html

* If your mail client supports setting the In-Reply-To header
  via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).