Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2006
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2006

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

Search the Archive

Re: Singularity-handling transformation employed by NIntegrate

  • To: mathgroup at smc.vnet.net
  • Subject: [mg71214] Re: Singularity-handling transformation employed by NIntegrate
  • From: "antononcube" <antononcube at gmail.com>
  • Date: Fri, 10 Nov 2006 06:38:29 -0500 (EST)
  • References: <eiut46$h63$1@smc.vnet.net>

NIntegrate uses the so called IMT Rule (the Iri, Moriguti, Takasawa
rule) that changes the integration variable with this transformation
over the interval [xmin, xmax]:

Clear[\[Phi]]
\[Phi][{xmin_, xmax_}, t_] := xmin + (xmax - xmin)*E^(1 - 1/t);
\[Phi][t_] := \[Phi][{0, 1}, t]


It is instructive to look at the plot of the transformation:

Plot[\[Phi][t], {t, 0, 1}, PlotRange -> All, AxesOrigin -> {-0.05,
-0.05}]

and at the plot of its derivative (Jacobian):

Plot[Evaluate[D[\[Phi][t], t]], {t, 0, 1}, PlotRange -> All]

The IMT rule "flattens" the singularity, but the points are clustered
closer to the singular point. That is why in many(most) cases extra
precision is required for the evaluation of the integrand.

Anton Antonov,
Wolfram Research, Inc.


Andrew Moylan wrote:
> Hi all,
>
> Further to a previous question: Can anyone tell me what transformation
> NIntegrate employs when a singularity is detected (i.e., when the
> number of recursive subdivisions specified by SingularityDepth is
> reached)?
> 
> Cheers,
> Andrew


  • Prev by Date: Re: linux font rotation in graphics & (strange) black graphics borders
  • Next by Date: Re: Simple Problem with Polyhedron::Hexahedron
  • Previous by thread: Singularity-handling transformation employed by NIntegrate
  • Next by thread: Re: Singularity-handling transformation employed by NIntegrate