Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Mathematica Mathematica: Speed up NIntegrate

  1. Jan 4, 2015 #1
    I have simplified my problem here, in my actual problem the matrix is much bigger which makes it impossible to find eigenvalues analytically. So, I chose standard BCS problem (2x2 matrix) to demonstrate my problem with Mathematica. But, with even this simple case, it takes forever to integrate. My code is the following:

    hamiltonian[k_,mu_,delta_]:={{k^2- mu, delta},{delta,-k^2+ mu}}
    eigens[k_,mu_,delta_]:=Eigenvalues[hamiltonian[k, mu, delta]]
    fermitotal[beta_,k_,mu_,delta_]:=Block[{ee = eigens[k, mu, delta]},1/(1+Exp[beta ee[[1]]])+1/(1+Exp[beta ee[[2]]])]
    nTotal[beta_,mu_,delta_]:=NIntegrate[fermitotal[beta, k, mu, delta] k^2,{k,0,20},Method->{"LocalAdaptive","SymbolicProcessing"->False},AccuracyGoal->4,PrecisionGoal->4,MinRecursion->10,MaxRecursion->300,WorkingPrecision->13]

    Any idea how to speed up the calculation?

    I have tried, using Compile for the hamiltonian matrix. Specifying the precision for the numbers,


    Did not help. This calculation is supposed to be fast, but I am completely stuck here and I don't know what causes this slowness. Any help would be very much appreciated. Thank you.
  2. jcsd
  3. Jan 4, 2015 #2


    User Avatar

    Staff: Mentor

    It's the WorkingPrecision->13 that is causing a slowdown. If you remove it, you'll get a solution quite fast.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook