MathGroup Archive 2006

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

Search the Archive

Re: Singularity handling in NIntegrate

  • To: mathgroup at smc.vnet.net
  • Subject: [mg70859] Re: Singularity handling in NIntegrate
  • From: "Andrew Moylan" <andrew.j.moylan at gmail.com>
  • Date: Mon, 30 Oct 2006 05:32:19 -0500 (EST)
  • References: <ehv7ph$ft1$1@smc.vnet.net><ei18qn$8nk$1@smc.vnet.net>

Thanks Dimitris,

I've had a look at those sources. Actually, that previous thread that
you linked to was also started by me ;-).

Anton's response does suggest that the answer to my first question is
just "yes". It's really the second question (what are the coordinate
transformations used when a singularity is detected?) that I'm most
interested in. (Anton's response describes the compactifying
transformation used for infinite regions of integration.)

Can anyone shed any light on this second question?

Andrew


dimitris wrote:
> Information[SingularityDepth]
> "SingularityDepth is an option for NIntegrate. SingularityDepth -> n
> specifies that n recursive subdivisions can be done before
> a change of variable is used at the endpoints of the interval of
> integration."
> Attributes[SingularityDepth] = {Protected}
>
> See
>
> http://groups.google.com/group/comp.soft-sys.math.mathematica/browse_thread/thread/11d40e17edbea9e9/b223cb3882539ac8#b223cb3882539ac8
>
> http://library.wolfram.com/infocenter/Conferences/5365/
>
> Hope these helps.
>
> Regards
> Dimitris
>
>
>
> Andrew Moylan wrote:
> > Hi all,
> >
> > I've been searching around for information about the particular
> > method(s) that NIntegrate uses to handle singularities at the endpoints
> > of a numerical integration. In particular,
> >
> > 1. How is the SingularityDepth option used? Is it just the number of
> > recursive bisections of the first/last subinterval of integration
> > before a singularity is assumed?
> >
> > 2. What is/are the coordinate transformation(s) used when a singularity
> > is detected?
> >
> > Thanks for any help you can give me with these questions.
> > 
> > Cheers,
> > 
> > Andrew


  • Prev by Date: Re: "unload" a package
  • Next by Date: Matching with a "Flat" function
  • Previous by thread: Re: Singularity handling in NIntegrate
  • Next by thread: Request for benchmark results ...