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

Re: Quick integral.

• To: mathgroup at smc.vnet.net
• Subject: [mg73489] Re: Quick integral.
• From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
• Date: Tue, 20 Feb 2007 06:14:11 -0500 (EST)
• Organization: The Open University, Milton Keynes, UK
• References: <erc23r$mv4$1@smc.vnet.net>

Jouvenot, Fabrice wrote:
> Thanks for your answers.
> 
> As some of you asked for the notebook here is it :
> http://www.onyrium.net/Mathematica/CoordScat.nb
> 
> At the end there is an integral named toto[], we launch it and have an
> evaluation time of how this implementation in a large notebook will
> affect time.
> 
> We want to minimize a really lot this time. We try to make more physics
> approximations, but there is certainly mathematica things we can do to
> quicker everythings.
> 
> One more things to know : we do not need a large accuracy on the
> results, large errors are ok.
> 
> Thanks for your help.
> 
> 
> 
> Fabrice.
> 
Hi Fabrice,

adding the option Method -> Trapezoidal to the definition of the 
function toto might slightly improve the performance on average.

toto[l_] := NIntegrate[dNdÏ?[l, angleÏ?[Ψ, thc, θsource], 
thc]*AcceptancePMT[angleLightOM[vOM, vlight[thc, Ψ, θsource, l]]], {Ψ, 
0, 2*Pi},
     AccuracyGoal -> 1, PrecisionGoal -> 0, SingularityDepth -> 1, 
MaxPoints -> 3, Method -> Trapezoidal];

HTH,
Jean-Marc