Function QCM::wk_integral

Function Documentation

void QCM::wk_integral(int dim, function<void(Complex w, vector3D<double> &k, const int *nv, double I[])> f, vector<double> &Iv, const double accuracy, bool verb = false)

Performs an adaptive integral over frequencies and the Brillouin zone.

Performs an integral over frequencies and wavevectors Uses the Cubature library Actually computes $\int {dw\over\pi}_0^\infty \int {d^3k\over (2\pi)^3} f(iw)

Parameters:

  • dim – spatial dimension (1 to 3)

  • f – function to integrate (may be multi-component)

  • Iv – value of the integral (adds to previous value: must be properly initialized)

  • accuracy – required absolute accuracy of the integral