Generic bandwidth and frequency resolution
Pre-release
Pre-release
-
Fixed an important bug affecting the D ≠ 1 runs.
-
Fixed some subtleties arising whenever the frequency array contains the origin:
- Modified the Fermi function to properly handle the beta = inf case, even when w = 0 (the exponential of inf*0 would give a NaN, which propagates dangerously downstream until breaking the dmft-loop when the norm of gloc_old-gloc is evaluated)
- Modified the Luttinger integral so to carefully exclude the nonanalytic pole at w = 0, which would totally break the sum-rule.
Everything carefully tested, we could be approaching a stable release.