MathGroup Archive 2008

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

Search the Archive

Probable bug: EvaluationMonitor with NIntegrate, "ExtrapolatingOscillatory",

  • To: mathgroup at smc.vnet.net
  • Subject: [mg85684] Probable bug: EvaluationMonitor with NIntegrate, "ExtrapolatingOscillatory",
  • From: Andrew Moylan <andrew.j.moylan at gmail.com>
  • Date: Tue, 19 Feb 2008 01:58:05 -0500 (EST)

For NIntegrate[BesselJ[0, x^2], {x, 0, Infinity}, EvaluationMonitor :>
bla], EvaluationMonitor is not being called with x set to the actual
values at which the integrand is being evaluated.

To see this:

dummy[x_?NumericQ] := (Sow[x, 1]; 1)
ListPlot /@ (NIntegrate[dummy[x] BesselJ[0, x^2], {x, 0, Infinity},
EvaluationMonitor :> Sow[x, 2]] // Reap // Last)

Presumably the reason for the bug is that method
"ExtrapolatingOscillatory" (which is being used automatically for this
integrand) is (correctly) integrating between the roots of BesselJ[0,
x^2], but EvaluationMonitor is being evaluated with x set to the roots
of BesselJ[0, x], i.e., with x being set to the value of an internal
transformed integration variable.


  • Prev by Date: Re: Garbage collection and mutable data structures
  • Next by Date: Re: [functional approach should give] an even faster way to normalize
  • Previous by thread: Re: Localizing Large Numbers of Variables in a DynamicModule:
  • Next by thread: Re: [functional approach should give] an even faster way to normalize