Speeding up this NIntegrate

Question

I'm trying out the code below and it takes quite a long time to compute and throws eincr and slwcon warnings. Is there a way to speed this up?
phase[u_, m_, [Theta]_, [Phi]_] := 
  u^3 + 2 m [Theta] (3/(4 m))^(2/3) Cos[[Phi]] u^2 + 
   m [Theta]^2 (3/(4 m))^(1/3) u - m [Theta] Cos[[Phi]];

integrand[u_, 
   m_, [Theta]_, [Phi]_] := (-2 (3/(4 m))^(1/3)
       u + [Theta] ) Exp[-I phase[u, m, [Theta], [Phi]]];

integral[[Theta]_, [Phi]_, m_, {ulo_, uhi_}] := 
 NIntegrate[integrand[u, m, [Theta], [Phi]], {u, ulo, uhi}, 
   MaxRecursion -> 20, WorkingPrecision -> 20] // 
  Timing

integraltab = 
  Monitor[Table[{[Theta], [Phi], 
     integral[[Theta], [Phi], 100, {-Infinity, Infinity}]}, {[Theta], 0, 
     Pi, Pi/32}, {[Phi], 0, 2 Pi, 
     2 Pi/63}], {[Theta], [Phi]}];

Ulrich Neumann · Answer

The integrand only depends on  Cos[[Phi], that's why you can decrease the integration range to  {[Phi], 0,  Pi} (- 50% evaluation time!)!
Include the option Method -> {Automatic, "SymbolicProcessing" -> 0} inside NIntegrate:
integral[[Theta]_, [Phi]_, m_, {ulo_, uhi_}] := 
NIntegrate[integrand[u, m, [Theta], [Phi]], {u, ulo, uhi} ,MaxRecursion -> 20, WorkingPrecision -> 20 
,Method -> {Automatic, "SymbolicProcessing" -> 0}]

the grid of 101x101 tablevalues is evaluated in nearly 75s!
integraltab = Monitor[Table[{[Theta], [Phi],integral[[Theta], [Phi],100, {-Infinity, Infinity}]}
, {[Theta], 0, Pi,Pi/10}, {[Phi], 0,  Pi, Pi/10}], {[Theta], [Phi]}]; // Timing
(*{74.7813, Null}*)

Without these modifictions the evaluation of 101x101 grid lasts around  1700s ( speedup  factor 23 ) !
Hope it helps!

Speeding up this NIntegrate

One Answer

Add your own answers!

Related Questions

Single expression containing the same pattern multiple times

How to get While to Return a Custom Value?

Manipulate mapping on lists

How to align a network graph to an image?

Vectorization of a function

Eliminate a variable within an integral

concept of tangent and residual in AceFem and automation

Contour plot a 2d array

Analytic continuation of the prime zeta function

Equations in text cell in wrong format

FEM: elastic cylinder under circumferential pressure

Control the display unit length of `PlotRange`

Most convenient way to set legend marker size?

Get file downloaded from click in WebSession?

center of mass of a 2D square

New Compiler Technology in V12

Applying Operators inside complex expressions to argument with Through

Why is the convergent product with GoldenRatio evaluated as ComplexInfinity?

Mathematica doesn’t simplify even providing the assumption

How to drag an overlay image over the reference image?

Ask a Question